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ABSTRACT 

TOPICS IN INFLATIONARY COSMOLOGY AND ASTROPHYSICS 

by 

Matthew M. Glenz 

The University of Wisconsin-Milwaukee, 2008 
Under the Supervision of Distinguished Professor Leonard Parker 



We introduce a general way of modeling inflation in a framework that is indepen- 
dent of the exact nature of the inflationary potential. Because of the choice of our 
initial conditions and the continuity of the scale factor in its first two derivatives, we 
obtain non-divergent results without the need of any renormalization beyond what is 
required in Minkowski space. In particular, we assume asymptotically flat initial and 
final values of our scale factor that lead to an unambiguous measure of the number of 
particles created versus frequency. We find exact solutions to the evolution equation 
for inflaton perturbations when their effective mass is zero and approximate solutions 
when their effective mass is non-zero. We obtain results for the scale invariance of 
the inflaton spectrum and the size of density perturbations. Finally, we show that a 
substantial contribution to reheating occurs due to gravitational particle production 
during the exit from the inflationary stage of the universe. 

The second part of this dissertation deals with a post-Minkowski approximation to 

a binary point mass system with helical symmetry. Numerical solutions for particles 

of unequal masses are examined in detail for two types of Fokkcr actions, and these 

solutions are compared with predictions from the full theory of General Relativity and 

iii 



with post-Newtonian approximations. Analytic solutions are derived for the Extreme 
Mass Ratio case. 

The third part of this dissertation discusses the detection sensitivity of the Ice- 
Cube Neutrino Telescope for observing interactions involving TeV-scale black holes 
produced by an incoming high-energy cosmic neutrino colliding with a parton in the 
Antarctic ice of the south pole. Parton Distribution Functions and the black hole 
interaction cross section are computed numerically. Our computation shows that Ice- 
Cube could detect such black hole events at the 5-sigma level for a ten-dimensional 
Planck mass of 1.3 TeV. 
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Chapter 1 
Introduction 

This dissertation is an exploration of space on scales that are small (quantum fluctua- 
tions, TeV-scale black holes, vacuum particle creation); scales that are big (anisotropies 
in the Cosmic Microwave Background, seeding of large-scale structure, the Hubble 
radius); and scales that are in between (Extreme Mass Ratio binary black holes, 
temperatures associated with horizons. Innermost Stable Circular Orbits). The first 
part of this dissertation is an outgrowth of methods developed by my thesis advisor 
in the following works [H El El H] . These methods are applicable to the creation of 
quantized perturbations of the infiaton field, which is the topic we explore in Part I. 
The new results that appear in this dissertation are based primarily on the work of 
three papers. The first of these papers, "Study of the Spectrum of Infiaton Pertur- 
bations," examines an exact calculation of the evolution of quantum fluctuations and 
the subsequent particle creation in a model of the early expansion of the universe 
that is relevant to a wide range of inflationary potentials consistent with observations 
and that does not depend on renormalization in curved spacetime [5]. The second 
of these papers, "Circular solution of two unequal mass particles in Post- Minkowski 
approximation," computes numerically a set of solutions to a helically symmetric 
binary system of point masses in a particular approximation to General Relativity 
and presents analytical formulas for the limit that the mass of the lighter particle is 



negligible with respect to that of the more massive particle [6J. The third of these 
papers, "Black Holes at the IceCube neutrino telescope," calculates the experimental 
sensitivity for observing TeV-scale black holes produced by a gravitational interaction 
between a cosmic neutrino and an elementary particle within the atomic nuclei of ice 
molecules [7]. This dissertation is divided into three main parts corresponding to 
these three papers. 

In Part I, "New Aspects of Inflaton Fluctuations," we begin with a brief summary 
of early universe cosmology. Two of the most important cosmological theories of 
the twentieth century are the Big Bang theory and the theory of Inflation. The Big 
Bang theory supposes that our universe was once much smaller and much hotter that 
it is today. It explains the expansion of the universe, the presence of the Cosmic 
Microwave Background Radiation, and the primordial abundances of light elements. 
Cosmological Inflation supposes that the early universe underwent an extremely large 
increase in size in a very small amount of time. This explains why the density of our 
universe today is so close to the critical density that separates a universe that expands 
forever from one that eventually recollapses, it explains the near homogeneity and 
isotropy of the universe, and it explains why we don't observe magnetic monopoles. 
Most importantly of all, however, inflation explains the origins of those anisotrophies 
that do exist in our universe. Although a key ingredient of the Big Bang theory is a 
high energy density in the early universe and a correspondingly high temperature, the 
classical theory of inflation predicts an extreme cooling of the universe as it expands — 
much like the air in a piston cools as it expands to do work on its surroundings. 
We consider Reheating, and specifically the energy density of particles created by 
an expanding universe, as a means of preserving both theories without sacrificing 
any of their successes. We give a general overview of the amplification of quantum 
fluctuations into large-scale density perturbations during inflation, and we describe 
some of the ways of relating theoretical predictions to observations. We then list some 
of the observational findings of experiments. 



We continue with the details of the method we use to model inflation. Instead of 
specifying an inflationary potential, as is usually done, we specify directly the change 
in the scale factor, which is a measure of the size of the universe, versus time. We 
consider a scale factor that accommodates several parameters, but its most important 
features are that it asymptotically approaches a constant values at early times, that 
it approaches a different constant value at late times, and that its first two derivatives 
with respect to time are continuous. The asymptotically flat regions of our scale fac- 
tor allow us to associate our model with Minkowski spacetime at early and late times. 
Identification with a Minkowski vacuum at early times leads us to initial conditions 
that contain no infrared divergences, and comparison with a Minkowski spacetime at 
late times leads us to an unambiguous measure of the frequency-dependent density 
of particles created by the expansion of the universe. That our scale factor is contin- 
uous up to its second derivative with respect to time ensures we have no ultraviolet 
divergences, in addition to the prevention of infrared divergencies mentioned before. 
We choose for our scale factor a composite of three segments. The initial and final 
segments are each associated with a particular form of asymptotically flat scale factor 
with different choices of parameters. The middle segment of the scale factor, where 
most of the expansion takes place, is a region that grows exponentially with respect to 
proper time. Such an exponential growth is indicated by experimental observations. 
We solve for the matching conditions necessary to maintain the desired continuity 
of our composite scale factor. For each of our scale factor segments we have exact 
solution to the evolution equation for fluctuations of a massless, minimally-coupled 
scalar field. We also describe two different approximations to the case of a constant 
mass. We match up our solutions to the evolution equation at the interfaces between 
the segments of our composite scale factor, and at late times we are able to determine 
the particle production due to the expansion of the universe. Prom here we discuss 
the dispersion spectrum. We note the scale-invariance of the scalar index, provided 
the requirement is met that each mode be converted into a curvature perturbation at 



a time related to when it crosses the Hubble radius, and that all modes not be con- 
verted at once after the end of inflation. Using a hybrid combination of our method 
with the slow roll approximation, we describe a way of calculating the density pertur- 
bations produced by inflation. Finally, we show how Reheating, or a return to the hot 
Big Bang conditions after the end of inflation, can accompany inflation. We discuss 
possible consequences of Reheating and its relationship to constraints on predictions 
for exotic particles and high energy physics. 

In Part II, "Binary System of Compact Masses," we examine a post-Minkowski ap- 
proximation to a helically symmetric binary system of point masses. The helical sym- 
metry is maintained through the presence of half-advanced and half-retarded fields. 
The equations of motion are given for one of two Fokker actions — parametrization- 
invariant and affine — by Friedman and Uryu in p|, and from their results we calculate 
numerically the solutions in the case of unequal masses. We also derive analytical for- 
mulas for the Extreme Mass Ratio limit where the ratio of the smaller mass divided by 
the larger mass goes to zero. This limit would be applicable to the inspiral of a solar- 
mass black hole into a billion-solar-mass black hole, such as is predicted to exist at 
the centers of many galaxies. For both the numerical computations and the analytic 
equations, we plot three graphs: the angular momentum versus the velocity of the 
lighter particle, the unit energy of the lighter particle versus the angular momentum, 
and the unit angular momentum of the lighter particle versus the angular momentum. 
These plots are given for four mass ratios and for both types of Fokker action. For the 
parametrization-invariant case we include one of two different correction terms that 
generates solutions that agree with the first post-Newtonian approximation, and we 
demonstrate this in the Extreme Mass Ratio limit. We discuss the locations of Inner- 
most Stable Circular Orbits, and we compare the predictions of this post-Minkowski 
approximation with both those of the post-Newtonian approximation and those of 
the full theory of General Relativity. 

In Part III, "Production and Decay of Small Black Holes at the TeV-Scale," we 



investigate the possibility of using the IceCube Neutrino Telescope to detect TeV- 
scale black holes. In the physics of the Standard Model, it is not impossible that 
a cosmic neutrino could come close enough to an elementary particle in the cubic 
kilometer of ice in the IceCube experiment to form a black hole. Such interactions 
involving gravity, however, are so much less likely than interactions involving the weak 
force, that IceCube would never differentiate their signal from the background noise of 
weak-interaction event rates. Many theories of physics beyond the Standard Model, 
such as string theory, require additional dimensions of spacetime beyond the 3+1 
dimensions of our common experience. These additional dimensions might not have 
been noticed before if they were compactified, or curled up, with a simple example 
being the topology of a higher-dimensional torus. At the compactification scales, 
then, gravity would be much stronger than in a 3-|-l-dimensional theory, whereas 
at macroscopic scales gravity would appear to be much weaker than the strong and 
electroweak forces. In addition, if only gravitons propagated into the compactified 
dimensions, then the scale of compactification could be anything small enough not 
to conflict with observations. On distances smaller than this scale, gravity would 
grow stronger with decreasing separation faster than an inverse-square law would 
predict. If the strength of gravity were equal to the strength of the electromagnetic 
force around energies of roughly one TeV, or 10~^^ meters, the scale at which the 
electromagnetic and weak forces unify into the electroweak force, then gravity could 
be sufficiently strong that the IceCube detector could observe the production of TeV- 
scale black holes in the interactions between cosmic neutrinos and partons, which 
are the fundamental particles — both quarks and gluons — that are found within 
nucleons in atoms. For the high energies of interest for this experiment, the nucleons 
cannot be treated as single particles, which is why we treat them as collections of 
partons. At any moment, a parton can have an energy ranging from nothing to the 
entire rest mass energy of the nucleon, and parton distribution functions describe 
the probabilities of finding each parton with a given energy. We develop simple fits 



to a specific model of tlie parton distribution function, and witli tliis information 
we are able to numerically integrate an expression giving us the cross section for 
the gravitational interaction. The black holes formed by these interactions would 
decay almost immediately via Hawking radiation, or particles produced by the strong 
curvature of spacetime outside of black holes. The Cherenkov hght of these events 
could be measured by the photomultiplier tubes of IceCube, and signals could be 
picked out from the background event rate by searching for muon-daughter particles 
with less than 20% of the total energy, which is sufficiently unhkely in Standard Model 
physics that we would be able to discern TeV-scale black hole events from interactions 
through the weak force. We find that the IceCube detector could measure TeV-scale 
black holes at a statistically significant 5a excess for a 10-dimensional Planck scale 
of 1.3 TeV. 

The relationship between space at the smallest and largest scales is, perhaps, 
nowhere so evident as the inflation of quantum fluctuations from below the Planck 
length to sizes beyond our observable universe in what follows: Part I - New Aspects 
of Inflaton Fluctuations. 



Part I: 

New Aspects of 
Inflaton Fluctuations 



Chapter 2 

Inflationary Cosmology 

At the beginning of the twentieth century, most scientists beheved that the universe 
was infinite and eternal. Such a situation is not compatible with cosmology governed 
by the theory of General Relativity, which predicts that a static universe would be 
unstable to perturbations. From this it follows that our expanding universe started 
from a singularity of infinite density and temperature. This Big Bang theory of the 
universe successfully explains several observational phenomena. One of these is the 
expansion of the universe and Olber's paradox, which asks — if the universe is infi- 
nite, then why do we not observe stars in every direction; why do we see dark space 
between stars? With help from Hubble, Einstein and others came to realize that the 
universe is not only expanding, but it must also have a finite age. Thus, not all of 
the light from stars in the universe has had time to reach us, and for distant stars 
this light is redshifted by the expansion of the universe. Another question resolved 
by the Big Bang theory is that of the primordial abundances of the light elements: 
hydrogen, deuterium, tritium, helium-3, helium-4, and lithium. Stars convert hydro- 
gen to heavier elements through nuclear fusion, but the light elements are found in 
definite ratios in galactic dust thought never to have been part of any star. This 
is explained by looking back to the high temperatures and pressures of the universe 
when it was much more dense, shortly after the Big Bang. The universe was hotter 



than any star, and a series of calculations involving the thermal-equilibrium ratio of 
protons to neutrons, the ratio of baryons to photons, the half-life for a free neutron, 
and the cross section for neutrons to become bound in nuclei PHQ]; predicts ratios of 
primordial abundances of the light elements that agree very well with observations. 
A final success of the Big Bang theory is the explanation of the observed Cosmic 
Microwave Background Radiation (CMBR) at a temperature of approximately 2.7 
Kelvin. This was first discovered by Penzias and Wilson in 1965 while they were 
working at Bell Labs, and for this discovery they were awarded a Nobel Prize in 
1978. This background noise is the red-shifted relic of the early universe's radiation 
dominance. Although the Big Bang theory explained some questions about our uni- 
verse, Cosmological Inflation was necessary to explain other observed properties of 
our universe. 

Inflation was originally conceived to explain three primary phenomena. The first 
of these was the flatness problem. The density of our universe is surprisingly close 
to the critical density needed to close the universe, above which a closed universe 
would eventually re-collapse into a Big Crunch and below which an open universe 
would expand forever — neglecting acceleration caused by the presence of dark en- 
ergy. Surprisingly close, because unless our universe's density is precisely equal to 
the critical density — and there is no reason to assume it must be — the ratio between 
the two drifts rapidly away from 1 in a Big-Bang-only universe. Inflation solves this 
problem by very rapidly driving this ratio exceedingly close to 1 during a short period 
of enormous growth of the universe. The second argument for inflation is that all the 
CMBR is, to excellent approximation of within about one part in ten thousand, in 
thermal equilibrium. Just as the resolution to Olber's paradox involves light taking a 
finite time to reach the Earth, so does this present a problem for early- universe light, 
emanating from different directions, that is just now reaching us. In a Big-Bang-only 
model, widely separated regions of the currently observable universe weren't previ- 
ously in causal contact, and that they should be in thermal equilibrium now is a 
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mystery. This problem is resolved by explaining how the space in minute regions of 
our universe that were once in thermal contact expanded sufficiently rapidly during 
inflation to remove the different parts of the equilibrated sections to causally discon- 
nected parts of the universe: the space between points within equilibrated regions 
of the universe grew much faster than signals could travel across the distance be- 
tween those points. Thus, the CMBR reaching the Earth today, even from different 
directions, has come from regions of the universe that were previously in thermal 
equilibrium. The third issue that motivated inflation is the observed absence of mag- 
netic monopoles, which may have been created in the very early universe. Inflation 
resolves this by showing how monopoles could be inflated away with the expansion of 
space such that — unless monopoles were produced after inflation — on average there 
shouldn't be any monopole close enough to us to detect after inflation. 

Inflation has come up with an unforseen prediction that has since turned out 
to be more important than any of the historical justifications for its existence: the 
creation of fluctuations during inflation that lead to the anisotropics of our present- 
day universe. For NASA-COBE's (Cosmic Background Explorer) 1989 detection of 
these anisotropics in the CMBR, Mather and Smoot were awarded a Nobel Prize in 
2006. In the most widely used models of inflation, this expansion is driven by the 
inflaton field, which is a scalar quantum field, and the perturbations of the inflaton 
field seed galaxy formation and are responsible for large-scale structure of our universe 
today. 

2.1 Cosmology in General Relativity 

In units of c = h = 1 Einstein's equation is [HI [12] 



Gab = Rab — TT-Rfl'ab 
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On large enough scales, our universe appears to be of a fairly uniform density in all 
directions. If the Earth is not in a privileged position in the universe, this implies 
that the universe is homogeneous and isotropic. Following the example of [T2l [T3]. if 
we assume no distinction between the spatial directions, we can write the Friedmann- 
Robertson-Walker-Lemaitre (FRWL) metric as 



ds^ = -dt^ + a{tf 



dr"^ 



1 — fcr^ 



+ r\d9^ + sin^ 9 dcj)"^) 



(2.2) 



where a{t) is the scale factor that relates the chosen coordinate scale to the proper 
time t, and the variable k describes the topology of the universe: k > corresponds to 
positive curvature (closed universe), k = corresponds to zero intrinsic curvature (flat 
universe), and A; < corresponds to negative curvature (hyperbolic, open universe). 
We then have 




a(t)V2 



gab 



-1 








1— fcr^ 

-^^^2-2 











a{tyr^ sin^ 6 



(2.3) 
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ab 



-1 













1— fcr^ g 

a(t)-V-2 



a(t)2 











(2.4) 



a{t)-^r-^ sin"^ e 

In this section, only, we will not use the Einstein summation convention. In the basis 
of {t, r, 6, 0}, the Christoffel symbols are given by 



ab 



^g"'^ {dagbd + dbOad - ddgab) 



(2.5) 



where VaV" = daV^ + T^kV'' and VaW^ = daW^ - T\jVb, with da the covariant 
derivative operator of the fiat metric [llj. For the metric given by Eq. (12. 2p . we see 
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that g'^'^ = 5fg'^'^ and Qcd = ^tdco where 5^ is the Kronecker delta, so we have 



(2.6) 



In the set of coordinates defined by {t, r, 6, 0}, we consider the four cases of a = 6 = c, 
a = h^c^a^h = c, and a ^ c (each of the indices is different in this last case) 
to get 



a = h = c 

a = h ^ c 

a =^ h = c 

a ^ h ^ c 



r 

r^afe = 



cc o 

ca 2^ ^aQm 



(2.7) 
(2.8) 
(2.9) 
(2.10) 



The non-zero derivatives are dtgrr, dtgee, dtg^^, drgrr, drgee, <9rfi'00, and deg^^. Thus, 
the non- vanishing Christoffel symbols are r*^^, T^ee, r*(/)</<5 T'^rr, T^ee, T'^t/x/), T'^rt = T'^ir, 

r^00; r^e* = r^te, r^er = r^re, r'^</>t = 1^*^*0, r'^</>r = r'^r<ji, and r'^^g = r'^e^. 

When we write the Ricci tensor as [T2l 



we find, using an underline to indicate terms that cancel, using an overline to indicate 
terms to be consolidated, and using a = a{t), a = da/dt, and d = dd/dt, that 



Ru = -dt {rrt + T%t + T%t) - {T\tr\t + r\tr\t + r'^^tT^t) 



a a 



a a 



a a 



a / \a I \a 



(2,12) 



+ 



( r rr^ rr ~1~ T ^^'^-T 0j- -\~ T (^r-T'^f^r- ~1~ ^i-T ^-^-^ 



aa + a 
1 — /cr^ 
kr 



+ 



+ 



aa a a a 



1 1 

1 — /cr^ V r 



1 — /cr^ V a 



- + 2- + 2- 
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(2.13) 



= (rMad + d^}-{l-3^r^} + {l + ^}) 



+ 



a a a 

vaa I 1 — 

a a a 



1 1 
- + - 
r r 



2 2 

a r 



a 



a 



+ 2— + 2- 



cos 



sin^^ 



(2.14) 



Rd 



= (a^rV^ + 9,r^^ + d0r%^) + [r*^^ (r,* + r\t + r^) + (rv^) 

= ({ad + d^} sin^ ^ - {l - 3kr^} sin^ 9 + {sin^ ^ - cos^ 6}) 
+ 



» I ' J I 

adr sin 9 [ — I 

a a a 



{r - kr^} sin^ 9 



kr 



1 1 



(2.15) 



The Ricci Scalar Curvature is 



R = g'^Rab 

ab 
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g^^Rtt + g" Rrr + g Ree + 



, a cP' k 
6( - + - + - 



(2.16) 



The most general stress tensor associated with homogeneity and isotropy is that of a 
perfect fluid [12], given by 



Tab = pUaUb + P {gab + UaUb 



(2.17) 



where p is the energy-density, P is the pressure, and in these coordinates f/" 
(—1, 0, 0, 0) is the four-velocity of a comoving observer, and 



Ub = Y,gabU\ 



(2.18) 



The time-time components of the Einstein Equation, Eq. (12. ID . give us the Friedmann 
equation: 



Gtt 



a 1 

-3 

a 2 



'da?' k 
6 ( - + — + — 



-1) = 3- + 3-^ = SttGp, 
a 



(2.19) 



or, 



mr = ■ 

where the Hubble constant is defined by 



iTiG k 



3 



(2.20) 



d a{t)/dt 
ait) ■ 



(2.21) 



Any same space-space components of the Einstein equation, for which we will use r-r, 
give us the Raychaudhuri equation: 



(j„ — g„. — h 2— + 2— 



^ , d k 
6 ( - + ^ + ^ 



g„ = ^iiGPg„, (2.22) 
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or, 



dak ^ 
2- + — + — = -MGP, 



(2.23) 



which, when we use H = H{t) and H = dH/dt = a d — a d , can be written 



2H + 3H^ + — = -8nGP, 



(2.24) 



which we rewrite, using Eq. (12.201) . as either 



H = -AnG{p + P) + 4 



(2.25) 



or as the Raychaudhuri equation, which is 



(2.26) 



We get the continuity equation by taking the time derivative of Eq. fl2.20p and then 
inserting Eq. ( 12.25^ to find 



8nG 



p = 2HH = 2H 



-47rG(p + P) + -^ 



(2.27) 



which becomes 



p = -3H{p + P) + 



3H k 



(2.28) 



In a fiat universe, where k/a^ can be neglected and the metric can be written as 
ds"^ = —dt^ + aity^dx^ + dy"^ + dz^), the continuity equation becomes 



p = -m{p + p). 



(2.29) 



A simpler way of deriving this equation would be to use conservation of energy in a 
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comoving reference frame to show, in units where E = mc^ = m, that 

/E\ M P 
\y) = -Y'^y - ydV^ (2.30) 

where M = pV and V (x a^. If there were no pressure, as is the case for what is re- 
ferred to as dust, then in the coordinates {t, x, y, z} this would reduce to conservation 
of a density current: 

a 
a 

= U'dtp + p{U'T\t + U'T%t + U'T%t) 

= -dtp-3Hp. (2.31) 

For dust, which is the term for matter that satisfies P = 0, such as cold dark matter 
and — to good approximation — galaxies, we can solve the differential equation 



^ = -3-, (2.32) 



by integrating both sides with respect to time to get 

Inp oc — 3 In a, (2.33) 

or 

p oc a-l (2.34) 
We combine this with Eq. (I2.20p to get 



(2.35) 
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which leads to 

a oc fl-^/^ (2.36) 

and (with = 0) 

adust(t) oc e'\ (2.37) 

We refer to this as a matter-dominated universe. For the case of a radiation-dominated 
universe, where radiation obeys the equation of state 

P = Ip, (2.38) 



we would have 



and (with k = 0) 



p oc a-^ (2.39) 

a oc a-\ (2.40) 

'^radiation 

(t) OC t^'^. (2.41) 



In the next section we will show that a slowly- changing scalar field displaced from its 
minimum potential energy obeys the equation of state 

P ~ -p, (2.42) 

for which we have from Eq. (12.251) 

i^inflation ^ 0. (2.43) 

We discuss inflation in more detail in the next section, but first we mention that 
with a time- invariant Hubble constant, we would have (in a flat universe) a de Sitter 
metric given by 

ds'' = -dt" + e2^*(da;2 + dy^ + dz^). (2.44) 
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Whether A; = in Eq. (12.201) . or not, we may define a critical density that would 
produce an equivalent Hubble constant if k were 0. This we define as 



(2.45) 



_ 3^ 

We define the density parameter as 

where a in a fiat universe {k = 0), we would have Q = 1. One of the primary 
motivations for infiation was reconciling observations that in our universe i7 ~ 1, 
when there was no reason to expect that it necessarily would be. In fact, in either a 
radiation- or matter-dominated universe (for both H oc when ~ 0), we should 
expect 

where k (x k. The Big Bang theory predicts — based on the presence of the ap- 
proximately 2.7K CMBR and the relationship between the current matter density 
and Hubble constant — that our universe was radiation-dominated until it was about 
300,000 years old and has been roughly matter-dominated (neglecting any recent ac- 
celeration of the universe due to dark energy) since then. Thus, Q in our universe 
should diverge rapidly from 1, unless the value of k was very nearly zero at early 
times in our universe. One mechanism for driving Q close to 1 is infiation. When 
a(t) = e^* and H = constant, we have 



Inflation very rapidly drives the value of towards 1. With enough inflation, an 
initial value of Q that may have differed from 1 by orders of magnitude, could have 
been driven close enough to 1 that it would still be approximately equal to 1 in our 
universe today. For the rest of this dissertation we will assume that the universe is 
flat, in the sense that we will take the curvature constant k to be zero. From now on 
we will not make use of this variable and will reserve k for other quantities, namely 
the Fourier mode-number. 



2.2 Inflation 

For the rest of this dissertation, we will adopt the Einstein summation convention. 
The Lagrangian density of a scalar field with metric signature of +2 is [T5[ [161 [13 



J^ = l\9\'^' i-9^'dM<P - n?<\? - iR<\?), (2.50) 



where g = det{gab)- A massless (m = 0), uncoupled = 0) field with a (^-dependent 
potential, where the potential may incorporate a non-zero scalar field mass, becomes 



^ = Igl'^'g^'da^pd^^-lgl'^'vi^p). (2.51) 



The origin of this potential depends on the various models being considered, but 
the main prerequisites are that (p initially be displaced from the true minimum of 
the potential, and that some portion of the slope of the potential must be relatively 
flat with respect to changes in during the slow roll approximation, for which see 
Sec. 12.3.11 If we were to retain the Ricci curvature scalar in Eq. (12.501) . then the 
variation of the action would lead to the Einstein Eq. (12.11) in the calculation below 
[I6][ni pp. 491-505]. The action is [l5j 



S = I d^x' C = / d^x' 



l\g\'^'i~g'''''da'(l)dk'(j)-2V) 



(2.52) 
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and the stress-energy tensor is 



ab 



2 6S 



(2.53) 



Using the identities [16 



leads to 



l\g\'/' g''%ab, 



ab 



2 6jd'x' [l\gix')\'^'i-g'^'''ix')da,<pd,,(j)-2V)] 
\gix)\'^^ Sg^K^) 

6jd'x'[\g{x')\'/' {-g'^'''{x')d,4db4 - 2V)] 

1/2 Xnabt 



\g{x)\ ' 6g''''{x 



, ,\g{x')\'^'Sg-'''ix') 



d X 



\g{x)\^'^ Sg^'ix) 



(2.54) 
(2.55) 



ga'b'{x'){--d''(t)dc(t) -V) + da'(t)db'(t) 



(2.56) 



Finally, using the delta function identity [16 



(2,57) 



the stress tensor is 



Tab = gab[--d"(t^d,(l>~V\+da 



(2.5^ 



and 



T\ = g\[-:^d^ 



V\+d'' 



(2.59) 
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The spatial slicing and coordinate threading of time is chosen such that = In 
absence of perturbations, space-time is homogeneous and isotropic: 



(2.60) 



where a dot represents derivatives with respect to time. Because of homogeneity and 
isotropy, the stress tensor is described by a perfect fluid, 



Ta 
h 



-p 














p 














p 














p 



(2.61) 



where p is the energy density and P is the pressure. It is now possible to solve for 
the energy density and pressure: the energy density is equal to minus the time-time 
component of the stress tensor; and the pressure is equal to any of the three diagonal 
space-space components of the stress tensor [10] . 



n - -T° 
P — -^0 



1 V 



V 



(2.62) 



P = T\ = T\ = T\ = -<p'~V{<P). 



The Friedmann equation, 



and the continuity equation. 



SvrG 



p = -'iH{p + P), 



become 



8nG (\ ■ 



3 V2 



-0^ + V(0) 



(2.63) 



(2.64) 



(2.65) 



(2.66) 
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and 



00 + ^0) 



j,2 



(2.67) 



+ 



d(j)/dt 



-3H(j), 



(2.68) 



0+ V' 



(2.69) 



where a dot represents a derivative with respect to time and a prime represents a 
derivative with respect to (p. The curvature term in the Friedmann equation is here 
set to zero. Whether or not this is precisely the case, soon after inflation begins the 
curvature of the universe will become negligible. 



Well after inflation has begun, the scalar field can be treated as a homogeneous, 
isotropic classical field with the fluctuations consisting of quantum perturbations. In- 
flation smooths out all other perturbations to the point that quantum fluctuations 
are all that remain. For models of inflation driven by a single scalar fleld, pertur- 
bations can be expressed as time-dependent, location-dependent fluctuations on a 
homogeneous, time-dependent background: 



2.3 Quantum Fluctuations of a Scalar Field 




(2.70) 



The Euler-Lagrange equation. 



dC 



0, 



(2.71) 




with Eqs. (12.501) and fl2.5ip . becomes 




(2.72) 
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or 

1 



da{V^g''%4>)-V'i<l)) = 0, (2.73) 
which is equivalent to [IHl p. 38] [171 p. 542] 

□0 - V'{(j)) = 0. (2.74) 

If we perturb this with Eq. f l2.70p . then we get 

□(0 + 50) - 1/'(0 + (50) = 0. (2.75) 

To first order in 50, we write this as 

□0 + 050 - [V'{(f)) + 6(f)V"{(f))] = 0, (2.76) 

and we then use Eq. f l2.74p to show 

□50 - 50y"(0) = 0. (2.77) 

We can see from Eq. 12.501) that for a free field we may make the association 

V"{(f)) =m^ + ^R, (2.78) 

where m is the scalar mass, is the coupling constant, and R is the Ricci scalar 
curvature. 

The perturbation of Eq. fl2.70p expanded in terms of creation and annihilation 
operators is [16] 

50 = (volume)-^/2 ^[ag^fc(t)e*^-"^ + H.C.]. (2.79) 

k 

where 

volume = [La(t)]^ (2.80) 
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which is the physical length found from multiplying the coordinate length times the 
scale factor. The time dependent part of the fluctuations is ipk, where 



'4'k = a{t) 2^fc, 



(2.81) 



and 



(2.82) 



The solution thus far has periodic boundary conditions, but in the limit that L oo, 
a volume even as large as the observable universe will not be affected by this choice 
of boundary conditions. Combining the metric 



ds^ = -dr + a{ty{dx^ + dy^ + dz^) 



(2.83) 



where 



and 



9ab 



-10 

a(t)2 

a(t)2 

a{ty 

-10 

a(t)-2 

a(t)-2 

a(t)-2 



V\Fm¥¥W¥m\ = aitf 



(2.84) 



(2.85) 



(2.86) 



with the massless, uncoupled scalar field equation [15 



n5(f> - 50F"(0) = _-5,(|^7r/V5fc50) - 50v^"(0) = 0, 



(2.87) 
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yields 



aity^'dt [a{t)'{-l)dt5(P[ + a{t)~'d, [a{t)' {a{t)-') & 
d^6(j) + 3H{t)dt6(f) - a{ty^d,d'6(f) + 6(f)V"{(j)), 



- 6<PV"{<P) 
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where 



H 



da/dt 



(2.89) 



With the spatial dependence given by Eq. fl2.79p . the evolution equation for mode-k 
becomes 



h2 

d^Scj) + 3H{t)dt5(l) + -rry^H + (50V^"(0) = 0. 

(X\T/ ] 



(2.90) 



Using the scale factor associated with the de Sitter universe given by Eq. (12.441) . 



a = e 



Ht 



(2.91) 



and assuming a constant value of V"{(j)) to simplify the calculation, leads to an evo- 
lution equation for mode-A; of 



h2 



(2.92) 



Combining this with Eq. (I2.79P leads to 



9 



+3H 



+ 



e-'^'^'gk 



e-'^'^'d^g, + ee-^'^'g, - ^-H'e-"^'''g, + V"e 



9 



(2.93) 



Using the change of variables, 



-m 



(2.94) 



which is k times the conformal time, we then have 
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<9t = — — = ke-^^du, 
dt du 



(2.95) 



and 



ke-^'duke-'^'du = ke-'^'du [-Hudu 



-Ht 



Hti 



-kHe-^'du - ukHe-'^'dl 



-Ht^2 



(2.96) 



so the evolution equation Eq. (12.931) for mode-fc in terms of u is 



rll„-%Ht. 



k 



k 



k' 



-2Ht 



9 



V" 



u d^gk + udugk + 



u 



9 _ 

4 IPJ 



9k- 



(2.97) 



Eq. (12.971) is Bessel's equation. The most general solution for a given fc-component, 
gk, is tlBj 



g,{t) = \,/^LH''l^{u) + C2H' ^„ 



r(i) 



r(2) 



[U 



We then have 



(2.98) 



(1) 



/9_ V" 
4 JF2 



(m) + C2H 



(2) 



/9_ V" 
4 JF^ 



(2.99) 



but for the k = mode of the massless, minimally coupled case in a purely de 
Sitter universe, a universe that has an infinite history and future that is at all times 
described by the metric of Eq. (12.441) . see also Refs. [IHl [20] . 

For sufficiently large fc-modes the solution should be asymptotically insensitive to 
the de Sitter curvature, as this corresponds to very small wavelengths. On a very 
small scale that locally appears nearly fiat, the curvature becomes negligible. For 
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these large fc-modes, the solution we expect — due to the rapid attenuation of matter 
and radiation in a de Sitter universe — is that of the positive frequency WKB vacuum 
solution [161 



Mi) ^ ^~ij^,{t')dt' ^ (2.100) 



where the frequency is 



h2 

Mt)^\l^ + ^'- (2.101) 



See also Sec. 13. 2[ To match our constants, ci and C2, when A; — > oo, we use the large 
argument expansion of the Hankel functions [21] 



hI^\z) ~ v/27Me-*(^-^--i-), (2.102) 
which means that, to within a phase, 

Mt) = , {cie^^ + c^e-™} . (2.103) 

To match to the positive-frequency, vacuum solution given by Eq. (12.1001) we choose 

iiafTa 



lim Ci(fc) ~ 0, 

fc— >oo 

lim C2{k) ~ 1. (2.104) 



The de Sitter metric and the physical volume are symmetric under the transforma- 
tion [16] 

t ^ t + to and f ^ e"^*«f . (2.105) 
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The Killing vector generating this isometry, 



e° = 1, C = -Hx\ (2.106) 

corresponds to conservation of energy. Since the vacuum fluctuations can be expected 
to share this symmetry of space-time, provided — as will be explained in Sec. 13.4.11 — 
there is an infinite expansion and the universe is de Sitter in the infinite past and 
infinite future, the variable u is thus invariant under 

t ^ t + to and ^ fce^*°. (2.107) 



Then, with k' = ke"^\ 



^k'{t + to)=Mt) (2.108) 



requires 



c^{k') = ci{k) and C2{k') = C2{k). (2.109) 
Thus, because to is arbitrary, we have 



4 h2 



We note for future reference that changing the sign of the argument in Eq. (12.981) 
also yields a linearly independent solution to Eq. (12.971) under the transformation 
u ^ u = —u, because the Hankel functions of the first and second kind form an 
orthogonal and complete set. The coefficients ci{k) and C2{k) will, in general, change 
under the transformation u ^ u, but the procedure outlined above for finding these 
coefficients in the /c — oo limit, leads to ci{k) = —i and C2{k) = 0. A simpler way 
of seeing this, once we have Eq. (12.1101) . is to change the sign of H. Although we 
will later take H to be real and positive, we have not yet made this assumption, so 
changing the sign of H should leave Eq. (12.1101) intact in the fiat-space limit of /c oo, 
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where again a mode should not see the curvature of space. Using Eq. (12.1021) . we see 
that this large argument limit of the Hankel functions takes — to within a phase — 
H^^\z) ^ -tHi^\-z). 

2.3.1 Relation to Observations 

In this section we will focus on defining the slow roll approximation, the slow roll 
parameters, the number of e-folds, the curvature perturbation, the spectrum of cur- 
vature perturbations, and the spectral index. 

In the slow roll approximation [MlESlESlEZlEHlESlEniEIlES] 

0' < V{(P) (2.111) 

and 

101 < IV^'I- (2.112) 
This means Eqs. (I2.66P and (12.691) become 

~ ^V(0) (2.113) 



and 

V 



(2.114) 



3i/ 

These conditions ensure that P ~ — p, which is the property of a space-time domi- 
nated by a cosmological constant, or de Sitter space; and that the kinetic term does 
not grow appreciably since the potential is assumed to be flat and H is large. During 
inflation, the slow roll parameters must satisfy [33] 



e ^ 1 and \ri\ -C 1, 



(2.115) 
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where the slow roll parameters are defined by 



;2.116) 



1 fV'V H 



^ = ^ ( ^ ) ' (2-117) 



Using the slow-roll equations (12.1131) and (12.1141) . we can express the number of 
e-folds of infiation as [3 



iVe^ln 



Ct(ifinal) 




.'^(iinitial) . 





"tfinal /"(/"initial 'U' ( A)) 

Hdtc^MG / T7T-r^d(P. (2.118) 



We define a mode to be crossing the Hubble radius when the mode's wavelength, 
a{t)/k, is the same size as the Hubble radius, H~^, which would be the horizon size 
in a purely de Sitter universe. During infiation, when the scale factor is growing 
exponentially and k and H are both constant, a mode exits the Hubble radius when 
k/[a{t) H] = 1. After infiation, when the scale factor is given by either a radiation- 
dominated a{t) oc t^^"^ growth or by a matter-dominated a(t) oc t^^^ growth, where for 
both cases H oc t~^, then k/[a{t)H[t)] = 1 defines the time when a mode re-enters 
the Hubble radius. 

We can apply the small argument limit of the Hankel functions [211 Eq. 9.1.9], 



\Hi'H,)f^m.)\'^(^) (ikl)"", (2.119) 

when the real part of the parameter v is positive and non-zero, to Eq. (12.1101) . to get 

In the massless, minimally-coupled case, v = 3/2, and we find 



|2 



k\ 



^a(tr^(^X('Y. (2.121) 

4:H ^ ' \ n J \2a{t)Hj ^ ' 
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Late enough into inflation for a given mode to be well outside the Hubble radius, we 
then have 

^ (2.122) 

which is approximately half the value of \ipk\^ at the time it exits the Hubble radius — 
see Sec. 13.51 Although this perturbation of the infiaton field is not a gauge-invariant 
quantity, there is a gauge-invariant quantity, a curvature perturbation that we call 
TZki that is approximately conserved outside of the Hubble radius, and we can use it 
to relate the infiaton fluctuations to density perturbations at the time of re-entry as 
follows: [9l[T0l[2l[27l[29l[3l[35l[36l[37] 

~ 7^fe,exit - '^fc,re-entry OC = — , (2.123) 

P 

where for re-entry into a matter-dominated universe 5k — §7?.^, and for re-entry into a 
radiation-dominated universe 5k — ^Tlk- The value of 50^ is usually taken (neglecting 
the coordinate length L) to be the unrenormalized value H'^/k^ obtained at the time 
of exiting the Hubble radius. The justification for using an unrenormalized value of 
5(j)k, when it is well known that the Bunch-Davies state given by Eq. (12.1101) leads to a 
divergent 5(j) when summed over all modes, is usually given as implicit large and small 
cutoff frequencies. It is often assumed that the infrared and ultraviolet divergences 
come from infrared and ultraviolet frequencies that do not affect the treatment of 
modes exiting the Hubble radius during inflation. Parker [38], however, has shown 
that the divergences affect every mode, and that neglecting a proper renormalization 
drastically alters the results that are obtained. 

We use the definition of a spectrum given by Liddle and Lyth 



Vf{k)^(^-^yArre{\fkf). (2.124) 
Thus, under the standard assumption that it is not necessary to renormalize the 
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inflaton fluctuations as they are exiting the Hubble radius, we could show 



from the super-Hubble radius behavior given by Eq. f l2.122p . The renormalization of 
[38], however, changes this: the renormalized spectrum of inflaton perturbations, at 
the time of exiting the Hubble radius when k/[a{t) H] = 1, is 

^'^-[2^) U' " ^ (rnj, + J' (2.126) 

where tjih = m/H and n = 1^9/4 — mn'^. The renormalized inflaton fluctuation 



depends critically on the mass and when the magnitude of the fluctuation is evaluated. 

2 

In the massless case, 



= 4/7r, and the renormalized fluctuation is precisely 
zero. Well outside the horizon, the renormalized Vs^ also goes to zero, but this is 
perhaps not a problem, as TZk is the conserved quantity, not 50^, and the value ofT^^ 
given by Eq (12.1231) is typically evaluated at the time a mode crosses the Hubble 
radius. Thus, renormalization has the potential to greatly alter the character of the 
spectrum of perturbations. 

The scalar spectral index, n^, is a measure of how the magnitude of density per- 
turbations changes with scale. A value of n<j = 1 indicates scale-invariance. A value 
less than one is called a red-tilted spectrum, and a value greater than one is called a 
blue-tilted spectrum. It is defined as 

(2.127) 

where the value of Us is given for a specific value of fc, called the pivot value, which is 
normally either of A; = 0.05 Mpc~^ [39] or k = 0.002 Mpc~^ [lO], relative to the value of 
the scale factor fixed to be such that a(tnow) = 1- There is little running, or change in 
ns{k) with changing scales, so the choice of fcpivot is somewhat arbitrary. We can relate 
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the scalar spectral index to the slow roll parameters given in Eqs. ( 12.1161) and (12.1171) 
[311 SI]- Because the curvature perturbations are evaluated at the time of Hubble 
radius crossing, when k = a(t)H ~ He^^, we see that with a nearly constant value 
of H during inflation d\nk = d[\n{H) + Ht] ~ H dt. This leads to, with Eq. (12.1141) 
rewritten as dt = —3H/V' d(j), 

d V d 1 V d , , 

(2.128) 



d\nk 3H^d(l) SttG V 

Again using Eq. (12.1141) . the spectrum of curvature perturbations given by Eq. (12.1251) 
becomes 

p..(l!V(iLV. 



V J \2tt 



With Eq. (12.1131) . this becomes 



where the observed value of Vn is typically listed for the specific value of = 
0.002 Mpc^^, which is different from the value of k used with the scalar spectral index 
|39] . In [ID], the pivot scale for the spectrum of curvature perturbations is chosen to 
be = 0.02Mpc"\ as this is a scale that puts tighter constraints on the magnitude of 
the curvature perturbation spectrum for a wider array of model assumptions. Within 
the assumptions of various models, there is still a relatively scale-invariant spectrum 
of curvature perturbations. 

With Eq. (1271281) . Liddle and Lyth find 

dlnVn 1 V d ^ fiSnGf 

m ' 



dink 87rGVd^ \ Utt"^ V'^ 

1 V d 



8nG V 



(31nV -21n V) 
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167rG \ V J 87iG\V ) 
~ -6e + 2r/, (2.131) 

where the slow roll parameters are given by Eqs. fl2.116p and (12.1171) . Thus, 

Us-l = -Qe + 2r]. (2.132) 

See also the end of Sec. 13.51 for a slightly different derivation. 

Finally, we note that when we define the mass by = (PV/dcjP', we find that 

mn = m/H = ^m^jE'^ = y/3V"/{8nGV) = (2.133) 

and thus the effective infiaton mass is related to the Hubble constant during inflation 
through the slow roll parameter 7]. 

2.3.2 Findings of WMAP and SDSS Experiments 

The Five- Year Wilkinson Microwave Anisotropy Probe (WMAP) data measures a 
scalar spectral index of 71^(0. 002/Mpc) ~ 0.96 [12]. The Sloan Digital Sky Survey 
(SDSS) measures a scalar spectral index of ?T,s(0.05/Mpc) ~ 0.95 ^3J. Because the 
WMAP experiment measures fluctuations in the CMBR, while the SDSS observes the 
locations of galaxies and large-scale structure in our universe, there is good, indepen- 
dent accord for the red-tilted spectral index measured by these different approaches. 
The Five- Year WMAP data finds a curvature perturbation spectrum of 

P7e(0.002/Mpc) ~ 2.4 x 10"^ (2.134) 

What follows in this section, where we apply these observations to two particular 
models, is based upon work done by |44]. The first model we consider, the quadratic 
chaotic inflationary potential [45t H6] . is in good agreement with the Three- Year 
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WMAP data [47j. The second model, a type of Coleman- Weinberg model 
in good agreement with the Five- Year WMAP data [IHl HQ]. 

The quadratic chaotic inflationary potential is given by [IHl US] 



IS 



(2.135) 



In Fig. l2.H we plot the potential given by Eq. (I2.135P versus 0. In chaotic inflation, the 
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Figure 2.1: Quadratic Chaotic Potential. 



value of is initially perturbed away from the minimum and rolls slowly — provided 
the slope of the potential is sufficiently gradual — down the potential to the minimum 
at = 0. To be contrasted with chaotic inflation is new inflation, in which begins 
near the maximum value of the potential located at = and rolls slowly to a 
minimum of the potential [50] . The Coleman- Weinberg potential, which was actually 
one of the earlier models considered for an inflationary potential that did not involve 
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tunneling through a potential barrier and its associated problems with bubbles of 
inflation not coalescing, is an example of new inflation. 

The one-loop Coleman- Weinberg potential is given in the zero-temperature limit 
by [291148] 



ln(</)Va^^ 



(2.136) 



In Fig. 12. 2[ we plot a dimensionless potential V{(j))/{Ba ) versus the dimensionless 
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Figure 2.2: Coleman- Weinberg Potential. 
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parameter 0/cr. In the low temperature limit, the stable minima of the potential are 
located at = ic. At the beginning of inflation ~ 0, where the slow roll conditions 
are satisfied, and rolls to either of two (in the low temperature limit) stable minima. 
Classically, inflation is a period of super-cooling, so the low-temperature limit should 
be justified, but see also Sec. 13.71 

For the quadratic chaotic inflationary potential, the slow roll parameters of Eqs. (12.1161) 
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and (12.1171) are equal to each other, and we have 

From Eq. (12.1321) and the Five- Year WMAP spectral index of n,, — 1 ~ —0.04, we 
find 

- 1 = -6e + 2r7 = -4r7 ^ -0.04, (2.138) 

or 

e = 7^-0.01. (2.139) 

From Eq. (12.1331) . we have 

m 



V(r03~0.2. (2.140) 
H 



From Eqs. (ITOTD and (127[39D . 



47rG0 



or 



^ 0.01, (2.141) 



^2.142) 



V0.047r 

where 0cmb corresponds roughly to that range of at which the modes observed by 
WMAP were exiting the Hubble radius during infiation. Using Eq. (I2.118p . we find 
the number of e-folds before the end of infiation at which these modes were exiting 
the Hubble radius: 

/</>cmb Y r4>cmh 1 d)^ 1 

—dcj) - 8nG J -(pdcp- SttG-^ - ^ ~ 50. (2.143) 

For the value of itlh — 0.2 given by Eq. (12.1401) . the renormalized spectrum of infiaton 
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fluctuations given by Eq. fl2.126p is 

i7\V 2.237\ f hV 



Vs^ - (^1.968 - - (0.019) - 0m0A7H'. (2.144) 

Using the relation given in Eq. (12.1251) and the slow roll approximation given in 
Eq. 02.1141) . we have 

Vn = (S-Yvs^ ^ (^)\m047H' ^ ^0.0042i7^ ^ (2.145) 



then, as a rough estimate of the general order of magnitude, we use 0cmb to get 

A.7H^ _f H 



We can equate this with the amplitude of the spectrum found in the Five- Year WMAP 
data to write 

H 



37(^ZT7j) -2•4xl0-^ (2.147) 



and 

H 



8 X 10^^ (2.148) 



Q-l/2 

Using the Planck scale, G"^/^ ~ 1.22 x 10^^ GeV, finally we have 

H ^7x 10^^ GeV, (2.149) 

which can be seen as an upper limit on H near the beginning of inflation, around the 
time the modes observed by WMAP were exiting the Hubble radius; as (p rolls down 
the potential towards zero, the size of H decreases. 

For the Coleman- Weinberg potential given by Eq. (I2.136p . we have 



6"^ 

V = 450^ In (2.150) 



9 ' 
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,3 a 

The slow roll parameters are 



4503 1^ <^ 



V" = 12S02(^ + ln^). (2.151) 



167rG \\B<7^ + [ln(0Va2) - \] 



167r(T2 \^ i + H [ln(r2) - 1] 
12502 (^l + lng 



' SvrG 1 \Ba^ + 50^ [ln(02/a2) - \] 

(G-i/2)2 / I2r2 (| + lnr2) 
87r(T2 l^i + r4 [ln(r2)-i] 

where r = 0/cr. We assume the values given by P, p. 292] of 



(2.152) 



(2.153) 



a ~ 2 X 10^^ GeV, 

B ~ 10"l (2.154) 

With those values and G~^/2 ~ 1.22 x 10^^ GeV, taking r ^ 1 we find 

e ~ (7.4 X 10^) 64r'^ (Inr^)^ (2.155) 
ri ~ (1.5 X 10*^) 24r2lnr^ (2.156) 

and we find in the limit <^ a, that e <^ rj. Using the WMAP value of 0.96 for the 
spectral index, this leads to — 6e + 2// ^ 2// ~ —0.04, or 

T] - -0.02. (2.157) 

Then we have ttih^ ~ —0.06, or 



niH — 0.2452 ~ i/A. 



(2.158) 
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An imaginary physical mass could lead to tachyonic behavior J5T\, however in this 
case, recall we are dealing with an effective mass. To find r, which we assume to be 
much less than one, we combine Eqs. fl2.156p and (12.1571) to get 



r ~ ±5.3 X 10~^ (2.159) 



With Eq. (I2T18D . we have 



'V. - I (^iJ^M I *• - 64. (2,60) 

Finally, Eqs. fl2TT4l) . fl2J25l) . (^Jm . fl27[34l) . and fl2:i58l) lead us to 

2.4 X 10-9 ^ (^] (0.012) ^ I -— ^ I (0.012) . (2.161) 



2ttJ ■ \A7T^{AB(f)Hn^)\ 
Then, using the values given in Eqs. (12.1541) and (12.1591) . we have 

= 4.7 X 10^ GeV. (2.162) 

This value of H listed here for the Coleman- Weinberg potential can be compared with 
that found in Eq. (12.1491) to see how discrepancies can arise when choosing between 
different models consistent with observations. 

The usual method of describing inflation by first specifying a potential and then 
calculating observable quantities is thus in some ways not very constraining in its pre- 
dictions for the early universe. In the next chapter we will discuss a means of modeling 
inflation in a potential-independent way by specifying the evolution of a scale factor 
consistent with inflation instead of attempting to discern between individual models 
of potentials consistent with inflation. 
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Chapter 3 

Spectrum of Inflaton Fluctuations 



In [38j, Parker showed how to renormahze fluctuations in the inflaton field in curved 
spacetime using adiabatic regularization, for which see also [52l [53l IS] . Other papers 
[55| [56] have since found similar disagreement with the standard treatment of the 
dispersion. The technique used in [SSj has been shown to give the same results in 
homogeneous and isotropic universes as other methods of renormalization, such as 
point-splitting, and to be related to the Hadamard condition in curved space time 
[SI EH EEl EHl [60], which states that the two-point function (0 |</)(x), 4>{x') \ 0), in the 
limit x' — > X takes the form of a Hadamard Solution [15], [59] 



where a is the proper distance of interval of spacetime between x and x', A = 
—det[dadba][g{x)g{x')]^^^'^ and reduces to [—gix)]^^^"^ as x' x, and 




(3.1) 



oo 



V 




(3.2) 



1=0 



oo 



W 




(3.3) 



1=0 
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As an additional check on adiabatic regularization, we examine the spectrum of in- 
flaton perturbations in spacetimes that asymptotically approach Minkowski space at 
early and late times. This is a method introduced and used in Parker's analysis of 
particle creation by an expanding universe [HEllS], and it requires no renormalization 
beyond that already known in Minkowski space. To make use of Minkowski space 
in the analysis of the spectrum of inflaton perturbations coming from inflation, we 
investigate a scale factor, which is a measure of the size of the universe, that is com- 
posed of different scale factor segments joined together, similar to the treatments of 
[6T1 [62] . We first tried evolving forward the inflaton perturbations using a fourth 
order Runge-Kutta numerical integration routine in C++ code, but we realized that 
we would need to use greater precision for our computation. We decided instead to 
use an analytical calculation by matching known solutions to the evolution equation 
at the boundary conditions joining the different segments of the scale factor. Our 
calculations were performed using 500 digit precision in Mathematica. 

3.1 Composite Scale Factor 

We consider the metric 



The time t will run continuously from —oo to oo. The scale factor a{t) will be 
composed of three segments. Our scale factor will generally be C^, i.e., a continuous 
function with continuous first and second derivatives everywhere, including at the 
joining points between segments. Briefly, we will consider scale factors that are only 
or at the joining points. The initial and final segments are asymptotically 
Minkowskian in the distant past and future, respectively. The middle segment is an 
exponential expansion with respect to the time t. We choose specific forms for a{t) 
in these segments that have exact solutions of the evolution equations for inflaton 




(3.4) 
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quantum fluctuations of zero effective mass. 

Fig. l3.1l shows an example our composite scale factor plotted versus dimensionless 
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Figure 3.1: Composite Scale Factor. 



time. This illustrative example summarizes our notation using a moderate expansion 
of ~ 2 e-folds. The scale factor, a{t), is continuous, as are a{t) and d{t). In this case, 
the parameters for the initial asymptotically fiat segment are an = 1, a2i = 2, and 
Si = 1. The free parameters of the final asymptotically fiat segment are = 9 and 
aif = 6. Both asymptotically fiat scale factors are given by different parameter choices 
of Eq. (13. 6p with the parameter b in both cases equal to zero. The asymptotically fiat 
scale factor of the initial region joins the exponentially expanding scale factor of the 
middle region at a time ti in t-time and in r-time. The exponentially expanding 
scale factor of the middle region joins the asymptotically fiat scale factor of the final 
region at a time t2 in t-time and rj- in r'-time of the final segment, where a prime is 
used to distinguish between the r-times of the initial and final segments. 
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The equation for the middle (inflationary) segment of our composite scale factor 
is given in terms of proper time by 



a(t) = a(ti)e 



(3.5) 



where ifinfi is the constant value of H{t) = a^^da/dt during the exponential expansion 
of the middle segment. 

We define the quantity of Eq. (12.1181) . N^, = In (a2//aij), in terms of a(tinitiai) = ciu 
and a(tfinai) = '22/- When there is a long period of exponential growth, is essentially 
the number of e-foldings of inflation. Typically, Ng will be about 60. Within the final 
asymptotically flat scale factor, the ratio of 02/ to ai/ determines how gradually the 
exponential expansion transitions to the asymptotically flat late-time region. (For 
example, this ratio might be 1 e-fold, which we would consider to be relatively gradual, 
or it might be 1.0001, which we would consider to be relatively abrupt.) 

3.1.1 Asymptotically Minkowski 

The initial and final asymptotically flat regions permit us to unambiguously interpret 
our results for free fields without having to perform any renormalization in curved 
spacetime. The final asymptotically flat region will not significantly affect the result 
obtained for the spectrum of inflaton perturbations created by the inflationary seg- 
ment of the expansion. The initial asymptotically flat region should have a negligible 
effect on the spectrum resulting from a long period of inflation, although we do find 
remnants of the early initial conditions in the late-time inflaton dispersion spectrum, 
which we will discuss in Sec. 13.4.11 

We base each asymptotic segment on a scale factor of the form. 



a(t(r)) = {a/ + e^/^[(a2'-a 



^)(e^/^ + l)+6](e^/^ + l)-'}% 



(3.6) 
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where r is related to the proper time t by 

dr = a{ty^dt. (3.7) 
See Fig. 13.21 This figure shows the asymptotically fiat scale factor, a (t^r)), and the 
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Figure 3.2: Asymptotically Flat Scale Factor. 

associated dimensionless Hubble parameter, sH{t{T)) = sa^^da/dt = sa'^^da/dr, of 
Eq. (13.61) with ai = 1, 02 = 2, 6 = 0, and s = 1. Note in the graph that the maximum 
of H occurs at a value of a(t(r)) closer to ai than to 02. In both the case where 
a2 ^ cti and the case where 02 — cti, -f^max occurs at a value of the scale factor where 
a(t(r)) ~ oi. 

The form of the scale factor in Eq. (13. 6p is based on the form of the index of 
refraction used by Epstein to model the scattering of radio waves in the upper atmo- 
sphere and by Eckart to model the potential energy in one- dimensional scattering in 
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quantum mechanics [S3l El] • It was first used in the cosmological context by Parker 
[H [651 ES] to model a{t). As can be seen from Fig. 13. 2[ this scale factor approaches 
the constant ai at early times and the constant 02 at late times, and the constant s 
determines roughly the interval of r-time for a{t) to go from ai to 02- A sufficiently 
large magnitude of b would produce a bump or valley in a{t), but unless otherwise 
noted, we will take the value of b to be zero. The parameters Oi, 02, b, and s are 
different in the initial and final asymptotically fiat segments. Where confusion would 
arise we will include subscripts i in the initial set of parameters and / in the final set 
of parameters. 

3.1.2 Continuity of Joining Conditions 

With our choices of a{t) in the three segments, we are able to join them so that a{t) 
and its first and second derivatives with respect to time are everywhere continuous. 
This requires that we join the exponentially expanding segment, in which H{t) has 
the constant value i/infl, to the initial and final segments at the times when H{t) is 
an extremum. This is a maximum value, when 6 = 0, and we equate this maximum 
value of H{t) with i^mfi. A simple power law form of the scale factor, such as that 
of a radiation-dominated universe, could not be used to simultaneously maintain the 
continuity of the scale factor and its first and second derivatives when matched directly 
to the infiationary segment of exponential expansion. An application of these methods 
of matching continuously to for the radiation reaction of the electromagnetic force 
is given in the Appendix El 

3.1.2.1 Matching Continuously to Second Derivative 

With bi = and 6/ = 0, we then find the following expressions. The time Xj at which 
the first segment joins to the exponential segment is 




(3.8) 
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The constant a(ti) in Eq. (13.51) is 



a{ti) 



(3.9) 



Because the maximum value of H{t) in the first segment must equal ifinfi, we find 
that 



infl 



23/4 (_a^^4 + 



'-21 



llaij"^ - 3a2i^ + Ci) Si 



X 



(-3ai/ - 3a2i^ + Ci)^^"^ 
X (3ai/ - 3a2/ + C,) , 



(3.10) 



where 



(3.11) 



Once we choose values for aif and a2/, the remaining constants are determined to 
have the following values: 



Co 7 (lla 



2/ 

X (-3 



1/ 



3^2/ + ^f) -^infl 



^1/ 



3^2/ + Cf) 



1/4 



X 



(3ai/ - 3a2/ + Cf) 



(3.12) 



We denote the parameter r of Eq. (13.61) as r' in the final segment. At the time rj- 
when the exponential segment joins to the final segment, we find that 



r} = Sf In 



3(2i / — 3(32 / ~^ Cf 



"1/ 



(3.13) 



^2/ 



The corresponding proper time t at which the exponential segment joins to the final 
segment is 



to 



AH, 



In 



infl 



-3(3j^j4 — "icLr,^ -\- Cf 



(3.14) 
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where 



Cf = \/9a-^f + A6a, /a,/ + 9a 



1/ "2/ 



'2/ 



(3.15) 



See Fig. 13.11 for a schematic diagram of how we match our segments of the scale 
factor together. 

Fig. l3.3l shows an example of our composite scale factor and a particular dimension- 
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Figure 3.3: Matching Boundary Conditions. 



less solution to the evolution equation, where both are plotted versus dimensionless 
time. This example shows our composite scale factor over a moderate expansion of 
~ 2 e-folds. The scale factor, a{t), is continuous, as are a{t) and d{t). The parameters 
for the first asymptotically flat segment are an = 1, a2i = 2, and Si = 1. The free 
parameters of the end asymptotically fiat segment are 02/ = 9 and oi/ = 6. We 



choose bi 



0. We plot the k = 2 Fourier mode of ^/Si'^pk alongside the scale 



factor to show how this representative evolution solution changes with respect to the 
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scale factor. The real part of ^fslipki "Re," the imaginary part of y^ipk, "Im," and 
the magnitude of ^/si'^|Jk, "Abs," are all plotted. 



3.1.2.2 Avoidance of Divergent Energy Density 

We have checked our method against known mathematical theorems. One such theo- 
rem is that in an oscillator with a changing frequency, the quantity E/uj is conserved 
if the changes in frequency are made continuously in all derivatives with respect to 
time; however, if any of the derivatives of the frequency with respect to time are 
discontinuous, then this introduces changes to the conserved quantity of order A^, 
where the A^-th derivative is the first discontinuous derivative [67] . It is also shown 
by [68] that for adiabatic changes, the changes to the conserved quantity fall off with 
increasing frequency faster than any power of the frequency. We find in this conserved 
quantity a close analogy with the average number of particles created per mode for 
high-energy particles, which are those particles whose wavelengths have not yet exited 
the Hubble radius before the end of inflation. It is found in Ref. [T], that when the 
scale factor is changed adiabatically, the amount of particle production falls off with 
frequency faster than any power of the frequency. The dependence of high-frequency 
particle production upon the continuity of the scale factor is also noted in [69]. The 
scale factor must maintain continuity in the zeroth, first, and second derivatives to 
avoid an ultraviolet divergence in the energy density. This is the reason why we 
choose matching conditions that are continuous in a{t), H{t), and H{t). We could 
in principle maintain continuity in higher derivatives of our composite scale factor, 
as well, which would further reduce the amount of high-energy particle production. 
This further reduction in the high-energy particles would not appreciably improve 
upon any of our qualitative or quantitative results. The need for matching con- 
ditions when trying to calculate a finite energy density was previously realized by 
|61] . In the work of [62l [70] upon the creation of gravitons during inflation, the 
scale factor is not C^, and both authors adopt a UV-cutoff frequency. The author of 
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[52] recognizes the dependence of high-energy particle production upon the transition 
from de Sitter space to a radiation dominated universe, and he attributes the entire 
amount of high-energy particle production to the instantaneous change in the Ricci 
scalar curvature given by Eq. ( ]2.16p from 12H^ during inflation to in a radiation 
dominated universe. In [3J, Parker has shown that massless gravitons satisfying a 
conformally invariant spin-2 field would not be produced for any a{t). However, an 
Einstein graviton that instead satisfied a weak field approximation such as Eq. (14. ip . 
which in vacuum would lead to Ohab = 0, is not conformally invariant. (We use here 
the definition hat = hat — \hrjab, and we work in the Lorentz gauge where hf^^ = 0, 
which means h^p = —h^p = [IZ].) This is analogous to a massless, minimally-coupled 
Klein-Gordon field equation of the form of Eq. (I2.87p . except for the two polarizations 
(/i+ and /ix) of gravitational waves [71], [721 [73l [71] . This means that for quanta of this 
linear field, we would expect the same results for average number of quanta created 
per mode for each polarization; therefore, l/JfclLstein graviton = 2 l/^fclLiar- 

3.2 Solutions to the Evolution Equation 

Consider an inflaton field composed of a spatially homogeneous term plus a first order 
perturbation, 



We investigate, in units of /i = c = 1, a minimally-coupled scalar field that obeys 
Eq. (I2.90p . which we will refer to as the evolution equation: 

3 

dt^6(j) + 3Hdt6(j) - a-\t) ^ di^6(j) + m{(f)^^^y6(f) = 0. (3.17) 

i=l 

The mass term is related to the inflationary potential by 



0(f,t) = 0W(t) + 50(f,t). 



(3.16) 




(3.18) 



rf(0(O))2' 



For simplicity, we take m(0*^°^)^ as a constant, m^. This is an effective mass, and 
from now on will refer only to this effective mass, which may or may not be the 
same as the mass of the scalar field, which we will call mscaiar- In Eq- (I2.78p . we show 
how could incorporate a scalar coupling to the background curvature. In what 
follows, we will assume the minimally coupled case of ^ = 0, even though the term 
could include a non-zero coupling term if the curvature were also constant. (In the 
asymptotically flat segments of our composite scale factor the Ricci scalar curvature 
is not a constant.) We note that the massless, conformally-coupled case of mgcaiar = 
and ^ = 1/6 (in a 4-dimensional spacetime) would be conformally-invariant. In the 
conformally-invariant case the metric tensor and field can be deformed continuously 
at all points as 

gabix) gabix) = n{xYgabix), (3.19) 
(f){x) 0(x) = Q{xy°''''(j){x), (3.20) 

where fi(x)^ is a continuous, finite, real, scalar function; in the conformally-invariant 
case, no particle production occurs [H El [3l [151 IE]- 

The quantized field 6(f) can be written in terms of the early time creation and 
annihilation operators. At and A^, as 

^^ = T.{'^kh + 4fl)^ (3.21) 

k 

where 

h = V~h''-^^Mt{r)). (3.22) 

We are imposing periodic boundary conditions upon a cubic coordinate volume, V = 
L^. In the continuum limit L would go to infinity. The function ipk{t) satisfies 

d^Mt) + ^HdtMt) + + m^Mt) = 0, (3.23) 
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where k = 27m/ L, with n an integer. Because the creation and annihilation operators 
in Eq. fl3.2ip correspond to particles at early times, we require that ipk satisfies the 
early-time positive frequency condition 

lim Mt{r)) ~ , I -^au^.uik)r^ ^3^24) 

where ujii{k) = a/ {k/auY + m?. 

At late times, this solution will have the asymptotic form 

lim 7/>,(t(r')) , \ ke--°^/'"^/(^)-' 

+/?fce*''2^'^2-f (3.25) 

where uj2f{k) = \/ {k/a2fY + m?. 
3.2.1 Joining Conditions 

Consider a spacetime composed of three segments of the scale factor, a(t), in a homo- 
geneous background metric given by Eq. (13. 4p . For an example, see Figs. 13. II and [3731 
The first and second segments are joined at the time ti, and the second and third 
segments are joined at the time t2- 

The quantities ipk and dipk/dt are continuous across the joining regions given a 
continuity of the scale factor of at least . Using Eq. (I3.47p . it is possible to show 
the conservation of the Wronskian. Multiplying Eq. (13.470 by its conjugate leads to 



dr'^ dr 

Integrating by parts shows 



■ ,,.dipk{t)* . ...^dipkit) 



Mt) = —T^Mt)*- (3.26) 



0. (3.27) 



boundary 
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Since the boundary conditions are arbitrary, it follows with Eq. (13.71) that the Wron- 
skian, 



a(t)= 



^k{t) — i^k{t) 



(3.28) 



dt dt 

is a constant. Using Eq. (I3.24p . we see that this constant is just i; and using Eq. (I3.25p . 
we see that iakdk* — iPkPk* = or pLj 



(3.29) 



We have two linearly independent solutions to the evolution equation in both the 
second segment, with solutions hi(t) and h2{t); and the third segment, with solutions 
gi{t) and g2{t)] for a total of four separate functions. These functions are multiplied 
by constant coefficients that we must determine. During the second segment, from ti 
to t2, we have: 



Mt) = Ahi{t) + Bh2{t), 
m) = ^Kit) + BKit). 



(3.30) 



For t > ^2, we have: 



Mt) = C9iit) + Dg2{t), 
^'k{t) = Cg[{t) + Dg',{t). 



(3.31) 



If we require that ipkif) and ip'^it) be continuous at ti and ^2- This imposes 4 matching 
conditions: 



A/li(ti)+S/l2(tl) =V^fc(tl), 

Ah\{tr) + Bh'Sr) = i;'^{ti), 
Cgi{t2) + Dg2{t2) = Ahit2) + 5/^2(^2), 



(3.32) 



Cg[it2) + Dg'^{t2) = Ah[it2) + Bh'^{t2). 
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Given the values of ■0^1 and i/j'^^, and the matching conditions 



Ahi{h) + Bh2{h) 

Ah[{h) + Bh'^iti) 
Cgi{t2) + Dg2{t2) 
Cg[{t2) + Dg',{t2) 



ipkih) = ipki, 

Ahi{t2)+Bh2{t2), 
Ah[{t2)+Bh'2{t2), 



(3.33) 



we wish to calculate the constant coefficients C and D in terms of the functions hi{t), 
h2{t), gi{t), and g2{t); and the values of ipki, i^'kn ^^'^ ^2- (Here a prime denotes 
derivative with respect to t.) Rearranging the first two matching conditions leads to 



B 
A 



-0fci - Ahi 

h2 



t=ti 



t=ti 



Combining these two equations leads to 



A- 
B 



h[h2 — hih'2 
h^hi — h2h[ 



t=ti 



t=ti 



At the time t2 we have: 



(3.34) 



(3.35) 



= { 



Ahi{t2) + Bh2{t2) 

h'-Ji2 — hih'2 
h'2hi — /i2^i 



hi{t2) 



t=tl 



+ 



t=tl 



:(^2)}, 



(3.36) 
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and 



+ 



h[h2 



hih'2 



Kih) 



t=ti 



t=ti 



(3.37) 



Let us also define -^^2 = 'ipkih) and z/'L = "^K^s)- In terms of 'i/'fc2 and ilj'k2 ^^^^ 
two boundary conditions in Eq. (13.331) become 



V 9'l92-9l9'2 )t=t2 
^ f ^fc2^1 - 



Substituting for ?/'a;2 and ?/'^2 yields 



(3.38) 



C- 
D 



[Ah[ + Bh'^]g2 - [A/ii + Bh2]g'2 

9'i92 - 9192 / t=t2 

[Ah[ + Bh',]g^-[Ah + Bh2]g[ \ 

929i- 929'i )t=t2 



(3.39) 



Finally, expressing A and i? in terms of the given values of ipki and ip'j^^ specified at 
ti leads to 



C 



1 



{.9'i92 - 9i92)t=t2 

^kih2 - i^kiK 



X 



+ 



h[h2 — hih'2 
h^hi — /i2^i 



{h\g2 - /ii^2)t= 



■■t2 



t=ti 



iK92 - /i2^2) 



t=t2 



t=ti 



(3.40) 



and 
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X 



+ 



h[h2 — hih'2 
h^hi — h2h[ 



{h[gi - hig[) 



t=t2 



t=ti 



(Kgi - h2g[) 



t=t2 



t=U 



(3.41) 



which are the combined joining conditions for ipk and ip'^,. 

We find ipki and ip'j^i from the solution to the evolution equation in the initial 
asymptotically fiat segment of the scale factor. In the massless case, this solution is 
given by Eq. (13.431) . The functions hi{t) and h2{t) are to be related to the evolution 
equation solutions in the infiationary middle segment of the scale factor. Comparing 
this with Eqs. (I3^ and (l3^ shows A = E{k) and B = F{k). Similarly, the 
functions gi{t) and (72 (i) are to be related to to the evolution equation solutions in 
the final asymptotically fiat segment of the scale factor, and we will later make the 
identification C = Ni{k) and D = N2{k), where the coefficients Ni{k) and N2{k) are 
defined through their use in Eq. (I3.46p . 

3.2.2 Exact Massless Solutions 



We will first consider the case, m = 0. Rewriting the evolution equation, Eq. (I3.23p . 
in terms of r instead of t leads to 



-k^a^i/jk- 



(3.42) 



For the first segment of our composite scale factor, the solution of (I3.42p having 
positive frequency form (13.240 at early times is the hypergeometric function [H 

EHIEE] 



1 



— ^|=au ^F{-ikau^s, + ika2iS^, 
y laii'^k 

—ikaii^Si — ika2i^Si\ 1 — 2ikaii^Si\ —e~), 



(3.43) 
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where -F(a, b; c; d) is the hypergeometric function as defined in see 15.1.1]: 



F(a, h] c; z) 



Tic) 



mm 



?i=0 



r(a + n)r(6 + n) 



r(c + n) 



77,! 



(3.44) 



For the exponentially expanding segment of the scale factor in the massless case 
(V" = in Eq. fICTD above) 



Mi) = <i)~' 



2 



k 



a{t)Hi 



infl 



F{k)H 



(2) 



A; 



a(t)i7i 



infl 



(3.45) 



where H^^^ and -ff*^^-* are the Hankel functions of the first and second kind. The vari- 
ables t and r are related by Eq. (13.71) . The coefficients E{k) and F{k) are determined 
by the matching conditions of the first joining point at t = ti. We note that the finite 
period of exponential infiation lacks the full symmetries of a de Sitter universe. In 
the pure de Sitter shown in [20] , the k = mode has to be chosen in a special 

way to avoid infrared divergences. For our a(t), infrared divergences do not arise (see 
Sec. 13X2]) . 

For the final segment of our composite scale factor, the solution of the evolution 
equation ( 13.42^ is a linear combination of hypergeometric functions [H [HI [Ml [66] : 



^fc(t(r')) = Ni{k)e-""''f'^'F{-tkai/sf + tka2f''sf, 

—ikaif'^Sf — ika2f'^Sf] 1 — 2ikaif'^Sf] —e^f ) 
+N2{k)e''"''f'^'F{ikai/sf + ika2/sf, 
ikttif'^s — ika2f'^Sf] 1 + 2ikaif'^Sf] —e'f ), 



(3.46) 



where the coefficients Ni{k) and N2{k) are determined by the matching conditions 
of the second joining point at t = ^2- An example of the evolution for a particular 
mode is plotted for a specific choice of parameters using our composite scale factor 
in Figs. [331 and [331 

Fig. 13.41 shows a dimensionless solution to the massless evolution equation, where 
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Figure 3.4: A Dimensionless Solution to the Evolution Equation. 

the k = 2 Fourier mode is plotted versus dimensionless time for the same composite 
scale factor used in Fig. 13.11 The real part of ^/Si'^pk, the imaginary part of ^/Si^pk, 
and the magnitude of ^/si4'k are all plotted. 

With joining conditions for the segments of the scale factor, the derived solution to 
the evolution equation can be matched up with the known solution for the exponential 
expansion of an inflationary segment by matching 6(f)k{t) and its time derivative across 
the boundary conditions. See Figure 13.51 for the evolution of modes in the middle of 
a long inflationary period for the massless case. The time t is taken to be zero when 
k = a{t)H (when the plotted mode exits the Hubble radius) and depends on the 
mode number k. Multiplied by ^^/-f^^fiation ^^"^ plotted against this mode-dependent 
time, all of the different fluctuation modes align along the same curve in this graph. 
This shows, in the massless case, the scale-invariance of the spectrum for those modes 
that exit the Hubble radius during a period of constant H{t). 
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Figure 3.5: Scale-Invariance of Inflaton Perturbations in Continuum Limit. 
3.2.3 Approximations to Massive Solution 



In the case of a massive scalar field, the evolution equation, Eq. fl3.23p . can be written 
in terms of r as 

^'^'^ -(/tV + m^a^)^^. (3.47) 



For the middle, inflationary segment of our scale factor, our solution given by Eq. (12.991) 
is 



V 



9 2 



k 



a{t)Hi 



infl/ V 



k 



, (3.48) 



where we define rriH in terms of the effective mass by 



rriH 



infl 



(3.49) 



We know the solution to the evolution equation for the region of the scale factor 
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given by Eq. (13.51) exactly, but we do not have an analytic solution for an asymptot- 
ically flat segment of our scale factor except for the trivial case of a constant scale 
factor. We instead use one of two different approximations that we find reduce to 
the same numerical solutions in their mutual realms of applicability: the effective-k 
approach and the dominant-term approach. 

3.2.2.1 Effective-k Approach 

In the first of these approximations, the effective-k approach, we choose our initial 
and final asymptotically fiat segments of the scale factor such that an ~ a2i and 
ai/ ~ 02/. The middle segment of our scale factor, under these conditions, is thus 
where almost all of the change in the scale factor occurs, and we make use of our 
exact solution in this region. In the beginning and final asymptotically flat segments 
we make the transformation k k^s, where k^s is an effective k defined in the initial 
region as 



kes = \ k^ + m^a^^, (3.50) 



kf^s = \k^ + m^ay^. (3.51) 



and in the final region by 



In the limit that 02 = ai in a given segment, the approximation becomes exact and 
reduces to the known Minkowski fiat space solution of 

Mt{r)) = la^e-''''-^ + Ae^'^'^^l , (3.52) 

where lo is given by 



cu = W^ + m2. (3.53) 



a2 



The closer the ratio 02/01 comes to unity in an asymptotically flat segment of the scale 
factor, the more trustworthy the effective-k approach becomes. If the two parameters 
are precisely equal, however, then the scale factor becomes a constant in time and 



derivatives of the scale factor are equal to zero. In such a case where 02 = Oi, we 
cannot join to the inflationary middle segment continuously in any derivatives of the 
scale factor. When ai/ ~ 02/ in the end segment of our composite scale factor, we 
observe ultraviolet particle production due to the rapid breaking, or deceleration, of 
the scale factor's expansion. This is true regardless of effective mass, because this 
"extended" region of particle production occurs where the mass is negligible and 
{k/a{t)f » w?. 

3.2.2.2 Dominant Term Approach 

The Effective-A; Approach works very well — especially for the case where the final 
asymptotically flat scale factor is parameterized such that ai/ ~ 02/. The Effective-A; 
Approach need not be as accurate when ai/ <S 02/, and for this situation we introduce 
an alternate massive approximation, that of the Dominant Term Approach. In this 
case we introduce a new asymptotically flat scale factor that yields an exact solution 
in the limit that ^ 0. For a fixed mass, this approximation becomes exceedingly 
close to the exact solution whenever \m\ ^ k/a{t). In the Dominant Term Approach, 
when k/a ^ |m|, we use the asymptotically flat scale factor given above along with 
the massless solution; and when \m\ ^ k/a(t), we use a new asymptotically flat scale 
factor and its associated zeroth Fourier mode solution. These two solutions can be 
matched up for the case of modes in the intermediary- 52 region, where we would use 
the massless solution for the initial asymptotically flat scale factor and the massive 
solution for the final asymptotically fiat scale factor. The Dominant Term Approach 
is suspect at the interface between the small- and intermediary-g2 behaviors and at the 
interface between the intermediary- and large-g2 behaviors, where the justification for 
neglecting either the m-term or the A;/a-term is weakest. Depending upon which term 
is neglected, however, this method provides tight upper and lower limits on the average 
particle production per mode even at these interfaces. When an abrupt transition 
from the exponential infiation of the middle scale factor segment to the asymptotically 
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flat final scale factor segment is taken to make a fair comparison, the Dominant 
Term Approach is in excellent agreement with the Effective-A; Approach — even at 
the interfaces of g2 — 1 and q2 ^ exp(— iVg). When the final transition between the 
second and third scale factor segments is not taken to be abrupt, the upper- and lower- 
limits place the results of the Dominant Term Approach very close to the Effective-A; 
Approach — even at the interfaces — and they differ only in their descriptions of the 
large-g2 behavior. This is because the Effective-A; Approach requires an abrupt end to 
inflation and is not a contradiction between the two approaches, but rather is a result 
of the previously mentioned fact that an abrupt transition at the end of inflation 
produces a high-energy region of residual particle production. 

Inflaton Field of Fixed Mass and Zeroth Fourier Mode 

In units of ^ = c = 1, the perturbations to the inflaton field satisfy the evolution 
equation for mode-A; 

54)k + 3H{t)6(t)k + -T7^S(pk + m^6(Pk = 0; (3.54) 
a[t) 

where a dot represents a derivative with respect to the proper time; where a{t) is 
the scale factor; where H(t) = a{t)/a{t) is the Hubble constant, which may vary 
with time; and where m is taken to be a constant effective inflaton mass, which is 
equal to the square root of the second derivative of the inflationary potential with 
respect to the homogeneous, background part of the inflaton fleld. With a change of 
variables from the proper time, t, to a new time variable that satisfies the relationship 
dr = a{t)~^dt; and examining the zeroth Fourier mode, where A; = 0, which can in 
fact can be taken to be approximately correct whenever k/a{t) <S m, the evolution 
equation becomes 

^ = -mMrrScj>o. (3.55) 
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Using an analysis patterned after that which Epstein used to model the scattering 
of radio waves off the ionosphere [63] and that which Eckart used to model potential 
energy in one- dimensional scattering in quantum mechanics [63j, we define a scale 
factor that is asymptotically flat in both the past- and future-time infinities as 



a(r) = jai^ + e^/'[{a^' - a,<'){e^/' + 1) + 6](e^/^ + l)-^}' 



(3.56) 



The form of this scale factor is modeled after the scale factor first introduced by 
Parker [H [161 [651 [66] which has four adjustable parameters ai, 02, s, and b that allow 
one to approximate a wide range of possible scale factors a(r). The field equation, 
Eq. (]3.55p . with this scale factor, a(r), has exact solutions in terms of hypergeometric 
functions [631 [6l]. With this scale factor, Eq. fl3.55p becomes 



-m 



{a,' + e"/^[(a/ - a,')ie^/' + 1) + 6](e"/^ + ly'} 500- (3.57) 



A change of variables to m = e'^/'^ leads to 



d^ 5(1)0 
d{s\\iuY 



-m^ {a^^ + u[{a^^ - a-^^){u + l) + h]{u + 1)"^} (50o- 



(3.58) 



With the chain rule, we use 



dH 



1 f d\nu 



(i(slnu)^ 



\ du 

u d 

du 
v^d^^ 

du'^ 



-'d_ 
du 



dlio.u\ ^ d 



du 



du 



u—d(po 
du 



du 



(3.59) 



to write, with a prime denoting a derivative with respect to the variable u. 



% + ^ + ^ {ai^ + M[(a/ - a^%u +l) + h]{u + 1)-^} 500 = 0. (3.60) 



u u 
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Without having yet made any assumption as to the reahty of r/s, the variable u may 
range from — oo to +00 on the complex plane. Portions of this evolution equation 
can be seen to become infinite at m = and u = —1. For the case of m = 0, where 
the evolution equation becomes 

+ M) + '^a,'6<Po = 0, (3.61) 

we use the chain rule to change variables to v = Inu, where du = u^^d^, to get 

e-^'d^ {e-''dJ(j)o) + e-^''dJ(j)o + e-^'' s'^m^a^^6(f)o = 0, (3.62) 

which simplifies to 

dv'^6(f)o = —s'^rn?a^5(t)Q, (3.63) 

the solution of which is, 

For the case of m = — 1, where the evolution equation becomes 



^0 - ^<Po + s^"^^ {oi® - b{u + 1)-^} (500 = 0, (3.65) 



we test the analog of the solution found in Eq. (13.641) to look for a solution of the 
form 

6(Po = {u + iy, (3.66) 
and insert this into the evolution equation for the case of u = — 1 to find 



x{x - l)(u + 1)^-2 - x{u + 1)^-1 + s^m^ai^u + If - s^m%{u + If-^ = 0. (3.67) 
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Because {u + 1) — 0, the factors with the lowest exponential power of {u + 
dominate this equation, and at the point oi u — —1 the evolution equation obeys 

x{x - l){u + 1)^-2 = s^m\u + (3.68) 



or 



with solutions 



so at 14 = —1 



x(x - 1) = s^m% (3.69) 



1±V1 + 45W6 
a;± = , (3.70) 



S(Po={u + iy^. (3.71) 



A second order differential equation has at most two distinct solutions; therefore, our 
test has found all the solutions for the case of = — 1. To write the u — case in an 
equivalent form, we define 

Pi = ismai^, (3.72) 

such that for the u = case 

5(1)0 ^u^P\ (3.73) 

and define for later use 

P2 = isma2^- (3.74) 
To find the general solution of 5(l)o{u), we write 



(3.75) 
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where the function f[u] is defined by this equation. We insert this expression for 500 
back into Eq. (13.601) to get 



= {{i + uy~u-p' 



u 



(3.76) 



f ^ {ai^ + M[(a/ - a,'){u + 1) + b]{u + 1)-^} (1 + uy-u-^^fiu], 



which, with s m ai 2 = —pi,2 and s m b = a;_x+ = a;_ — x_, becomes 



— PlX_(l + m' 



/M + PliPl + 1)(1 + n)^-^i-P^-VM - + ur-u-P^-'f'[u] 



+a;_(l + uy—\-P'f'[u] -pi{l + u)"-n-^i- + (1 + m)"-^-^^"^ 
x_(l + m)^-~^m-pi/[m] - Pi(1 + m)^-m-pi-7[m] + (1 + uY~u-P'f'[u] 



u 



+ - {-Pi' + u[{-p2^ + Pi2)(m + 1) + x.x+]{u + l)-2} (1 + uT-u-P-f[ul (3.77) 



multiplying by (1 + u) ^-+i^pi+i produces 



x_(x„ — 1)(1 + m) ^m/M — pix^f\u] + a;_n/'[M] 

+Pi(pi + 1)(1 +u)u~^f[u] -pi{l + u)f[u] 
+x_uf'[u] - pi(l + u)f'[u] + (1 + u)uf"[u] 
x^uf[u\ — pi(l + u)uf[u] + (1 + u)uf'[u\ 



u 



u+l 



f {-Pi' + n[(-p2' +Pi')(m + 1) + x_a;+](^i + 1)-^} /M, (3.78) 



which can be simplified to 



= m(m+ 1)/"+ [2x_M-2pi(l +m) + (1 +m)]/ 
+ 

n, I 1 

{-Pl^ + U[{-P2^ +Pl^){u+l)+X.X+\{u+iy^] f. 



(3.79) 



X„(x_ — 1)(1 + m) — 2piX_ + X_ + [j9i(pi + 1) — ]9i](1 + «)« ^ 
M + 1 



which can be further simphfied to 
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= u{u+l)f" +[2x^u-2pi{l + u) + {l + u)]f 



+ 



— + a;_x+)(M + 1) ^ + (pi^ — pi^)u ^ 



+ (-2pix_ + pi^ - P2^ + Pi' - pi^ + x_) / 



(3.80) 



then to 



= u{u + l)f" + [2x-u-2pi{l+u) + {l + u)]f 
+ (-2pix_ + pi^ - p2^ + xJ) /, 



(3.81) 



and finally to 



= M(u + l)/" + [(2a;_-2pi + l)M+(l-2pi)]/' 

+ (X_ -pi +P2){X- -pi -P2)f- 



(3.82) 



This is a hypergeometric equation and can be solved in terms of the hypergeometric 
function / = F{x- — pi + p2,x^ — pi — p2;l — 2pi; — m), using the notation of [21] . 

Joining Scale Factors Continuously to Second Derivative 

To achieve a finite energy density we must maintain the continuity of the composite 
scale factor to at the matching points of the individual scale factor segments. 
Sec. 13.3.11 discusses further the need for joining conditions. See Figure 13.61 for 
an example of the asymptotically fiat scale factor described in the previous section 
joined to a region of inflation where the scale factor grows exponentially with respect 
to proper time. This graph shows how an asymptotically flat region could be joined 
onto the beginning or end of an exponential region. 

To join these different scale factors continuously to the second derivative, we 
note that an exponentially growing scale factor, of the form a{t) = aoexp{Ht), 
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Figure 3.6: Joining Segments of Scale Factor Continuously to C^. 
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has a time-independent Hubble constant. To find a point in the asymptotically 
flat scale factor described above where H = 0, we must find a local extremum of 
H(t). When 6 = 0, there is a unique maximum value of H{t). In a simpler scale 
factor of the form a{t) oc t", which describes a radiation- or matter-dominated uni- 
verse, no such point would exist. Using the relationship dr = a{t)~^dt, the Hub- 
ble constant is H{t) = a{t)~^{da/dt) = a{T)~'^{da/dT), and its time-derivative is 
H(t) = aiTy^dr[aiT)-^dra{T)] = aiT)-^dr^a{T) - 4a(r)-8[9,a(r)]2. This is zero 
when {d^ a{T) / dr"^) = 4:a{T)~^ {da / dr)"^ ] in other words, when 



{i - 1(1 + e-^l')e^l'] (1 + e-^'')-^{a^^ - a^^)e-'^^/' 



+ 
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s2 [a,6 + (i + e-/-)~i(a26-ai6)] 



[a^e + (1 + e-^/^)-i(a/ - a^^)] " 
4 [a,' + (1 + e-^^T' 



(3.83) 



an — a. 



:i + e--/^)-2( 



ttr, — ai e 



s [a^s + (1 + e-^/")-i(a2^ - ai^)]^ 



where the parameter b in Eq. ( 13.60p has been taken to be zero so that there might be 
a unique maximum value of the Hubble constant. To simplify this, we multiply both 
sides of the equation by 12s'^a{T)^ (a2^ — ai^)~^{l + e~'^/'^)'^e^'^/'* to get 



({4e"/^ - 2(1 + e-"/^)e2-/^} [a/(l + e""/^) + 



On — a 



3 



(3.84) 



which can be expressed as 



(3.85) 



This is a quadratic equation with two roots for e^^^. The ratio r/s is now taken to 
be real, which means e'^^'^ is non-negative; this leaves only the positive root solution 



70 



of 



(3.86) 



Once that is found, the matching conditions for r, a('r), and if are 



r 



:ln 



a(r) 



X (fli^ - 02*^ + Vaii2 + 1401602^ + 02^2) 



(3.87) 
(3.88) 



X 



12 



(3.89) 



3.3 Particle Creation 

At late times, our solution to the evolution equation will have the asymptotic form 
given by Eq. fl3.25p . The early- and late-time vacua are related through a Bugoli- 
ubov Transformation [Ij (alternately Romanized in the literature from the Cyrillic as 
Bugolubov or Bugolyubov or Bogoliubov), where the early-time creation and anni- 
hilation operators (At and are related to the late-time creation and annihilation 
operators (at and a^) through 

a^, = auA^, + f3lAl (3.90) 

where au and [3k are the Bugoliubov coefficients given by Eq. (13.251) and satisfying 
Eq. (13.291) . Because our scale factor is asymptotically Minkowskian, the meaning of 
particles at early and late times has no ambiguity. At late times, the number operator 
is 

{N^)t-,oo = {^\a\a^\Q) = m\ (3.91) 



where |0) is the state annihilated by the early-time annihilation operators Aj:. For 
the rest of this chapter, the notation \6(j)kf is defined as (0 0) = l/gP- In 

the continuum limit, this reduces to (27r)"^|'?/'fcp. Thus, \Pkf is the average number 
of particles in mode-fc created by the expansion of the scale factor from a state that 
initially has no particles [U [3] . 

3.3.1 Dependence on Mode, Expansion, and Mass 

In the absence of units, the magnitudes of k, a, H, and m have no inherent sig- 
nificance. The ratio of the Hubble radius, H^^, to wavelength, a/k, however, does 
have significance. This combination of k/aH is what we call q2 when we take the 
particular values of a = 02/ and H = ifinfl. The other relevant dimensionless ratios 
are ttlh = m/ifinfl and Ng. Transformations that simultaneously leave the values of 
k/{a{t)H{t)) and rriH intact do not change the arguments of any of the evolution so- 
lutions used in our composite scale factor. See Eq. (13.481) for the infiationary middle 
segment of our composite scale factor. For an asymptotically fiat scale factor of either 
the form described by Eq. (13.61) or the form described by Eq. (I3.56p . no matter how 
we scale a = a(r/s, ai, 02), the ratio 02/01 remains a constant; furthermore, when 
keeping the particular value of t/s fixed, H oc l/(sa^^) oc l/(sa2^). For example, 
if we multiply A; by a constant and multiply a{t) by that same constant, we don't 
change the wavelength of our mode. If we don't alter H, this rescaling won't change 
\(3g,\. When 6 = 0, we see that this transformation is 

k k * X 

ai — >■ ai * X 
02 — >■ a2 * X 

s*x'\ (3.92) 
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For a second example, rescaling k, H^^q, and m by the same factor is equivalent to 



k * y 



m 



s * y 
m*y. 



(3.93) 



This second example won't change the average number of particles created per mode, 



either. We note that in the massless case the coefficient 1 / \/'lax^k from Eq. fl3.43l) 
may change in invariant transformations, but l/Jgal^ does not change because Eqs. (13.941) 
and (13.951) contain factors that compensate for the change in A^^i. The same is true in 



the massive case under the transformation k/a{t) — > a/ {k/a{t)Y + ra^ 



In the massless case we find the following: 



For our choice of the final asymptotically fiat segment given by Eq. (13. 6p . where we 
use Eq. (I3.46P to define our functions giif) and g2{t) in terms of the relationship 
iljkif) = Nigiitij)) + N2g2{t{T)), we find the coefficients ak and (3k of Eq. (I3.25P from 



the large argument asymptotic forms [H [161 ESI [66l [2T]. With hf = 0, Ci = iksfC 



1/ ' 



and C2 = iksfa2^, we have 



cr(i - 2ci)r(-2c2) _^ Dr(i + 2ci)r(-2c2) 



r(i - ci - C2)r(-ci - C2) r(i + ci - C2)r(ci - C2) 



and 



(3.94) 



cr(i - 2ci)r(2c2) 



D r(i + 2ci)r(2c2) 



r(i - ci + C2)r(-ci + C2) r(i + c, + c2)t{c, + 02) 



(3.95) 



Recall that C and D and the functions gi{t) and g2{t) were defined in Sec. 13.2.11 A 
useful check of our method is the test of whether Eq. (I3.29p is vahdated, which we 
find to be true in all our numerical calculations. 



The variable \f3kf is the average number of particles created in the mode k, as 
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measured at late times, from the expansion of the scale factor through N^. number of e- 
folds, starting from a universe that is initially in a vacuum state that is asymptotically 
Minkowskian. We use the dimensionless variable 

q2 ^ -4-' (3-96) 

where k is the wave number, 02/ is the asymptotically flat late-time scale factor, and 
i^infl is the constant value of {a{t)/a(t)) — where the dot represents a derivative with 
respect to proper time — during the exponential expansion of the middle segment. We 
express our results using q2 instead of the wave number, fc, because we find that 1/9^2 P 
is an invariant quantity (see Fig. I3.7p . whereas depends on the arbitrary value of 
the scale factor. By 1/5^2 we refer to the average number of particles created in the 
mode given by A; = g2-f^infiCt2/- See the end of Sec. 13.4.31 for a discussion of invariant 
transformations. 

We define three regions of q2- Values of q2 ^ exp(— A^^g) are in the small-g2 region. 
Values of exp(— Ai'g) < 92 ^ 1 are in the intermediary-g2 region. Values of 1 < g2 are 
in the large-g2 region. 

Fig. 13.71 shows the average late time particle number per mode (1/5^2 versus 
I2 = k/{a2fHi^fi) for 60 e- folds of inflation. Two cases are plotted for the massless 
case based on the behavior at the matching conditions: the scale factor continuous 
in 0th, 1st, and 2nd derivatives (C^); and the scale factor continuous in 0th and 1st 
derivatives (C^). Note that in the case, 1/3^2 1^ transitions from a ^3"^ dependence 
at the end of the intermediary-g2 region all the way to a gg"^ dependence, temporarily 
parallel to the large-g2 regime, before settling down into its ultraviolet gg"^ behav- 
ior. For the wiggles near the transition from the small-g2 region to the intermediary-g2 
region at q2 = e~^'=, compare with the graph of the dispersion spectrum in Fig. 13.111 

When a{t) is or C^, i.e. when i^infi is continuous, we find numerically that the 
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particle production per mode in the small-g2 region, (g2 ^ e"^""), is 

= smh[Ne]. (3.97) 

We also find this to be the case, analytically, by taking the limit k —>■ 0. This 
analytical limit can be seen as follows. Eq. fl3.42p . in the k ^ limit tells us that 



dipk{T)/dT is constant. From Eq. (13.241) . we see that at early times V'fc(T) = l/\^2kaij^ 



and dil)].{T) / dr = —i^Jkaii^ /2 ^ in the A; ^ limit. Because dil)i.{T)/dT is both 
constant and zero, so must ipki^) be constant. Matching ipki^) and dil)k{'T)/dT with 
the late-time conditions — which do not make any assumptions about the changing 
scale factor before the late-time asymptotically flat region of spacetime is reached — 
leads to two boundary conditions: 

l/^/2k^={ak + ^3k)/^/2k^, (3.98) 



V^«iiV2 = (-^«fc + ^/?fc) Uka2fVA ■ (3.99) 



The solution to this is 



This leads to 

«fc + /3fc = e^% (3.100) 
«fc-/5fc = e-^^ (3.101) 

ak = cosh Ne, (3.102) 

(3k = smhNe. (3.103) 

In the limit of — 0, both coefficients happen to be real, and we can see that 
Eq. (13.291) is naturally satisfied. Although this result was derived in the ^ limit, 
it is valid in the massless case whenever A;/(aiiifinfl) ^ 1. This small-g2 limit holds 
for arbitrary expansions, besides those described by our parameterized composite 



scale factor, provided they initiate from a Minkowski vacuum state. We find that the 
requirement for an ahernative to the Bunch-Davies state for the k = mode in de 
Sitter space would be a consequence of taking A^^e — oo in this analytical limit. 
For at least a moderate number of e-folds, this simplifies to 

l/9..^^|e2^^ (3.104) 

The dependence in the intermediary-g2 region (e~^= ^ ?2 ^ 1) for the or 
massless case is 

\P,f^\q2-'- (3.105) 

When Ne is finite, with our composite scale factor there are no infrared diver- 
gences. For infinite inflation, where iVe ^ oo, we find the infrared divergences of a 
de Sitter universe. This problem is resolved for a true de Sitter universe in pO]. Our 
composite scale factor is different from a purely de Sitter universe in that our initial 
conditions are specified by our initial asymptotically flat region of the scale factor. 

Discontinuities in the derivatives of the scale factor at the matching points in- 
troduce additional particle production for modes in the large-g2 (or ^2 ^ 1) region. 
For the case, where the scale factor and H = a{t)/a{t) are both continuous, the 
large-g2 region goes like 

|/3^^|2 = ^452"^ (3.106) 

For the case, where the scale factor and H = d{t)/a{t) and H{t) are all continuous, 
the large-g2 region goes like 

\P,,\' = neq2'. (3.107) 

Here and uq are constant coefficients, with ~ tiq ~ (9(1/4) for a gradual end 
to inflation. For a sufficiently abrupt end to inflation, n4 and ng can be made to be 
arbitrarily large. See Sec. 13.71 

In the C° case, H{t) is not continuous, and we find quite a different behavior. 



The evolution equation, Eq. fl3.23p . may be written [3j 



77 



+ rn? 



3 / a(t)\ ^ 3 a{t) 



a{tf 4 \a{t) J 2a{t) 



Mt) = 0. (3.108) 



At the discontinuity in a{t), if we express the jump as a step function, then the form of 
d{t) picks up a delta-function contribution. Thus, there is a finite jump in dil)k{t)/dt 
across the discontinuity. The Wronskian is still conserved. In the C° case, \l3q^'f is 
proportional to gg"^ iii the small- and large-g2 regions, and it is proportional to gg"^ 
in the intermediary-g2 region. A C° scenario would suffer from both infrared and 
ultraviolet divergences, hence we will not consider it further. 

For a non-composite scale factor composed of one asymptotically flat scale factor 
defined by Eq. (13. 6p . at large values of q2 the value of 1/3^2 P frills off faster than any 
power of g2, and in terms of k we have: [U [161 ESI EG] 



sin^ (i[l - Vl + 4A;2s26]) + sinh^ [71/^5(0^2 



sinh^[7r/cs(a;^2 _|_ Q_^y^ _ sinh^[7r/cs(a;L2 _ Q^^y^ 

In the limit that A; ^ for the case of the scale factor of Eq. (13.61) . which is asymp- 
totically fiat at early and late times and has no exponential segment, we find that 
limfc^o \l3k\^ = sinh^[A'"e], where in this case Ne is In (02/01). This is the same small- 
q2 limit for the average number of particles created per mode as we found above in 
Eq. (I3.97p . The analog of the intermediary-g2 region extends over a range of lng2 
equal to 2Ne, as opposed to Ne for the particle production associated with our com- 
posite scale factor. Thus, a graph of the average number of particles created per 



2^ 



mode for a single asymptotically fiat scale factor would look similar to Fig |3.7[ except 
the region analogous to the intermediary-g2 region would be twice as long and would 
have half the slope relative to a scale factor dominated by an exponential expansion. 

In the massive case we find the following: 

In Fig. the dependence of particle production (l/^gjl^) on mass is shown for an 
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Figure 3.8: Particle Production in the Effective-/c Approximation. 



10 



expansion of 60 e- folds. The beginning and end segments are defined by = aij(l + 
10"^^), 02/ = aije^°, and aif = 0.9999a2/. The massless case can be compared with 
the plot in Fig. 13.71 which is continuous up to the second derivative of the scale factor 
to see that the two graphs are the same for 5'2 ^ 1- In this graph, however, there is 
an "extended" region of IPq^f oc ^2^^ shortly after g2 — 1 that lasts until q2 — 10^ 
before the ultraviolet behavior of {(^q^f oc ^2 ^ is seen. The term "extended" is defined 
in Sec. 13. 7[ This is due to particle creation caused by the rapid transition from the 
inflationary region to the asymptotically fiat scale factor. The two approximations, 
the effective-k approach and the dominant-term approach, give the same results with 
this particular parameterization of inflation. Both of the massive cases shown here 
produce more red-shifted particles of low momentum than the massless case. The 
case of mfj = 1/10 produces many more low momentum particles than the case of 
m'jj = 1. See also Figs. 13.91 and [3.131 

In Fig. 13.91 the dependence of particle production (l/^gal^) mass is shown for an 
expansion of 60 e-folds. This graph is different from Fig. 13.81 in that the transition 
from exponential expansion to the final asymptotic segment of the scale factor is 
more gradual, happening over about an e-fold. Thus, we use the dominant-term 
approximation. The effective-k approach, in spite of the gradual transition to an 
asymptotically fiat scale factor, overlaps with the dominant-term approach in this 
graph except very close to g2 = 1- For values of g2 ^ 1, this graph is identical to that 
of Fig. [311 

3.3.2 Limit of Negligible Mass with Respect to H 

In Fig. I3.10[ particle production as a function of q2 is plotted for 60 e-folds for both 
the massless case and the case of m = lO^^'^i/infi, labeled as m << H. This graph 
was made using the dominant-term approximation. The effective-k approach would 
overlap on this graph except very near to q2 = mn = 10^^'^. It is always the case 
that [k/alt))"^ ^ for q2 > tuh and in this region the plot of itlh = 10~^^ overlaps 
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Figure 3.10: Non-Zero Mass, Negligible with Respect to H. 



with the massless case. For q2 < rriH e^p{—Ng), relative to the mass we can take 
k = 0, and in this region of q2 in the tiny mass case of niH ^ 1, the value of l/Jgal^ 
approaches the constant (1/4)^2^^^- In the region of ttih exp{—Ne) < q2 < f^H, we 
have {k/a{t))'^ ^ in the initial asymptotically flat region and {k/a(t))'^ <^ m? in 
the final asymptotically flat region. Between q2 — m/fexp(— A^e) and q2 — exp(— A^e), 
we see \l3qj^ oc q2~^] and between q2 — exp(— A^^g) and q2 — mn, we see {(^g^f oc q2^^- 
In light of these characteristics, a comparison of Eqs. fl3.115p and fl3.116p can be made 
with consideration to where {k/a{t)Y ^ fn? and to where {k/a(t)Y ^ m^. Such an 
analysis shows that in the tiny mass limit of rriH <^ 1 , the dispersion spectrum reduces 
to the massless dispersion spectrum. The tiny mass limit bridges the transition from 
the massless case to the case of small, non-negligible rriH such as rriH = 0.01, and the 
dispersion spectra as a function of q2 for all cases changes continuously when going 
from massless to tiny mass to small mass. This is a successful check on our method. 

3.4 Dispersion Spectrum 

The dispersion spectrum is [HI [75] 



We will first consider the massless case where m = 0. See below in Sec. 13.4.31 
for the massive case. We subtract off the late-time Minkowski vacuum contribution, 
which is that part of the unrenormalized dispersion which would be present in a 




(3.110) 



Minkowski vacuum without any particles {\Pkf = for all k), to get the dispersion 
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which in the continuum hmit becomes 



' ' ' ' a2/(27r)3 
Spherical symmetry, where d^k = Ank'^dk, gives us 



^^d^k. (3.112) 



{\^<l^^\) = TT^-^ / k\Pk\''dk. (3.113) 
zir Jo 



With k = q2a2fHi^{{ and dk = (ig2a2/-f^infl, we have 



„ 2rr2 .oo 
"2/ -"infi / \fJ \2 



^'n'^a^f JO 

tt2 

infi / „ I |2 



q2\Pq2\ dq2 



2n 



/■oo 

/ q2\P,,\'dq2. (3.114) 



In the massless case, the dispersion spectrum amphtude is thus 



Z^^^IMIS^, (3.115) 



We plot Z/H^^Q in Fig. 13.111 We see that in both the case where a{t), d(t), and 
d{t) are all continuous; and the case where a{t) and d(t) are continuous; ( | | ) is 
finite without the need for any renormalization beyond subtracting off the Minkowski 
vacuum terms. When none of the derivatives of the scale factor is continuous, then 
the dispersion spectrum does not converge. 

Fig. 13.111 shows the dispersion spectrum Z / Hf^^ given by Eq. fl3.115p for our 
composite scale factor continuous in a(t), d(t), and d{t) over an expansion of 60 e- 
folds. The y-axis, Z/ H^^^, is shown multiplied by a factor of 6"^"= ; and the x-axis, q2, is 
shown multiplied by a factor of e^^. When using this scaling, the region plotted in this 
graph would look identical for an expansion of 10 e-folds, and it would look identical 
for an expansion of 80 e-folds. In the case of a2i = an + w, where w = lO^^^aij, 
we see marked peaks in the dispersion spectrum. When we change the parameters 
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in the initial asymptotically flat region to = lOaij, these peaks are damped as 
shown. The ending conditions of the final asymptotically fiat segment do not affect 
these peaks. 

A calculation of the dispersion spectrum in the massive case leads to an equation 
analogous to Eq. fl3.115p : 



3.4.1 Sensitivity to Initial Conditions 

We take for our initial conditions a quantum state to be asymptotic at early times 
to that of a Minkowski vacuum spacetime for all modes. This is a consequence of 
our asymptotically flat scale factor and our assumption that no particles are initially 
present. It is more common in the literature to take instead the Bunch-Davies state 
for quantum fluctuations, that is to assume a de Sitter spacetime. As pointed out by 
[62j . this leads to an infrared divergence of the two-point function, where the two-point 
function is another name for our dispersion spectrum, and the cause of this divergence 
is correctly diagnosed as being due to the choice of initial conditions in [70] . Both of 
the authors of [621 [70] handle these infrared divergences with a cutoff frequency that 
omits modes that are currently outside the Hubble radius of our observable universe. 

The use of de Sitter initial conditions is equivalent to supposing an inflationary 
period that extends over an infinite number of e-folds, or iVg — oo. If we assume 
a finite Ng, and if we assume that in the future our universe will be approximately 
matter-dominated for all times, which means neglecting any dark energy or cosmo- 
logical constant, then eventually every mode that exited the Hubble radius during 
inflation would eventually re-enter the Hubble radius of our universe after inflation if 
it has not already done so. 

Both Figs. 13. Ill and [3TT2I show additional peaks after the primary peak, where the 
primary peak roughly indicates the interface between small-g2 and intermediary-g2 
behavior. These minor peaks are caused by phase differences between modes with 



Z = 




2 tt2 



(3.116) 



similar wavelengths as they exit the Hubble radius near the beginning of inflation. 
The modes that exit the Hubble radius with a large amplitude — either a positive real 
amplitude, a negative real amplitude, a positive imaginary amplitude, or a negative 
imaginary amplitude — quickly have this large amplitude translated into a near con- 
stant value outside of the Hubble radius. Those modes that exit the Hubble radius 
with relatively small amplitudes are frozen into evolutions of relatively small mag- 
nitudes outside of the Hubble radius; these relatively low-amplitude modes have a 
relatively high change in amplitude with respect to time, but this initial excess in 
the derivative of the amplitude with respect to time is rapidly redshifted away dur- 
ing inflation. With an abrupt transition from an asymptotically Minkowski vacuum 
to an exponential inflation of the scale factor, by which we mean that a(ti) ~ a^, 
where a{ti) is the scale factor at the transition from the initial asymptotically flat 
segment to the exponentially growing segment of inflation, and where an is the scale 
factor at early times, we see that the minor peaks are more pronounced. With a more 
gradual transition from the initial asymptotically flat segment of the scale factor to 
inflation (when a{ti) ^ 1.2aij), these minor peaks are damped out. If these modes 
were observable in our universe, that is if they have already re-entered our Hubble 
radius, their measurement might tell us something about initial conditions before the 
beginning of inflation: whether there had been a phase transition from the very early 
universe to inflation, how rapidly the very early universe had been expanding (or 
contracting) relative to the expansion of inflation, and what the dominant contribu- 
tion to the evolution of our universe might have been before the start of inflation. 
Because measuring the contribution of these minor wiggles to the scale dependence of 
large-scale structure would be experimentally challenging (if not impossible), this is 
in some sense speculation, but that does not change the fact that the two dispersion 
spectra shown in Fig. 13. Ill are different, and this difference — if observed — would tell 
us about our pre-inflationary universe. 



3.4.2 Sensitivity to Sub-Planck Length Physics 

Consider a quantum fluctuation of tlie particular mode that, at the beginning of 
inflation, has a wavelength equal to the Planck length. By the time this wavelength 
has been stretched to the point that the mode is exiting the Hubble radius, it will 
have a wavelength the size of the Hubble radius. For this to happen, the scale factor 
must increase by a factor of -ffj^fl/^pianck- 

With n = c = 1, the Planck length is £pianck = VG = 8 x lO^^o (GeV)"\ Using 
the value of Hi^a = 7 x 10^^ GeV given in Eq. dmoD . we find f^^fl/^ianck ^ 2 x 10^ 
which corresponds to a mode exiting the Hubble radius In(lO^) ~ 12 e-folds after 
the start of inflation. All higher frequency modes, that is for q2 > e^^"^^"^, will have 
originated from trans-Planckian modes during inflation. With Ne = 60 e-folds of 
inflation, if we use the estimate of the number of e-folds before the end of inflation 
in which the observable modes of the CMB are exiting the Hubble radius given by 
Eq. fl2.14:-il) (50 e-folds) or by Eq. (l37[44D (53 e-folds), then it might be possible to 
observe the difference in amplitudes between those modes that were initially super- 
Planckian quantum fluctuations and those that were initially sub-Planckian quantum 
fluctuations before the start of inflation. With either a smaller value of the Hubble 
constant during inflation or with a larger number of total e-folds of inflation, the 
re-entry of the first trans-Planckian modes back into our Hubble radius after inflation 
could be postponed to epochs of our universe much later than recombination. 

3.4.3 Model in Terms of Expansion and Mass 

Fig. l3.12l shows a comparison of the dispersion spectrum (Z/H^^q) given by Eq. (13.1161) 
and normalized to 1 for our composite scale factor continuous in a{t), a{t), and a(t) 
over an expansion of 60 e-folds for various masses. The values of Z/ Hf^^ were divided 
by the maximum value of the primary peak for each located at q2 — exp(— Ai'e). To 
normalize these peaks, Z/ Hf^^ was divided by the following factors: 1.3 x 10^^ for the 
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massless case, 2.3 x 10^^ for = 0.1, and 2.2 x 10'* for = 1. 

The dispersion spectrum is plotted for three different cases of tuh in Fig. 13.121 
The effective-k approach is useful for this approximation. This approach demands 
that in the initial asymptotically flat segment of the scale factor, a{t) must always be 
approximately equal to an. Because a{ti) ~ ai if either a2 — ai or 02 ^ ai, however, 
this approach can be used with a wide range of initial conditions. Specifically, when 
02 > ai in Eq. ([M]), we have a{ti) ~ (7/3) ^^/'^^ Although we are not at the moment 
considering the case of 02 ^ cli in Eq- (13.561) . for comparison we note that it would 
lead to a(ti) ~ 3*-^/^^ In both cases a(ti) — 1.2ai. We have found that, even in 
the massive case, the observed humps are dependent only upon the initial conditions. 
In the region shown in this figure, the graph would not be significantly altered by 
using joining conditions instead of our matching conditions. The effective-k 
approximation plotted on this graph would overlap with the exact solution, if an 
exact solution were available. The shapes of the curves are fixed above a moderate 
number of e-folds. We define the variable J such that the maximum value of Z / H"^^^ 
for the major peak, which is the peak located nearest to q2 = e~^'=, is J e^^'^^^" 
in the massless case and is J e^^'"^^^" in the massive case. Then, the normalization 
factor scales like 6*^-^"^)^= in the massless case, as can be seen from Eq. (13.1151) : and 
the normalization factor scales like e^^'"^^^" in the massive case, as can be seen from 
Eq. (13.1161) . where we define the exponent P in the following way: 

\^.?^\q2~'' (3.117) 

in the region of intermediary-g2 (e~^'= ^ q'2 ^ 1), and 

1/3,/ -^e^^^ (3.118) 

in the small-g2 region (^2 ^ e~^=). The exponent P is well described by a q2- 
independent value in the case of m = and in the case of 0.01 ^ '^jl ^ 9/4. 



In the massless case, P = 2, so the height of the major peak in the graph of the 
massless case in Fig. 13. 12] grows with an increasing number of e-folds as e^=, while the 
widths of the peaks narrow with an increasing number of e-folds as e"^". The area 
under an individual peak in the massless graph therefore does not change appreciably 
when changing the number of e-folds of expansion, provided there are at least a few 
e-folds of inflation. For the massive cases, we see that P = 2.93358 when m^l = 0.1, 
and that P = 2.23607 when niH = 1. 

Fig. 13. 131 shows the dependence of the variable P, as defined in Eq. 03.1170 . upon 
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Figure 3.13: Inflaton Spectrum Characterized in Terms of Infiaton Mass. 

ttlh = m/ifinfl. The calculated data points shown lie on the curve 
Outside of the region plotted, however, P does not have a constant, g2-independent 
value. For mn > 1-5, the argument, a/(9/4) — m^, of the Hankel functions becomes 
imaginary, and \Pq2f oscillates with changing g2- For an example of a non-zero mass 
much smaller than -ffinfl, see Fig. 13.101 
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Table 3.1: Approximation 



0.01 <mH< 1.49 



1 J ,ip-m. ^ x^.^-^^^T^ 



niH = 



I J + J^-N, dq2 J ql 



0.01 <mH< 1.49 




mH = 



We wish now to approximate the dependence of the configuration space dispersion 
( I (50^ I ) / H^^Q with regard to the number of e-folds. In our approximation we neglect 
the minor peaks; we assume that the major peak is located exactly at q2 = e'^", that 
Z/H^ increases linearly with q2 up to the major peak, and that Z/H^ decreases as 
{q2e^^Y~^ in the massive case — or as (^26^")^"^ in the massless case — until the 
onset of large-g2 behavior at q2 = 1, which effectively serves as a cut-off point. The 
maximum of the major peak is given by height = J ei^-"^)^^ ^ in the massive case; and 
height = J e^^'^^^", in the massless case. We find J ^ 0.01 for all three cases. In this 
simple approximation, the configuration space dispersion is given by TABLES 13.11 
and l3.2[ The small mass limit of P ^ 3 and the massless case of P = 2, both reduce 
to the same limit of 



For further discussion of the small mass limit reducing to the massless dispersion 
spectrum, see Fig. 13.101 

3.5 Scalar Spectral Index and Scale Invariance 

We define a given mode of S(j)k to be crossing the Hubble radius when k/ {a(t)H(t)) = 
1. We define a mode of k to be inside the Hubble radius when k > a{t)H(t), and 
we define a mode of k to be outside the Hubble radius when k < a{t)H{t). Modes 
in the intermediary-g2 range exit during infiation to eventually re-enter the Hubble 




(3.119) 
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radius at some time after inflation has ended, provided any cosmological constant or 
dark energy can be taken to be negligible. Using our composite scale factor, we note 
that after a few e-folds of inflation, the quantum perturbations that are exiting the 
Hubble radius are found numerically to satisfy 

where \ipkf is the time-dependent part of as given by Eqs. (13.211) and fl3.22p . 

The variable D{mH) ~ (1 + |mj|)^ is a constant of order 1 that we have evaluated 
numerically to be 

D{mH = 0) = 1.00, 
D{mH = Vol) = 1.04, 

D{mH = 1) = 1.45. (3.121) 

Thus, our spectrum of l^f/'fcl^, if evaluated at the time of exiting the Hubble radius, 
is scale- invariant, regardless of effective mass. By Eqs. (I2.124p and (I2.125P we have 

Vn^k^\S(pkf. (3.122) 
The scalar spectral index given by Eq. (12. 1270 is 

= 1 + ^Ynk ' (3.123) 

We see that when taken at the time of crossing the Hubble radius, the spectrum, 
which is proportional to k'^ \6(f)kf, has no k-dependence because we have shown the 
spectrum is proportional to k^ Hf^^/ D{mH)- Evaluating the scalar spectral index at 
the time of exiting the Hubble radius thus leads to = 1, which can be used as the 
definition of a scale-invariant spectrum. 



The modes that exit the Hubble radius at the very beginning of inflation, however, 
along with those that exit the Hubble radius before reaching the middle segment of 
our composite scale factor where a{t) begins to grow exponentially with respect to t, 
are not described by Eq. ( ]3.120p . These modes in the small-g2 region are not scale- 
invariant; therefore, well after the end of inflation, long- wavelength modes that are 
not scale-invariant would eventually re-enter the Hubble radius of a matter-dominated 
universe. If the total number of e-folds of inflation is sufficiently small, it would be 
possible to observe a transition from the scale-invariance to a scale-dependence of 
large-scale structure. See Figs. 13.71 13.81 and 13.91 Because the small-g2 modes of large 
enough wavelength exit the Hubble radius before evolving away from the early-time 
conditions specified by Eq. fl3.24l) . we would expect a massless infiaton to generate 
a spectral index of = 3 in the small-g2 region, and we would expect a massive 
infiaton to generate a spectral index of n<j = 4 in the small-g2 region. If scale- 
invariance continued indefinitely for large wavelength modes, the dispersion would be 
infrared divergent, so this eventual end to scale-invariance is not an artifact of our 
initial conditions. The modes responsible for the galaxy-size structure of today left 
the Hubble radius approximately 45 e-folds before the end of infiation P p. 285], so 
if Ne were not too much larger than this, we would expect it to be possible to measure 
the end of scale-invariance in our observable universe. 

In Fig. 13.141 we plot six scenarios depicting the behavior of iV'fcl^ after exiting the 
Hubble radius. The first two cases, A and B, are for m^/ = 0. In both of these cases, 
an expansion of 20 total e-folds is plotted. In case A, there is a gradual end to infiation 
spanning one e-fold; and in case B, there is an abrupt end to infiation. Because the 
mode has exited the Hubble radius, neither of these end conditions changes 
and the two lines overlap. Here, and in general for the massless case, reaches a 

constant value a few e-folds after exiting the Hubble radius, and this constant value 
is close to the value at the time of exit. In the massless case, we find a scale-invariant 
spectrum even when the spectrum is defined in terms of the value of at the 
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end of inflation. This can be found by noting that at the time a particular mode is 
crossing the Hubble radius, the value of \ipkf given by Eq. fl3.45p is approaching a 
constant value as the argument k/[a{t)Hi^f[] becomes much less than 1. 

For cases labeled C, D, E, and F; we use tjih = VO.l. In the massive cases, I'lpkf 
never reaches a constant value, although it changes much more slowly after exiting the 
Hubble radius. Cases C and D are the massive analogs of cases A and B, respectively. 
In case E, we end inflation gradually over the length of one e-fold, starting just as our 
specific mode crosses the Hubble radius. In case F, we end inflation abruptly just as 
our specific mode crosses the Hubble radius. 

From j2l|, {Eq. 9.1.9}, we see that in the small argument limit of the Hankel 
functions 

^ (.)|^^ (iMy (3.124) 

when the real part of the parameter v is positive and non-zero. In Eq. (13.481) . we find 
numerically that E(k) ~ —{i/2)^J'n/Hi^p^ and F{k) ~ for modes of intermediary-g2 , 
which are the modes that exit during the exponential expansion of our composite scale 
factor. Long after these modes have exited the Hubble radius, when k/[a{t)H\^Q\ <^ 



2 



oc a where z 



1, we expect Eq. (I3.48P to approach \ilJk\^ ~ a ^ 
/c/[a(t)-?/infl] and v = a/(9/4) — m^. Using this small argument approximation with 
Eqs. (KVm and (1371231) leads to 

dlnVn 



dink 



3-i/9-4mj|, (3.125) 



under the assumption that we evaluate |(5</)a;|^ at late times, in which case 



= 4- ■y/9-4m^. (3.126) 



In determining the spectrum, is often evaluated at the time when a mode 

exits the Hubble radius [SU [10]. In this case, we would get = 1, exactly. The 
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WMAP results [42j find rig ^ 0.96, which would suggest a value of tuh — i/^ to go 
along with our assumption of a constant value of -ffinfi, provided the end of inflation 
is the appropriate time to evaluate Vn- This value of mn — i/^ is what we found 
in Eq. f l2.158p for the Coleman- Weinberg potential. An imaginary mgcaiar would lead 
to tachyonic behavior [51j, but here m is an effective mass, so this need not be a 
problem. 

A different method of calculating Eq. (13.1261) . involves combining Eq. (I3.117P and 
q2 = k/a2fH[^f[ with Eq. (I3.132p . To renormalize Eq. (I3.132P would involve dropping 
the Minkowski vacuum term such that (1 + 2 \Pkf) 2 \Pkf; although, in the case 
of the intermediary-g2 modes, the term to be subtracted off is already negligible 
compared with the particle number per mode. To simplify the massive case, we treat 
|m| ^ k/a{t) in the intermediary-g2 region of modes that exit the Hubble radius 
during inflation. Then Eqs. (I3.122P and (13.1230 . together with the relationship given 
in Fig. 13.131 of P = a/9 — ^m^, or P = 2 in the massless case, give us the same result 
as in Eq. fl3326|l . 

If ifinfl were not constant, but were slowly decreasing during inflation, then we 
would find a red-tilted spectrum. We could incorporate this effect into our exact 
calculation by taking the adiabatic approach and using the value of ifinfl(ti) for our 
first matching conditions and the value of i/infi(i2) for our second joining. Combining 



Eqs. f l2:T25D . ([MI]), dMl, fl3:T20D . and (133231) . we find 




(3.127) 



which, with Eq. (12.1140 . becomes 



d 



(61n[i7i,fl(t)]-21n[y']), 



(3.128) 



ns = l + 



din A; 



where a dot denotes a derivative with respect to time, and a prime denotes a derivative 



with respect to 0. Then, using d/d\nk = H ^d/dt, we have 
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(3.129) 



which, through the chain rule and with Eq. fl2.114p . V = (pV' 



= -1/'1/"/(3ifinfl(t)), so that, finally, with Eq. (|2TT3|) . we have 




(3.130) 



which we write in terms of Eqs. ( I2.116P and fl2.117p to get 



= 1 - 6e + 2r], 



(3.131) 



which is equivalent to Eq. fl2.132p first shown by |^T]- 

3.6 Density Perturbations 

Fig. 13.141 shows that the maximum difference between the late-time values of \S(f)k\ 
in all six of the cases plotted is about 40%. We conclude that when uih <^ 1, the 
value of at late times is a reasonably good indicator of the value of 

at Hubble radius exit. For the rest of this section we will adopt the assumption 
that the late time value of |50a;|^ is indicative of the value of at the time of 

exiting the Hubble radius. This assumption allows us to extrapolate our method of 
late-time renormalization in Minkowski space to a time of curved spacetime in lieu 
of applying a more rigorous analysis that would require a more complex method of 
curved spacetime renormalization such as in |38j . 

The final conditions do not affect the value of much once a given mode has 

crossed the Hubble radius. Thus, we could end inflation just after a mode has exited 
the Hubble radius to flnd that the value of \6(j)kf will be very close to its late-time 
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value. At late times, Eqs. fl3.22l) and (13.251) show that the time averaged expectation 
value — 



2L'^a2f^uj2f 



{l + 2\(3kf). (3.132) 



This value of {\6(j)kf) obtained from Eqs. (I3.2ip . (I3.22p . and (I3.120p . however, is 
un-renormalized. To use the renormalized values, we take (1 + 2 l/^gj^) — * 2 1/5^2 

Although [9l p. 285] identifies the scale factor, agai = ^2/ e""^^, as the one in which 
the k-modes responsible (by seeding the density perturbations) for the formation of 
galaxies are exiting the Hubble radius; we note that when rriH <^ 1 there is a relative 
constancy of \S(j)kf after a mode crosses the Hubble radius, and thus our subsequent 
method is widely applicable to the range of intermediary-g2 modes. In our assumption 
described above, a mode defined by g2 = 1 at late times is an excellent indication of 
the state of any mode just after crossing the Hubble radius when uih ^ 1. We can 
assume for the moment that inflation ends abruptly just as the mode k = a2/i^mfi 
exits the Hubble radius. This abrupt ending does not change 1/3^2 1 ^ the q2 = 1 
mode, because we have found that the ending conditions do not affect modes of 
q2 < 1. In this case, 150^2 1^ isn't changing from its value at Hubble radius crossing 
(or is roughly equal to the late-time value it would have reached a few e-folds after 
crossing the Hubble radius), the late-time value of IPq^f for the g2 = 1 mode isn't 
changing (because there is no more inflation and the mode g2 = 1 is insensitive to 
other factors), and the scale factor isn't changing; therefore the renormalized value of 
|50g2l^ is not changing. This argument wouldn't hold for modes of large-g2, because 
they are sensitive to the time-derivatives of the scale factor, but we find that the 
late-time dispersion spectra for the mode g2 = 1 is a good approximation to the 
renormalized value of 1 50^2 1^ the time any mode exits the Hubble radius. For 
an analysis of the instantaneous renormalized value of 6(f)k that does not rely on a 
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late-time argument, see [55] . 

We next consider the curvature perturbation given by Eq. fl2.123p and defined at 
the time of Hubble radius crossing as 

TZk = -^5^k. (3.133) 

The variable is the rate of change of the homogeneous background scalar field. 
The quantum perturbations we have considered so far, 50^, are assumed to be much 
smaller in magnitude than the zeroth-order field. 

3.6.1 Hybrid Combination with Slow Roll Approximation 

So far our method in this chapter has not been linked to any particular potential or 
model of infiation. In what comes next, we choose a simple potential, which is found 
to be in good agreement with the 3- Year WMAP data and we use a hybrid 
combination of our method and the slow roll approximation. The remainder of this 
section is intended to be of a more speculative nature than the rest of this dissertation. 
For our example, we use the Linde quadratic chaotic-inflation potential [ISl HH] 



V = -mV^- (3.134) 
2 



From Eqs. (12.1131) and (12.1141) . the two slow roll conditions are 



^ ^V, (3.135) 



and 



• dV/d(f) 
d) ~ 



(3.136) 
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Table 3.3: Comparison of 5h for V = im^0^ 





H = 10^2 GeV 


H = 10^4 GeV 


H = 10^6 


mn = 0.0001 


2.263 X 10-4 


2.263 X 10-2 


2.263 X 10° 


tuh = 0.01 


2.263 X 10-6 


2.263 X 10-4 


2.263 X 10-2 


rriH = 0.1 


2.258 X 10-^ 


2.258 X 10-5 


2.258 X 10-3 


tuh = 0.25 


8.917 X 10-^ 


8.917 X 10-6 


8.917 X 10-4 


rriH = 1 


1.903 X 10-s 


1.903 X 10-6 


1.903 X 10-4 



We combine these two sfow roll equations with the potential specified in Eq. ( 13.1341) 
to find 

fo 1 

(3.137) 



'2 1 
~ —m\l 



3 VS^' 



We rewrite this as 



n niH 



2 1/V8^ 



3 -f^infi 

In our notation, with 5(j)k taken from Eq. (13.1161) . 



7^. = - 



infi 



2 I/VSttG 



q2 1 



27r2, 1 + 



(3.138) 



(3.139) 



therefore, with l/v^S^ ^ 2.436 x 10^^ GeV, 



7^i 




m^^l + m^ V 2.436 x lO^s GeV 



(3.140) 



The magnitude of the curvature perturbation has been shown to be a conserved 
quantity outside of the Hubble radius [35l|36], and the curvature perturbation can be 
related to the amplitude of density perturbations at the time of re-entry, when once 
again k/[a{t)H{t)] = 1. In a matter-dominated universe this relationship is [3] 



= Ok — -K-k- 

P 5 



(3.141) 



See TABLE 13.31 for sample values of 5h, the density contrast defined in |34j, at 
the time of re-entry into the Hubble radius and for the potential given by Eq. (I3.134p . 
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3.6.2 Relative Constancy of Modes Outside Hubble Radius 

In Fig. 13.141 the value of \ipkf, in units of Hl^^/k^ and for the mode k = {a{ti)e^^)Hi^f[, 




J I I l_ 

5 10 15 20 




Figure 3.14: Modes Exiting the Hubble Radius. 

is plotted versus dimensionless time H\^f[t. The graph shows the relative constancy of 
for modes that have exited the Hubble radius during inflation. For additional 
discussion of this graph and the difference between cases A-F, see Sec. 13. 5[ 

The inflaton perturbations only approach a true constant well outside the Hubble 
radius for the massless case. During an exponential expansion in the massive case, 
mil ^ ^/[a'(t)-ffinfl], and we may rewrite Eq. (13.231) as 



(3.142) 



The two linearly independent solutions to this are 
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ip± oc exp 



-^(3±P)i/inflt 



(3.143) 



where P is defined as in Fig 13. 131 In the small mass limit of P — > 3, well outside the 
Hubble radius one of these linearly independent solutions approaches a constant value 
with respect to t, while the other solution decays exponentially. With itih of order 1, 
both linearly independent solutions decay exponentially outside the Hubble radius. 
Although the magnitude of these massive perturbations are constant when the scale 
factor is constant in time, the rates of their decay well outside the Hubble radius 
depends on how the scale factor is changing. For our composite scale factor with an 
abrupt end to infiation, where almost all of the expansion occurs in the exponentially 
growing segment of our scale factor, there is more small-g2 particle production than 
in our composite scale factor with a relatively gradual end to infiation, where more 
of the total expansion of the scale factor takes place in the final asymptotically fiat 
segment of our scale factor. We turn now to tracing a particular mode as it exits the 
Hubble radius until it re-enters our observable universe. 

Consider, as an example, the fc-modes responsible for large-scale structure for- 
mation. As an approximation, take the following three epochs to be simultaneous: 
recombination (the time light was emitted from the surface of last scattering), the 
transition from a radiation-dominated universe to a matter-dominated universe, and 
the re-entry of the modes that would provide the density perturbations to seed galax- 
ies. Furthermore, also as an approximation, assume a transition to a radiation- 
dominated universe, where a{t) = Ct^^"^, immediately after the end of infiation such 
that H{t) is continuous. Call the time of the end of infiation tj, and call the time 
of re-entry and recombination tr- Turner and Kolb give the temperature of infiation 
and the temperature at recombination as 10^^ GeV and 1 eV, respectively [9]. In a 
radiation-dominated universe, the energy density — neglecting particle production — 



is related to the scale factor as prad oc a(t)~^, and the temperature is related to 
the energy density as T oc pl[^, so the temperature is related to the scale factor as 
T oc a{t)~^ after radiation and matter have decoupled and are no longer in ther- 
mal equilibrium. Thus we know that a{tr) = 10^'^a(t/). The radiation-dominated 
scale factor then gives us {tr/tff^ = 10^3, or tr = lO^Hf. The hubble constant 
in the radiation-dominated universe is H{t) = a{t)-^da{t)/dt = (2t)-^ Because 
H(tf) = i/infi, we have H(tr) = lO^^^ifinfl. When a mode exits the Hubble radius 
during inflation, we have k/{a(t)H(t)) = k / {a{t f)e~^^ H^-^q) , where the variable is 
the number of e-folds before the end of inflation at which a mode exits the Hubble 
radius. When our example mode re-enters the Hubble radius after inflation, we have 
k/ {a(t)H(t)) = k/ {a(tr)H(tr)) = 1. By equating the relations for exit and re-entry, 
we have k/ {a{tf)e-^-^ H;^^) = k/{10'^^a{tf)10-^^Hina), or e"^-^ = 10-^3. This means 
in our approximation 

Ke - 53 e-folds. (3.144) 

Turner and Kolb find, with a more detailed calculation, a value of 45 for this number 
P, p. 285]. The simplification of treating recombination, matter-radiation equality, 
and galaxy seeding as concurrent is a relatively useful approximation. The radiation- 
dominated universe transitions to a matter-dominated universe at a temperature 
roughly one order of magnitude higher than the temperature of recombination, which 
means the mode that will later re-enter the Hubble radius at the time of radiation- 
matter equality exits the Hubble radius during inflation roughly 2 e-folds later than 
the mode that will later re-enter the Hubble radius at recombination. The exact rela- 
tionship between these two events with the Hubble crossing for the modes responsible 
for seeding galaxy formation depends on the nature of dark matter: the current size 
of galaxies does not lead to a simple estimate of their size in the past, because their 
size does not scale with the size of the universe once they have become gravitation- 
ally bound. Baryonic matter will clump to structure initiated by cold dark matter, 
but not until after recombination, when radiation pressure overcomes gravitational 



attraction; dark matter will start clumping earlier than this, at the epoch when it 
decouples from the dominant radiation background [Qj [TOl [31] . The approximation 
of an immediate transition from inflation to a radiation-dominated universe is less 
certain, as the validity of this approach could vary based on the specific inflationary 
potential being considered. 



3.7 Reheating 

Our analysis of this particle creation reveals a mechanism for Reheating, which is 
a return to the temperatures and densities that are responsible for the successes of 
the Big Bang model. We find that the energy density present after inflation depends 
on how abrupt the transition is from the inflationary middle segment of exponential 
growth to the final asymptotically flat region of the scale factor. 

3.7.1 Energy Density from Abrupt End to Inflation 

Our scale factor can be made to be continuous to the scale factor and two of its 
derivatives, but no more, so we see additional particle production caused by discon- 
tinuities of higher derivatives. When we maintain continuity of a{t), a{t), and d{t), 
the particle number is proportional to q2~^ for large-g2- The energy of a particle of 
mode-fc at late times is uj2f = \J {kla2ff' + m^- The energy per mode in the large-g2 
regime is then proportional to {k / a2f)q2~^ = q'2~^-^mfl- 

With a gradual transition between segments of a{t), the large-g2 behavior in which 
\f3q2f falls off as q2~^ starts around q2 — 1. With an arbitrarily abrupt transition from 
the end of inflation to our final asymptotically flat scale factor, however, this transition 
can be prolonged to an arbitrarily high value of q2, which we denote by g2cut-ofr- We 
find empirically that g2cut-ofr — «2//(«2/— ^i/). We define the region between 1 < ^2 ^ 
5'2cut-ofT as the "extended" region. In the "extended" region the fall off of \Pq2f oc q2^^ 
is extended from ^2 — 1 to larger values of q2, such as the value of g2 — 10^ shown in 
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Fig. 13.81 in which 02/ — cii/ <^ 02/ as a{t) makes a rapid transition to flatness. This 
extension is caused by the production of particles of higher momenta by the rapid 
change in H{t) after inflation. When the transition of a{t) is gradual, one finds, as 
in Fig. 13. 9[ that beyond q2 — 1, the quantity \Pq2f falls off more rapidly, eventually 
going as q2~^ if the function a{t) is C^. With sufficient extension, the particle number 
per mode in the "extended" region is proportional to ^2"^) regardless of the value 
of P in the intermediary-g2 region, so for both the massless and massive cases the 
contribution to the total energy density is dominated by these "extended" modes, and 
we neglect both the red-shifted modes and the ultraviolet modes. When a(t) is and 
there exists a significant "extended" region, the contribution to the energy density 
from values of q2 > q'2cut-ofr is negligible. The energy density associated with the 
"extended" region, which dominates the total energy density when 02/ — ai/ <^ 02/, 
is 




(3.145) 



When g2cut-off ^ 1, we find for the energy density 




(3.146) 



Because g2cut-ofr — ('^2// [(^2/ ~ c^i/])) we can see that an abrupt end to inflation can 
lead to energy densities large enough to produce reheating. For particle production 
as the cause of reheating, see also [HI] . 
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3.7.2 Associated Temperature 

In units oi h = c = ks = I, the temperature is T = {{E) / a)^/^ ^ where a is the 
Stefan-Boltzmann constant. Then the energy density attributable to an abrupt end 
to inflation given by Eq. (13.1461) leads to an effective temperature of 



T^J^^lNL^. (3.147) 
^' a2f-aif 20F ^ ' 

This approximation holds for any relatively abrupt transition and does not depend 
on any discontinuities of the scale factor. 

In an expansion governed by the asymptotically flat scale factor of Eq. (13.60 with 
no exponential middle segment, the large-fc behavior — in both the massless case and 
the effective-/c approach — follows a thermal spectrum given by [H [161 ESI |66] 



4:7Tsa?a2 



When ai ~ 02, we use 



4 

1 _ % 



i^max ^ ^7-^, (3.149) 

16a2 s 

to show that in our notation this is equivalent to 

vr(l-^) 

for a single asymptotically fiat scale factor with ai ^ a2. In the large-g2 regime of 
our composite scale factor with ai/ ~ 02/, we would expect to find the temperature 
approaching this same value, regardless of mass, of P, and of the number of e-folds; 
but only if we were able to maintain continuity with the previous segments of the 
scale factor across an infinite number of derivatives. 

With a gradual transition between segments of a{t), the large-g2 behavior in which 
\Pq2f falls off as q2~^ starts around q2 — 1. For such a gradual transition, we find a 



late-time temperature — which is red-shifted after the end of inflation by the expan- 
sion of the final asymptotically flat segment of the scale factor — that is comparable 
to the Gibbons-Hawking temperature of H/{2tx) [76]. 

It is tempting to imagine the temperature varying continuously from the Gibbons- 
Hawking temperature describing a de Sitter state — or from an approximate Gibbons- 
Hawking temperature associated with the approximate de Sitter state in our case — 
to the near Gibbons- Hawking temperature equivalent at late times in our asymptot- 
ically flat space, but this is perhaps unwarranted. At late times, the average number 
of particles created per mode from an early-time vacuum is well defined. This is not 
necessarily the case during inflation, when a choice must be made whether to make a 
measurement rapidly or slowly. If the measurement were made quickly, then by the 
time-energy uncertainty relationship, particles would be created through the act of 
measurement; if the measurement were made slowly, then the size of the scale factor 
would change appreciably during the measurement process, which could change the 
outcome [1]. Just as an observer accelerating through a Minkowski vacuum mea- 
sures particles [H], [771 EHl [79], [HO], [HI], so would a temperature-measuring device be 
excited in de Sitter space; however, unlike a thermal bath in flat spacetime, a mov- 
ing observer in de Sitter space would register no red-shifting in any direction. In 
fact, the authors of [82j find that for a massless, minimally-coupled scalar field in 
de Sitter space, no particles would be produced, and the associated effective tem- 
perature from these particles would be zero. With our composite scale factor, and 
using our late-time evaluation method alone, it is difficult to say whether particles are 
present during the exponential expansion, or whether they are created by the chang- 
ing Hubble constant at the end of inflation. It is likely that during the expansion, 
the long-wavelength modes that have exited the Hubble radius correspond to real, 
low-energy particles, while the high-frequency modes that have not left the Hubble 
radius correspond to virtual particles whose promotion to real particles depends upon 
the future evolution of the universe — such as our matching conditions — but to say 



conclusively whether particles exist during inflation would require a quantum field 
renormalization in curved spacetime, such as the adiabatic method given by [38] . 

By showing that the particle production of certain predicted particle species would 
cause conditions incompatible with observations in our universe, high-energy particle 
physics may be able to constrain the amount of reheating. Because we have shown 
how reheating — subject to ending conditions — is general to large- i^infi infiationary 
models, this can similarly be used to place mo del- dependent constraints on predictions 
for new particles, such as theorized supersymmetric partners of observed particles, 
under particular values of -ffmfi- In one such analysis [HS], if the gravitino G is the 
lightest supersymmetric particle (LSP), then this constrains the maximum reheating 
temperature to be less than 10^ GeV. If the G is not the LSP, and if its mass might 
be expected to be ~ 100 GeV, then the maximum reheating temperature may still be 
less than or about 10^ GeV [84j. Another example of a constraint on reheating is for 
the particle creation of scalar moduli, which may be present in supersymmetry and 
string theories: if the magnitude of the effective mass of the moduli field is less than 
i^infl, then the upper limit on the reheating temperature could be as low as 100 GeV 
[85] . This constraining works both ways. If evidence were found for the existence of 
such a reheating-constraining particle, this could eliminate those models of inflation 
that predict a large, nearly constant value of i^infi along with a rapid end to inflation. 
Those models that would be in agreement with such a low reheating temperature 
would be those with either a relatively small value of i/mfi, or those with a flnal 
period of inflation at which the inflationary potential has reached a near-minimum 
value, but at which it remains the dominant influence on the evolution of the scale 
factor, so that the initial high-energy particle production is greatly red-shifted and 
so that any unwanted relic particles are sufficiently attenuated such that they do not 
interfere with later early-universe processes, such as Big Bang Nucleosynthesis. 
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Part II: 

Binary System of 
Compact Masses 
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Chapter 4 

Unequal Mass Binary Solution in a 
Post-Minkowski Approximation 

In [8], Friedman and Uryu investigate a particular system of binary point masses 
that acquires a helical symmetry by taking the half-advanced plus half-retarded fields 
from the linearized Einstein equation. This time-invariant system in the co-rotating 
reference frame provides for an action at a distance theory, as has been previously 
discussed by [86l [871 EH]- It allows for a single action integral that depends on the 
dynamical variables and trajectories of each particle, without requiring a description 
of the force field acting on the particles. Such an action is called a Fokker action 
[8j. The Fokker action is not a true action, as the variation of the Fokker action 
integral depends on the boundary conditions and it involves integrals over each point 
mass's parameter time. When, however, a limit is taken after the variation of the 
Fokker action, in which its endpoints are taken at times of —oo and +00, the variation 
yields the correct equations of motion. The conserved energy and angular momentum 
associated with the Fokker action remain finite, even though energy and angular 
momentum of the field are infinite due to radiation from the system occurring over 
an infinite amount of time. 

In the post-Minkowski (PM) approximation, the metric is assumed to be fiat with 



small perturbations of the form gab = rjab + hab, where to linear order hab is the half- 
advanced plus half-retarded field of each particle. Unlike the post-Newtonian (PN) 
approximation, however, v/c '^1 need not be the case [89]. For the rest of this chapter 
we will use units of c = G = 1. Friedman and Uryu note that in zeroth order PM 
approximation T"''' = pu"-u^ = 0, and particles travel on flat space geodesies. A naive 
first order perturbation would then lead to 6T"'^ = 6 pu"'u°' + p6u°'u'' + pu°'6u^ = 6pu°'u°', 
which, because is the unperturbed straight-line motion, does not allow for bound 
orbits. In [8], this is avoided by considering a parameterized family of solutions to 
T°*(s) = p{s)u°'{s)u^{s) + p{s)[g°'''{s) + u°'{s)u^{s)] that corresponds to fiat space for 
s = 0. In a radiation gauge, 

- 2G1;) = D{hab - lvab)h = -levrTi;), (4.1) 

where h is the trace h°'a, the first-order stress-tensor is constructed from the first-order 
u", from the first-order p, and from the fiat-space metric. In the binary solution, to 
first order the motion of each mass is given by the linear field of the other, and the 
self-force serves only to renormalize the mass as a self-energy. Furthermore, Friedman 
and Uryii note of their post- Minkowski solution that it is correct to Newtonian order 
(OPN), the radiation field of the linearized metric is correct to 2.5PN, and a correction 
term to the equations of motion is necessary to have the orbits agree with the IPN 
solutions. For the case of the electromagnetic force, a specific example, in which 
the self force and radiation reaction are calculated, is given in greater detail in the 
Appendix [Al For the case of gravity, instead of photons the radiation takes the form 
of gravitational waves. The measurement of the energy loss due to this radiation in a 
particular binary system which contained a pulsar earned Hulse and Taylor a Nobel 
Prize in 1993. 

Although in linearized gravity non-linear terms are dropped that are of the same 
PN-order as linear terms that are kept, which means the next highest PM-order will 
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have terms of equal magnitude to those used at hnear PM-order, the post-Minkowski 
approximation may be helpful in evaluating solutions involving the full Einstein equa- 
tions in General Relativity that use helically symmetric initial data sets. Such initial 
conditions neglect the radial velocities associated with the radiation-reaction force, 
but a second-order post-Minkowski framework might lead to a better understanding 
of requirements for initial data in full-GR simulations. 

Fig. 14.11 shows the two point masses, m and m, with respective velocities v and 




Figure 4.1: Binary in Circular Motion. 

V. The radial parameters can be expressed as a = v/Vt and a = v/Q, where Q is the 
angular velocity shared by both point masses. Accounting for relativistic velocities. 
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the radial parameter is not equal to the 1 / (2-7?) times the circumference observed in 
the particle's co-moving frame. The position vectors are x"" = tt°' + aw"- and x'^ = 
it"' + dzu". The trajectory of m is tangent to the helical Killing vector k"" = t"" + Vlacf)"', 
and the trajectory of rh is tangent to the helical Killing vector k"" = i"" + fia0", where 
7 = dt/dr. 

In Fig. 14.21 the Law of Cosines relates t"^ = + — 2aacos(7r — f), or (ip/Q)'^ = 

m 




Figure 4.2: Retarded Angle (p. 

{v/VtY + {v/VlY + 2{v /VL){y /Vt) cos((y9), so that the retarded angle, which is equal in 
magnitude to the angle associated with the advanced position, is given by the positive 
root of the transcendental equation ipP' = v"^ + v"^ + 2vvcosip. 
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From two types of Fokker action, a parametrization invariant action with a post- 
Newtonian correction and an affinely parametrized action, the equations of motion 
and expressions for conserved energy and angular momentum are derived following the 
variational calculation of Ref. [90] . In the Affine case, we parameterize the trajectories 
using the perturbed flat-space metric as {rjah+hab)x'^i^ = —1 and {rjab+hab)x"'x^ = — 1, 
where the dots represent derivatives with respect to the parameter times of x(r) and 



x{f). This leads to 7 = (1 - t;^ - habk^'k^)'^/'^ and 7 = (1 - t;^ - habk^'k^)'^/'^ . In the 



parameter-invariant case, we parameterize the trajectories using the fiat-space metric 
as rjahi'^x^ = —1 and r]abx"'x'' = —1. This leads to 7 = (1— f ^)^^/^ and 7 = (1— -u^)^^/^. 
The affine parameterization is characterized by the following: the parameter times of 
geodesies are the proper times of the perturbed metric; the PM-form of the geodesic 
equation {r]ab + hab)x^ + Cabci^'x" applies, where Cabc = 0-/'^)iy bhac + ^ Ma-^ ahc)] 
and, finally, the 4-velocity is orthogonal to the 4-acceleration, or U"'VaU'' = 0, that 
is the particles travel along geodesies. The linear post- Minkowski approximation is 
not at this point accurate to IPN order, but Friedman and Uryii give two different 
adjustments to the parametrization-invariant case: the simplest correction consistent 
with IPN (called PN where confusion will not arise) and a correction that is both 
parameterization-invariant and special-relativistically covariant (SPN), where results 



are given in |8j for the deDonder gauge. They show also that for both of the Fokker 



actions the form of the first law of thermodynamics dE = VtdL holds, and this law 
can be used to check for the presence of an Innermost Stable Circular Orbit (ISCO). 

We find a solution describing a helically symmetric circular orbit in the post- 
Minkowski approximation (with post-Newtonian corrections) that is analogous to the 
circular solution of two charges obtained by Schild for the electromagnetic interaction 
|91j . In [6j we report results supplementing those of [8j: numerically computed solu- 
tion sequences for unequal mass particles, and analytic formulas in the extreme mass 
ratio limit. The latter results agree with the first post-Newtonian (IPN) formulas; 
hence a consistency of our model is confirmed in this limit. 
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We present a set of formulas governing the helically symmetric circular orbits of 
two point particles, {m, v} and {m, v}, and derive analytic expressions in the extreme 
mass ratio limit g = m/m — 0. The set of algebraic equations is solved numerically 
for a fixed binary separation to specify each circular orbit. The result for the unequal 
mass binary orbit is presented in Sec. 14.11 

We compute the solution to the equation of motion numerically for three mass 
ratios: q=1.0, q=0.1, and q=0.001. We solve the equation of motion for each mass 
ratio in the PM+PN model, the PM+SPN model, and the affine model. We also 
calculate the solution for the g ^ limit analytically in each of the three models, 
for which see Sec. 14. 2[ Whenever the analytical solution is plotted along with the 
q=0.001 numerical solution, the two lines overlap in the graphs given here. 



We discuss solutions to the post- Minkowski approximation in the case of parametrization- 
invariant plus IPN correction terms, and then we discuss solutions in the affine case. 
For the analytical solution in the g — > limit, see Sec. 14.21 

Parameter Invariant Circular Solution 



We first list the result from |8j for the parametrization invariant model with IPN 
correction terms. After integration, the equations of motion for particles m and m 
are written in terms of the velocities, v and v, of particles m and m, which are related 
to the orbital radius hj a = v/Q and d = v/Q, through the equations: 

— m-y'^vQ = —mThj'^'yQ'^[F{ip,v,v) + {m + Th)QFi{ip,v,v,'y,j)], (4.2) 
—fh'j'^v^l = —mmj'j'^^l^[F{(f,v,v) + {m + Th)^lFi{(p,v,v,'y,'y)]. (4.3) 



4.1 Numerical Solutions 
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As shown below, {yj, t>, t), 7, 7} are not independent. The functions F{ip,v,v) = 
F{ip,v,v) are the post- Minkowski terms, while -u, f , 7, 7) = f , 7, 7) 

is either of two alternative IPN correction terms that agree at IPN order: Fi = 
-Fpn(v^, V, V, 7, 7) derived from a non-relativistic correction, or Fi = Fspn(v^, v, v, 7, 7) 
derived from a special relativistically invariant correction. These are 



F{ip,v,v) = —4 



Fpn(</?,w,w,7,7) 
^spn(¥',w,w,7,7) 



{if + vv sin 



(1 + vv cos ip)v 



X {ifi cosif — f ^ siny?) + ^'^{^ -~ '^^)(v^ + "^"^ sin (p) 

— -[vsin(f{(f + vvsinif) + (1 + vvcos!f){v + vcos(f) 



V ^ 

-{ip + vvsm^Y](^{^,v,v) \, 

1 — ) 



^2^2^^ _|_ ^^3 
1 



(v? + vv sin ipY 



1 + -7^f (t; + v) 



-7 + 

}■ 



w sm 



(77)^^^ (V5 + wsinv?)^ l4 (/} + OTsin(/} 

(1 + vv cos (f + w cos 99) 



(4.4) 
(4.5) 



(4.6) 



The function $(y9,f,f) is defined by 



(1 + wcosv?)^ - \{l - v'^){l - v^) 
if + vv sin ip 



(4.7) 



For the parametrization invariant models, 7 and 7 are derived from a flat-space 
normalization of the four-velocity. 



7 = (l-^')~^: 



7 = (1 -t;2)-2. 



(4i 



The retarded angle ip is the positive root of <y9^ = f ^ + -u^ + 2ft; cos ^9. 

In Fig. 14.31 the angular velocity, in dimensionless form fiM, is plotted against 
the velocity of the lighter particle for 3 mass ratios and the g limit for the 
parametrization invariant model with PN correction. Curves of the analytic solution 
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for g ^ and that of q = 0.001 overlap each other in the plot. The inflection 
displayed in the logarithmic plot changes near the cutoff velocity for the small mass 
ratio cases. In Fig. 14.41 the angular velocity, in dimensionless form QM, is plotted 



a 




Figure 4.4: Parametrization-Invariant (SPN) Omega versus Velocity. 

against the velocity of the lighter particle for 3 mass ratios and the q ^ limit for 
the parametrization invariant model with SPN correction. Curves of the analytic 
solution for g — > and that of g = 0.001 overlap each other in the plot. The inflection 
displayed in the logarithmic plot changes near the cutoff velocity for the small mass 
ratio cases. 

In this notation Kepler's Law, (Tperiod)^ = 47r^a^/m, may be written as {rfiQ)'^ = 
{fh/aY [6]. In this form it may be compared with Eq. (14. 2p . when written as. 



{^mf = (^^) ^ {v^^ [F+{m + m)n Fi] } 



(4.9) 
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To see how this post-Minkowski approximation is related to Newtonian gravity in the 
non-relativistic {v « 1) hmit, see Sec. I4.2[ 

Parameter Invariant Energy and Angular Momentum 

In Figs. 14.51 (PN) and 14.61 (SPN), the unit energy of the lighter particle, in dimen- 
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Figure 4.5: Parametrization-Invariant (PN) Energy versus Omega. 



sionless form E /m, where E = E — rh, is plotted against QM. In the limit that the 
particle approaches becoming unbound, v ^ 0, the unit energy of the lighter mass 
approaches 1, its rest mass energy. As it becomes more tightly bound, its energy 
decreases below the rest mass energy it would have in flat space. 

The conserved energy and angular momentum for the parametrization invariant 
model are written 

E = EpM + ei, and L = Lpm + ^i, (4.10) 
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where -EpM and LpM are the post-Minkowski terms 
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^PM = - + - 4.11 

7 7 



^PM 



2mm77 f, f), (4-12) 



and ei and ii are the parametrization invariant IPN corrections ei = epN and ii = ^pn, 
or those of the special relativistically invariant model ei = cspn and ii = £spn given 
by 



epN = ^^^^PN, (4.13) 
esPN = ^^^4pn, (4.14) 



^PN 



-SPN 



mm{m + ffi)VL 
77(f + vY ' 
mm{m + m)Q 



(^-y)3/2 [ip + vvsimpy 



(4.15) 
(4.16) 



In Figs. 14.71 (PN) and 14.81 (SPN); angular momentum, in dimensionless form 
J/ {mM), where M is the total mass of both particles and m is the mass of the lighter 
particle having velocity, v; is plotted against QM for 3 mass ratios and the g 
limit. There is no ISCO, but at the maximum value of v for each mass ratio, beyond 
which there are no further solutions, there is an Innermost Circular Orbit (ICO). 
One possible explanation for the termination of solutions can be found by looking at 
Eq. (14. 9p . which may be written as a quadratic equation in terms of {rfiQ). Beyond 
the maximum value of v, the solutions for (mf2) become imaginary. 
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AfRne Circular Solution 



In Fig. 14.91 The angular velocity, in dimensionless form f2M, where M = (m + m) 

10° 
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Figure 4.9: Full Affine Case of Omega versus Velocity. 

is plotted against the velocity of the lighter particle for 4 mass ratios in the affine 
model. The behavior for solutions existing beyond fisco? the velocity at which the 
minimum energy and angular momentum occur, is most prominently displayed for 
the q=1.0 case. In the q = 1 case the ISCO occurs at v ~0.184. Fig. 14.101 shows the 
same data as Fig. 14.91 but it only shows solutions up to the ISCO. 

For the affinely parametrized post- Minkowski model, analogous forms of Eqs. (14.21) 
and (14. 3p are written 



(4.17) 
(4.18) 
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where the function F^{ip,v,v) = F^{ip,v,v) is defined as 



-4 



{(f + vv sin (fY 



(1 + vv cos ^)v{ip cos V9 — sin Lp) 



+ -f (1 — f^)(v5 + ff siny^) — - [vsiinp{ip + ff sin(^) 
+ (1 + vv cos V9)(t> + cos (v^, f , "w) [• 



(4.19) 



For the affinely parametrized world fine, 7 and 7 satisfy 



7 (1 - V ) + 4m7 7^] V, v) 
-7^(1 - -u^) + 4m77^fi f , t;) 



(4.20) 
(4.21) 



AfRne Energy and Angular Momentum 

In Fig. 14.111 the lighter particle's unit energy per mass, in dimensionless form E/ra, 
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Figure 4.11: Full Affine Case of Energy versus Omega. 
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where E = E — m, is plotted against QM for the affinely-parameterized case. In 
Fig. 14.121 angular momentum, in dimensionless form J/{mM), where M is the total 

10^ 



10 












0.05 0.1 



0.15 



Figure 4.12: Full Affine Case of Angular Momentum versus Omega. 

mass of both particles and m is the mass of the lighter particle having velocity, v; is 
plotted against QM for 3 mass ratios and the g — > limit in the affine model. Minima 
of each curve corresponds to the ISCO. 

The conserved energy and angular momentum for the affinely parametrized model 
are written 



mm _ 
E = h — + 4mm77iZ f , f ), 

7 7 

= m7(l — f^) + 777,7(1 — v'^) — Amrh'j'jil $(</?, v, v), 

L = 277777777 <!>((/?, f, -u), 

(1 + wcosv?)^ - i(l - f^)(l - v'^) 



277777777 



(p + vv sin ip 



(4.22) 



(4.23) 
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where the form of ^{<f, v, v) is the same as that of the parametrization invariant model 
( 14. 7p . Using Eq. fl4.20p and fl4.2ip . the energy can be rewritten 



hi = 1 — m7 1 — f H 1 — m7 1 — f . 

2 7 2 ^ 2 7 2 ' 



(4.24) 



This can be compared with Eq. (14. lip , noting the different definitions of 7 in the 
parametrization-invariant and affine models. 

Figs. 14. 131 and 14. lil show the same data as Figs. 14. 1 ll and 14. 121 respectively, except 
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Figure 4.13: Truncated Afiine Case of Energy versus Omega. 



that only the solutions where < f < v-^sco are plotted. 

In the affine case, for any mass ratio q G [0, 1] we find a simultaneous minima in 
the energy and angular momentum which corresponds to the ISCO. The values of the 
normalized angular velocity, angular momentum, and energy that occur at the ISCO 
in the afiine model vary monotonically from q=l to 1=0. With q ranging from 1 to 
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Figure 4.14: Truncated AfRne Case of Angular Momentum versus Omega. 
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0, fiM decreases from ^ 0.0521 to ^ 0.0440, L/{mM) increases from ^ 2.0558 to 
^ 4.2617, and E/m decreases from ^ 0.9775 to ~ 0.9593. 



Numerical Solutions for Unequal Mass Circular Orbit 

A circular solution is calculated from algebraic equations given in Eqs. (14.21) and (14.31) 
for the parametrization invariant model or Eqs. (I4.17P and (14.181) for the afinely 
parametrized model. One method of solving for a fixed ratio q = m/rh is — (1) 
assume a ratio of velocities v/v and determine the corresponding mass ratio from the 
equations of motion, then (2) change the velocity ratio to adjust the value of the mass 
ratio to a fixed value (using the bisection method, for example). The mass ratio, q, can 
be determined by multiplying both sides of either Eqs. (14. 2 p and (14.30 or Eqs. (14.170 
and (14.181) by v'^'^/{v)'^) and then dividing both sides of the first equation listed in 
either of these pairs of equations with the second equation to yield an expresion for 
the mass ratio q. 

Another method involves solving the relationship MVt[q^ v, v]=M Q[{l/q),v, v] by 
varying the parameters v and v. For a description of an efficient means of sampling 
a parameter space and fine-tuning the optimal result, see Sec. 15.2. 1[ 

4.2 Analytical Formulas for Extreme Mass Ratio 

For the Extreme Mass Ratio, the mass of the lighter particle is negligible relative 
to that of the more massive particle. This would be appropriate for a test mass 
orbiting in the spherically symmetric gravitation field of a much more massive ob- 
ject. Whereas the signal from a merger of identical black holes, each with a mass 
on the order of one solar mass, would fall within the sensitivity of LlGO's (Laser In- 
terferometer Gravity-wave Observatory) frequency band; the inspiral of a solar-mass 
black hole into a billion-solar-mass black hole, such as those predicted to be at the 
centers of many galaxies, would fall into the most sensitive frequency band of LISA 
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(Laser Interferometer Space Antenna) [92j. The Extreme Mass Ratio is an excellent 
approximation to this latter scenario of a mass ratio on the order of 10~^. 

Extreme Mass Ratio Limit 

The extreme mass ratio limit g = m/m^Ois identical to the limit v with Q 
fixed. In the limit w — > 0, we may assume that v and rh remain finite. Consequently, 
we have 7 — 1, <y9 — f , and rh M, where M = m + m is the total mass. With v 
and Q regarded as independent variables, Eq. (14.21) is a quadratic equation for QM, 
whose q = form is 

Fi {VLMf + F {VlM) - t; = 0, (4.25) 

with physical solution 

nM=^ [~F + ^F^ + AFiv^ . (4.26) 

The functions F (the post- Minkowski term), Fi = Fp-^ (the simplest IPN correction), 
and Fi = -Fspn (the special-relativistically covariant IPN correction) for q = become 

1 - S7;2 

nv>,v,v) = ,y (4.27) 

Fp^{^,v,v,j,j) = -^(1-^^')' (4-28) 

(I _ ^;2U/4 / 1 \ 



where $ has the form 



H^,v,v) = (4.30) 



Without the IPN correction, the parametrization-invariant post-Minkowski model is 
given by setting Fi = 0, and therefore QM = v/F. In the g — > limit, this is 

QM = (4.31) 
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Parameter Invariant Solution Sequence in g ^ Limit 

In the g — > limit, the conserved energy and angular momentum normalized by 
the mass remain finite. Subtracting the mass of the heavier particle from the post- 
Minkowski energy, EpM = EpM — rh, and taking the limit v with rh —>■ M, we 
have 



(4.32) 
(4.33) 

(4.34) 
(4.35) 

QM, (4.36) 



-E'PM 




m 




Lpu 




mM 


77(1-^2)1/2' 


epN 
m 


2mM 


esPN 
m 


= ^^^™f]M, 
2mM 




(1_ ^2)1/2 


mM 




^SPN 


(1 - i;2)3/4 


mM 


y2 



QM. (4.37) 



In [8], it is shown that the first law of thermodynamics that relates the changes in 
the conserved energy and the angular momentum, dE = QdL, is satisfied by binary 
solutions derived from the parametrization invariant Fokker action. This relation 
is used to cross check both the analytic formula in the g — > limit as well as the 
numerical solutions obtained in Sec. 14.11 by calculating dE/dv = VtdL/dv, where 
E = EpM + ei. 

In the parametrization invariant post- Minkowski model without a IPN correc- 
tion, the normalized angular velocity, f2M, is defined in an interval < f < 1/ for 
g — 0, and f2M becomes infinite at f = 1/^/3- With the IPN correction Fi = -FpN, 
the range of finite QM is approximately < v < 0.361598, and with the special rela- 
tivistic invariant IPN correction Fi = -Fspn, it is < < 0.36166. Newtonian point 
particles have no innermost stable circular orbit (ISCO), but adding a IPN correc- 
tion to the Newtonian orbit recovers the ISCO that is present in the exact theory of 
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general relativity. In the post-Minkowski framework, we find that the existence of an 
ISCO depends on our choice among actions that are equivalent to first post-Minkowski 
order. In particular, we find that the parametrization-invariant action leads to se- 
quences with no ISCO even when IPN terms are included. This is plausibly due 
to the fact that the sequences associated with the parametrization-invariant action 
terminate before reaching the angular velocity of an ISCO. For the IPN formalism 
given in [93], an Extreme Mass Ratio ISCO occurs at the unrealistically high value 
of QM = 0.544. The 2PN and 3PN values for the g = ISCO are QM = 0.124 and 
0.0867, respectively [93J. Below we show that sequences associated with the affinely 
parametrized action do have an ISCO; however, the g ^ ISCO of the affine case 
occurs at an unrealistically small value of QM. 

In Eq. f l4.26p . an expansion of QM in the small v limit becomes QM = + 3v^ + 
0(f for both PN and SPN models, and this is inverted to write v in terms of small 



Substituting this into the energy and angular momentum formulas, the leading two 
terms agree with the post-Newtonian formulas (see e.g. [93]) up to the IPN order for 
the extreme mass ratio q ^ 0, 



VtM as 



V = {QM) 



(4.38) 



E 




(4.39) 



m 



L 1 



l + ^(fiM)2/3 + 0((nM)^/=^) 



(4.40) 



mM (fiM)V3 



AfRne Solution Sequence in g — > Limit 



Eq.f l4.17l) implies 



-1 



(4.41) 



where 7 is evaluated from Eqs. (14.201) and (14.211) . 
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In the limit of g — » (or more directly f ^ 0), 



l-v 



2 



F^'i^.v.v) = (4.42) 



From Eq. fjiTTTl) and fj420l) . we have 



1 2 \ 1/2 



while in Eq. (14.211) . taking t> — >• and m — >• yields 7 — >• 1. As a result we have in 
the extreme mass ratio, 

3 

= (4.44) 
1 — f ^ 

In the g — limit, the energy without the rest mass of the heavier particle, 
E = E — rfi, and the angular momentum become 

^ (l-3t;2) 

(4.45) 



m [(l-w2)(l_4^;2_^4)]l/2' 

/2 

(4.46) 



L 1 + i; V 1 - t;2 ^ 



The first law 6E = Q6L is also satisfied for the affinely parametrized model, and 
hence one can cross check formulas in the q limit using the relation dE/dv = 
QdL/dv. Although the lighter particle's normalized angular velocity, QM, is finite 
in an interval v G [0,1), the redshift factor 7 as well as conserved quantities E 



and L become infinite at v = vv5 — 2 ^ 0.485868, which corresponds to QM 
{V5 - 2)^^V(3 - v^) ~ 0.150142. 



In this interval, v G [0, v v5 — 2), the energy and angular momentum have a 
simultaneous minima at v = a/(1 + 2^/^ — 2^/^)/3 ~ 0.339136, which corresponds to 

(1 + 2^/3 _25/3^ 3/2 

QM = ^ ^ 0.0440743. (4.47) 

2v^(l -2V3 + 22/3) ^ ^ 
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The Schwarzschild ISCO occurs at fiM = G'-^^^^-'v^ ^ 0.096. In terms of this exact 
solution, the ISCO of the affine parametrization has an error of 54%, whereas the 
ISCO of the IPN approximation given by [93j has an error of 465%. 




Radial Parameter in g — Limit 

With the definition a = f we can write a/M = v/{MQ), where M = m + m in the 
g ^ limit is just M = m. Then we insert into Eqs. (14.311) (OPN parametrization- 
invariant without IPN correction term), fl4.26p (PN and SPN cases), or fl4.44p (OPN 
affine case) the maximum velocity (parametrization-invariant) or the ISCO velocity 
(affine). These cutoff velocities are— 1/^3 (PM), 0.361598 (PN), 0.36166 (SPN), or 
0.485868 (Affine). This leads to a minimum radial parameter for a circular orbit. In 
units of M~^, these minimum radial parameters are as follows: (PM), 2.67 (PN and 
SPN), and 3.24 (Affine). Note that the ISCO of the affine case occurs atv = 0.339136, 
which corresponds to a/M ^ 7.69. 

The parametrization-invariant model without a IPN correction has no minimum 
radius, which is the same as Newtonian gravity. The affine case has an ISCO on the 
order of the 6M that is predicted by the full theory of general relativity. In the case 
of the PM+IPN correction term model, and in the affine case without a correction, 
the minimum radial parameter occurs on the order of the 2M event horizon for a 
Schwarzschild black hole. 

IPN Energy and Angular Momentum 

For comparison with our post- Minkowski analysis, we list the IPN equations of j93j . 
where Blanchet's u is our g(l + g)~^. In our notation, 




M 



E 



1 q 
2(g + i) 



-(M^])2/3 1 



( 



3 1 q 

4 ^ 12(1 + g)2 



) 



(4.48) 
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or, 



and 



- = 1 - -^{unf^ 

m 2 g + 1 



1-1?+^ ' 



M2 (1 + g)2 (M^])V3 



4 12(l + g)2 



. ,31 q 

1 + + 77- 



{Mn) 



2/3 



2 6(l + g)2 



2/3 



or 



1 



mM 



q (Mfi)V3 



2/3 



(4.49) 



(4.50) 



(4.51) 



2 ' Q{l + q)\ 

Thus we show exphcitly in the hmit g ^ that in the parametrization-invariant case 
with a first-order post-Newtonian correction term, the energy and angular momentum 
Eqs. f lCT]) and flOO]) agree with Eqs. fICT]) and fOT]) . 

Agreement between the energy and angular momentum formulas of the IPN cir- 
cular solution, and those of the parametrization invariant post-Minkowski model with 
post-Newtonian correction, is exhibited explicitly for the extreme mass ratio limit. 
For an arbitrary mass ratio one needs to expand the retarded angle ip to the next or- 
der in the velocities, v and a.s (p ^ {v + v){l — vv/2), and the rest of the calculation 
closely parallels that of the g = case. 
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Part III: 

Production and Decay 
of Small Black Holes 
at the TeV Scale 



Chapter 5 

TeV-Scale Black Hole Production 
at the South Pole 

We discuss the possibility of observing TeV-scale black holes produced at the 
IceCube Neutrino Telescope [7]. After giving a brief summary of the IceCube exper- 
iment, we explain what TeV-scale black holes are. We then examine a gravitational 
interaction between a neutrino and a nucleon. Because a nucleon is not a point 
particle, we rely on the parton model, which describes the nucleon as a collection 
of quarks and gluons. Following this, we describe our method for modeling Parton 
Distribution Functions (PDFs), we evaluate the cross section for the interaction of 
neutrino+nucleon— i>black hole, and then we calculate IceCube's detection sensitivity 
for observing TeV-scale black holes. 

5.1 IceCube Neutrino Telescope 

In the Standard Model a neutrino can interact with a nucleon through both charge 
current (CC) interactions and neutral current (NC) interactions [SI [95], [96]. In a 
CC interaction, a neutrino (anti-neutrino) interacts with a quark to become a lepton 
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(anti-lepton), conserving electron-, muon-, and tau-lepton number. In this interaction 
a particle is exchanged with a down-quark (up-quark), which becomes an 

up-quark (down-quark). In a NC interaction, a neutrino exchanges a Z° with a quark 
and neither the neutrino nor the quark changes flavor. 

Because neutrinos experience only gravity and the weak force, they may travel as- 
tronomical distances without interactions. Thus, they preserve information about the 
environment in which they were produced. The corollary to this is that a sufficiently 
large detector must be used to observe these cosmic neutrinos here on Earth. 

The IceCube Neutrino Telescope is composed of approximately one cubic kilometer 
of Antarctic ice ranging from 1400 meters in depth to 2400 meters in depth below the 
surface near the Amundsen-Scott Station located at the geographic South pole [97] . 
IceCube is already taking data, and it is scheduled to be fully operational by 2009- 
2010. At that time, it will consist of 80 strings, each a kilometer long, of 60 evenly 
spaced PhotoMultiplier Tubes (PMT) each, for a total of 4800 PMT. The strings 
are 125 meters apart, and each interior string will be surrounded by six equidistant 
neighbors [98| . 




When high energy charged particles move faster than the local speed of light 
through the ultraclear Antarctic ice, in which the absorption length of the relevant 
wavelengths is greater than 100 meters [99], they emit Cherenkov radiation. This 
radiation, within a range that includes visible light and some UV light, can be detected 
by the PMT used in IceCube, and the time at which this happens — including the time 
for the signal to register — can be recorded within an accuracy of a few nanoseconds 
|100] . The paths of these charged particles may be dominated by jets from a high 
energy muon or tau. They may also be diffused throughout a shower. With sufficient 
data, these paths can be used to reconstruct the particle interactions that have taken 
place. This requires the energy of the incident neutrino to be greater than 100 GeV. 
In the case of a series of interactions caused by a single incident particle, the total 
energy — provided that it is contained within the volume of IceCube and is less than 




10^° GeV so that it does not saturate the detector — can be measured 

When measuring neutrino interactions, one must contend with a background event 
rate of charged particles, such as muons produced by cosmic rays hitting the atmo- 
sphere |101j . Examining upward going tracks, or particles that have passed through 
a significant fraction of the Earth, effectively restricts the progenitor particle of an 
interaction to a neutrino, which, because it is only weakly interacting, is able to eas- 
ily penetrate the Earth, whereas charged particles are not. A horizontally traveling 
neutrino passing through the center of IceCube travels through 150 kilometers of the 
Earth |102j . There is also a background trigger rate for IceCube's PMT of less than 
one kilohertz [97[[99]. That is, in the absence of a signal a PMT will discharge on av- 
erage no more than once every millisecond. This is not a problem, because the transit 
time across IceCube for those particles that produce Cherenkov radiation is on the 
order of a few microseconds and the PMT recording time is accurate to within a few 
nanoseconds. Although the volume of IceCube is one cubic kilometer, the effective 
volume for detecting neutrinos is larger, because muons may be produced outside the 
IceCube volume and still travel inside to be measured [103] . 

Another useful veto is the IceTop surface array of 160 Cherenkov detectors of 2.7 
meter diameter tanks of ice spread out over one square kilometer of area [104j . IceTop 
helps reject background events and is also useful for calibration. 

Amanda, the prototype of IceCube that proved the viability of detecting neutrinos 
in polar ice caps, is still running. Because its volume overlaps with the volume of 
IceCube, it can either contribute to IceCube's sensitivity, or it can serve as a check 
on IceCube detections, depending upon whether data from the two experiments is 
examined collectively or independently [lOOj . 
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5.2 Black Hole Production in Higher Dimensions 

In the standard model (SM), gravity is by far the weakest of the four fundamental 
forces. It has been theorized that this weakness is due to the presence of extra 
dimensions beyond the 4 familiar dimensions of our spacetime [IDB]- If gravitons 
propagated into the extra dimensions while SM fields were confined to our brane of 
3+1 dimensions, then gravity thus diluted would appear much weaker than the other 
forces. In this case, gravity might become much stronger at small distances than a 
4-dimensional theory would predict. 

We will investigate the possibility that the distance at which gravity and the 
electromagnetic force have the same strength is at ~ 10^^^ m, the distance at which 
the electromagnetic and weak forces unify as the electro-weak force. This would 
mean that for the small distances at which gravity matched the electro-weak force in 
strength, there would be a fundamental D-dimensional Planck mass of about 1 TeV, 
in which case our 4-dimensional Planck mass would just be an effective Planck mass 
over macroscopic dimensions. 

The strength of TeV-scale gravity at small distances could potentially make it 
easier for interacting particles to form black holes. This can be qualitatively un- 
derstood via Gauss's Law \106\ I107j . The surface area of a sphere in D dimen- 
sions, where there is 1 time dimension and D — 1 spatial dimensions, is propor- 
tional to r^^^. The magnitude of a D-dimensional Newtonian gravitational force 
acting between two masses would be proportional to Mim2G£,, where G^i is the 
D-dimensional gravitational constant. Spread evenly over the surface area of a 
sphere, this force would be proportional to Mim2GDr'^~^ ■ The value of the po- 
tential energy at a separation r between the two masses would be proportional to 
Mim2GDT^~^ ■ Taking Mi as the primary mass and m2 as a test mass, then using a 
non-relativistic argument to relate the maximum kinetic energy of a test mass moving 
near the speed of light (mc^/2) to the potential energy, places an event horizon at 
r oc (MiG£)/c^)(^/[^^^l). The dimensionality of Go is length^"^ mass~^ time~^. The 



D-dimensional Planck mass, Md, is then proportional to {h^ ^c^ ^Gr, ^)(^/[^ ^D. 
In units of /i = c = 1, then Gd oc M|,~^, and the event horizon would be r oc 
(MiG^)(i/[^-3]) oc (MiM2-^)(i/P-3]) oc (1/Mz5)(Mi/M^)(i/[^-31). 

It has been suggested that the Large Hadron Collider (LHC) could easily produce 
such black holes in this scenario [108[ I109^ 1110] . If the LHC would be powerful enough 
to detect this sort of black hole interaction, then cosmic rays would also be energetic 
enough to produce this interaction. In particular, we will discuss the possibility that 
neutrinos produce black holes in the ice of the south pole and can be detected by the 
IceCube Neutrino Telescope. 

To model the gravitational interaction between a neutrino, which is a point par- 
ticle, and a nucleon, which is an object of finite extent and which has an internal 
structure attributable to constituent point particles, we turn to Parton Distribution 
Functions. 

5.2.1 Modeling Parton Distribution Functions 

In high energy interactions between a neutrino and a nucleon, the neutrino interacts 
primarily with a single parton, a quark or a gluon. For these collisions, the proton 
and neutron are not just an up-up-down and an up-down-down, but are composed of 
these and other, virtual particles that are continually created and annihilated through 
the time-energy uncertainty relationship. 

Similarly, for low energy interactions, a nucleon acts as a single particle of rest 
mass energy in its rest frame. For high energy interactions between a neutrino 
and a parton, in the nucleon's rest frame the parton will have have some fraction 
of the total energy rest-mass of the nucleon. This fraction is denoted as x, where x 
ranges from to 1, or from none of the nucleon's total energy to all of it [111] . 

A Parton Distribution Function (PDF) describes the probability of finding a given 
parton — up (m), anti-up (m), down (rf), anti-down (J), strange (s), anti-strange (s), 
charm (c), anti-charm (c), bottom (6), anti-bottom (6), or gluon {g) — with a fraction 
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X of the total rest energy of m^v- The contributions from the PDFs for the super- 
massive top and anti-top within the nucleons at rest are neghgible, and we neglect 
them. Thus, the probability that an ith species of parton exists with a fractional 
energy between xi and X2 is 



X2 



V= / h{x,Q)dx. {h.l] 



The variable Q is the momentum transfer, where we choose Q ~ ^ [112] . and the 
PDFs are somewhat insensitive to changing Q |113] . We thus use 



g = min{r;\ 10 TeV}. (5.2) 



The PDFs cannot be calculated analytically from first principles in the Standard 
Model. They must be fitted to experimental data. We use the CTEQ6D PDFs [114j . 
The largest uncertainty in the PDFs exists for large-x gluons, where /(x, Q)^\aon may 
be off by more than a factor of 2 [115] . At small x, where the PDFs are much 
more certain, the gluon quickly comes to dominate the neutrino-parton interactions 
through its high probability of being available for a collision. 

Fig. 15. II plots log(a;) versus log(a;/j(x, Q)) for a representative quark, the up, and 
for a gluon, both of which for the relatively low Q of 10 GeV, or = 100 (GeV)^ 
[116] . Fig. 15.21 plots log(x) versus \og{x fi{x,Q)) for a representative quark, the up, 
and for a gluon, both of which for the relatively high Q of 10 TeV, or = 100,000,000 
(GeV)^ [IISJ. The variable changes by six orders of magnitude between these two 
cases, but the PDFs shown only change by about an order of magnitude. For x less 
than about 10^'^, the graphs of the PDFs are nearly linear in these log-log plots. For 
this reason we use different models of these PDFs for small x and large x. We also 
use a different modeling of PDFs for the ranges Q > 10 TeV, 10 TeV > Q > 1 TeV, 
1 TeV > g > 100 GeV, 100 GeV > Q > 10 GeV, 10 GeV > Q > 1 GeV, and 1 
GeV > Q. These different regimes of PDFs lead to the almost imperceptible bulge 
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Figure 5.1: Parton Distribution Functions: Lower Momentum Transfer. 



between = 10^° GeV and = 10" GeV in Fig. 

For these different regions of x and Q, we make use of simple approximations to 
tlie PDFs by fitting tlie CTEQ6D data to a form of 

Mx,Q)=Ax^, (5.3) 

where, for example, in the small-x and large-Q regime, n ~ —1.4 for all the partons 
and A ~ 0.2 for quarks and A ~ 3.5 for gluons. Because we use different PDFs for 
the different regions, this form of Ax"' is a good approximation that is simple to use 
when we integrate the cross section for a black hole interaction. 

What follows is a brief description of our numerical method. To accurately fit both 
the variables A and n, we refine our best guess and also sample the two dimensional 
parameter space. At a fixed value of Q and given a two dimensional array relating 
fi{x) to X, we start with a reasonable guess for the variables A and n and a reasonable 
value for our step variable. At each iteration, we compare our previous lowest result 
for the sum of the squares of the difference between the given data points and Ax"' for 
all the points in the array with new values of the variables A and n. We try altering 
our current best values of A and n by increasing or decreasing one or the other or 
both in tandem or opposition for a total of 8 different combinations. If one of these 
combinations results in a better fit, then we store these new values of A and n as our 
new current best values, and we retain the new sum of the squares of the difference 
between the given data points and the new Ax" as the new best target, and then we 
repeat the eight combinations. If we do not find a better fit, then we decrease the size 
of our step variable and repeat the above algorithm. If we reach a sufficiently small 
step variable, we do not yet give up: there are local quasi-minima in the parameter 
space that are not good fits, such as A = and n <^ —1. We instead pick a new 
value of our step variable that is large enough to jump to unexplored, and potentially 
rewarding, areas of the parameter space. To prevent getting stuck with the same 



poor choice iteration after iteration, we choose a random number between and a 
reasonable maximum for our new step variable. At this point it is better to choose 
a step variable that is too big, or big enough to jump to a new region of parameter 
space, than too small to be effective. The subsequent shrinking of the step variable 
will take care of any initial excess. After a set number of loops of this entire process 
where the best fit remains unchanged, we take the resulting best values of A and x 
as our best fit. 

A simpler method than the one just described involves simultaneously changing 
both variables (or more, if a problem requires sampling a higher dimensional param- 
eter space) with different random steps weighted towards zero. This method was 
not used for the project described here, but I have used it elsewhere with success. 
Each parameter is changed by a different step value, and each step value involves 
two layers of randomness. The first layer determines the order of magnitude of the 
step, with a sizable probability it will be insignificantly small or zero, and the second 
layer determines the coefficient and sign associated with the order of magnitude. This 
simpler method effectively encapsulates the whole of the important parts of the above 
method in very few lines of code and is faster at searching for and homing in on the 
best solution. 

To check our results, we plot our best Ax"" against the CTEQ6D PDF data to 
ensure our data is a good fit. For an appropriately nearby value of Q, the benefit to 
using a simple function over an array of data is that we can calculate fi{x) for any 
given X, and we avoid both having to maintain in our program memory a complicated 
series of PDF arrays and having to interpolate between data points in these arrays. 
Our fits are an excellent approximation to the PDFs being used and are certainly well 
within the uncertainty of the PDFs, themselves. 

For examples of these PDF fits for the u-quark, see Figs. 15.31 and 15. 4[ The first of 
these shows why the PDF fits are broken up across the range of the parton momentum 
fraction, x, and the resulting excellent fit. The second of these figures shows how the 
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Figure 5.3: An Accurate Fit for a Small Range of Data. 
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PDF fit can drift from the exact PDF wlien overextending the fit over too large a 
range of x. This may also result from overextending a particular fit over too wide 
a range of Q^. In both of these figures a value of = 10^ (GeV)^ is used for the 
CTEQ6D PDFs. Our fit in Fig. is given by Eq. with A = 0.201144 and 
n = -1.391611. Our fit in Fig. [Elis given hy A = -0.391611 and n = -1.447867. 

5.2.2 Cross Section of Black Hole Interaction 

In its simplest form, calculating the cross section for the interaction of neutrino + 
nucleon —>■ Black Hole involves using the Thorne Hoop Conjecture [11 7j and checking 
to see if the neutrino and parton come close enough together to be within the radius 
of the Schwarzschild black hole that would be formed from their combined center-of- 
mass energy. At this stage of our simple approximation of checking to see if the impact 
parameter, b, is smaller than the Schwarzschild radius, r^, for our cross section, we 
would simply have the area vrr^ of a disk. 

To find the center-of-mass energy, Eqm, we use the conservation of relativistic 
4-momentum, Ptotai" = Pu'' + Ppa.rton"' ■ We define our metric as rjab = diag(— 1, 1, 1, 1). 
In the lab frame of IceCube, P,," = {Ei,,pi,^a:,Pu,y,Pu,z) and Pparton" = (xmA?, 0, 0, 0), 
where the variable x is the fraction of the total rest-mass energy of the nucleon present 
in the parton at the time of the interaction. Squaring the 4-momentum, |118j 

Ptotal Ptotal a (Pparton ~l~ Pf )(Pparton a ~l~ Pua) 

— P -I-9P a p , p p a A\ 

^ parton ^ parton a ~r ^-i parton ^ i/a i va^ v i v"'^/ 

and using 

n2 



we have 



PaP"" = -E'cM, (5.5) 



P'CM = ~^ + ^{-xmNEi, + ■ pij) - m^. (5.6) 



Because we are interested in energies where E,y ^ and m,^, we find 
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(5.7) 



We denote this quantity by 



s = 2x171^ El,. 



(5. 



In terms of the variables V^, the neutrino-parton center-of-mass energy; D, the 
total number of dimensions of spacetime; and M^), the D- dimensional Planck scale; 
we express the Schwarzschild radius as |106l I107j 



,{VS,D,Mn) 



M, 



D 



D 



1 



2D-4^(I3-7)/2p(^D_l^ 

D-2 



1 

D-3 



(5.9) 



From here on we will work within the assumption of string theory that D = 10, and 
we will have M^) = Mio, which we will eventually take to be near 1 TeV [119j . In 10 
dimensions, we then have 



,(71, M: 



10 J 



10 



10 



Stt^/^ r(9/2) 



1/7 



(5.10) 



for the Schwarzschild radius. 

The actual radius of the black hole will differ from the Schwarzschild radius r^, 
due to factors such as angular momentum and the geometry of spacetime, and we 
will call this corrected cross sectional area -Fvrr^, where the variable F is a prefactor 
used to correct for differences from an exact Schwarzschild metric. We define the 
inelasticity as |120] 



M 



y 



BH 



(5.11) 



which is a measure of how much of the center-of-mass energy is available to the black 
hole for Hawking radiation [12H 11221 1123j . The energy difference, the deficit between 
the final mass of the black hole after its ring-down phase |124[ I125[ 11261 1127j and 
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the center-of-mass energy initially present in the collision, is carried off via incoming 
shock wave multipole moments radiating gravitational waves |1281I129] 11301 [T3T|ll32j . 
The inelasticity y depends on the impact parameter b, and we define 

^ (5.12) 



5 



where 6max = vFts. The values of F and y{z) calculated depend upon the slicing of 
spacetime used to determine whether or not an apparent horizon is present. In the 
work of |1331I134|I135] . it is found for D = 10 that F = 1.819 and we approximate their 
findings for the inelasticity as y{z) = 0.59 — 0.572;^. In the later work of |136j . in which 
a slicing on the future light cone is used, it is found for D = 10 that F = 3.09 and 
we approximate their findings for the inelasticity as y{z) = 0.59 — 0.59^;^ + 0.234^^. 
We will refer to these two different slicings as the "old slice" and the "new slice," 
respectively. 

The prefactor F and the inelasticity y{z) were derived using classical general 
relativity. Since we don't yet have a quantum theory of gravity, we need to make sure 
we stay within a semi-classical regime. We expect a thermal distribution of Hawking 
radiation [1371 USHl [139] for 

Mbh > XminMio, (5.13) 

where Xmin = 3 ensures a well-defined resonance not dominated by the 3-brane tension 
[120[ I140j . and thus Mbh > 3 TeV. The thermal distribution of Hawking radiation 
is a Planckian spectrum, where the emission rate per degree of particle freedom i of 
particles of spin s with initial total energy between uj and u + duj is 



duj {D-2){27i)^ 

where 
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is the instantaneous Hawking temperature, 

= r[(D-2)/2] ^^-'^^ 

is the volume of a unit — 3)-sphere, and as{uj) is the greybody factor that accounts 
for the backscattering of part of the outgoing radiation into the black hole [142] . Note 
that a rough estimate of the instantaneous Hawking temperature can be found from 
the first law of black hole thermodynamics (which is analogous to the combined first 
and second law of thermodynamics): T = dE/dS ~ {dA/ dM)~^ [107] . Combining 
Eqs. (ESI), (EUD, and fl513l) shows that 

(XminMio)^ /cr 1 

where to find the cross section we integrate the PDFs over the parton momentum 
fraction x and use x GUI' lower limit of integration. 

In addition to integrating the PDFs over the parton momentum fraction, we also 
integrate over z for an impact parameter-weighted average over parton cross sections. 
The area of a thin ring of inner radius z and thickness of dz is proportional to zdz. 
We multiply this by a factor of 2, so that when we integrate zdz alone, we get 
a factor of 1; therefore, if y{z) did not depend on z, this weighted average could be 
neglected. Because the value of y{z) does in fact depend on the weighted average 
ensures we use the correct lower limit of integration, x, when integrating over the 
parton momentum fraction, x. 

The final expression for the vN —>■ BH cross section is |143j 



(7= [ 2zdz [ dxFTTrl{y/I,Mio) y^M^^Q)^ 
Jo Jx 



(5.1^ 



where i labels parton species, and the fi{x,Q) are PDFs. 

Fig. 15.51 shows log(cr) plotted versus log(-E'i,) for both the case of apparent horizons 
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Figure 5.5: Cross Section: New and Old Slicing. 
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Figure 5.6: Cross Section: Varying Semi- Classical Regime. 



on the "old slice" (dot-dash line) and the "new slice" (solid line). The cross section 
is given in units of picobarns (pb) and the energy of the incoming neutrino is given 
in units of GeV. We use Xmin = 3 and Mio = 1 TeV. 

Fig. 15.61 shows log(cr) plotted versus \og{E^) on a log-log scale for different values 
of Xmin using the new slicing. The cross section is given in units of picobarns (pb) 
and the energy of the incoming neutrino is given in units of GeV, where Miq = 1 TeV 
and Q = min{r~^, 10 TeV}. 

The cross sections were integrated with a variable step size with respect to the 
parton momentum function x. The dominant contribution from the PDFs comes 
from the small-x region, which is only probed when the lower limit of integration x 
is sufficiently small. This happens with large enough values of the incoming neutrino 
energy E^,. We keep the step variable of integration smaller than x/100 for x < 10^^ 
and equal to 1/100 for x > 10^^. This gives us excellent accuracy and a fast numerical 
calculation of the cross section. 



5.3 Detection Sensitivity 

One of the major outstanding questions that IceCube is hoped to be able to answer 
is — what is the flux rate of cosmic neutrinos? A good estimate involves a consid- 
eration of the number of neutrinos expected to be created in association with the 
observed flux of charged cosmic ray particles: this is the Waxman-Bahcall (WB) flux 
[m] of 

(j)^ - 6.0 X 10-^{E^/GeV)-'^ GeV"^cm-^s"^sr-\ (5.19) 

including all species of neutrinos. Another estimated flux assumes that extragalactic 
cosmic rays dominate the spectrum at energies above ~ 10^'^ GeV and that additional 
neutrinos are to be expected from sources opaque to ultra-high energy cosmic rays; 



Table 5.1: Probability of Signal. 



Mbh 


^sig 


3Mio 


0.078203 


4 Mio 


0.122514 


5 Mio 


0.161455 


6Mio 


0.196967 


7 Mio 


0.230733 



this is the AARGHW flux [HS] of 



^ 3.5 X 10-=^(E^/GeV)-2-54 GeV 



""^cm ^s ^sr 



(5.20) 



including all species of neutrinos. 

To confirm the existence of black hole interactions amidst the background noise of 
standard model (SM) interactions, we pick out a signal that has a high likelihood for 
the relatively democratic Hawking radiation and a low likelihood for charge current 
(CC) interactions: we search for soft muons, or muons with less than 20% of the 
incident neutrino energy. In SM CC interactions, a produced muon will generally 
carry away at least 80% of the incident energy. We only consider interactions with 
at least 4 secondary particles, where at least one of them is a muon ^l08j . The cross 
section for the SM CC interaction producing a soft muon is |7j 

(T^>°-'~1.2 (^,/GeV)°-=^^%b . (5.21) 

For incident neutrino energies larger than 10'' GeV, the background number of SM 
CC interactions meeting these criteria for the AARGHW flux, which produces more 
events than the WB flux, is 10 events over the 15 year lifetime of IceCube. For 

> 10^ GeV, the expected event rate for the SM CC interaction over IceCube's 
lifetime is less than 1 event. 

The probability that a black hole interaction produces the criteria we propose 
to search for depends on the mass of the black hole formed [7]. See Table 15.11 for 



some values of the signal probability versus the size of the black hole created from 
a neutrino-parton interaction. In this probability we neglect the gravitons radiated 
into the bulk of the compactified dimensions, but these are thought to carry away 
less than 15% of the radiated energy when D = 10 [1461 114711148] . 

With the probability of signal given as a function of black hole mass, we need a 
way to determine the expected number of TeV-scale black holes formed within a given 
mass range. We do this by dividing the expected number of black holes produced into 
bins at 0.1 Miq mass intervals. We vary our value of Xmin, and repeat our calculation 
for the expected number of black holes created at IceCube. For example, should 
a rate of 235 TeV-scale black holes be created at IceCube for Xmin = 3.1, and 246 
created for Xmin = 3.0, then we could assign 11 black holes to the Mbh = 3.0 TeV 
bucket. Each of these eleven black holes has a probability of about 0.078 to produce 
our signal, so this means our expected detection rate for this bin is approximately 
0.86 TeV-scale black hole signals over the 15-year lifetime of the IceCube experiment. 
When we calculate the rates and associated signal probabilities for all of our buckets 
in bins of Xmin > 3, we get our cumulative totals. 

We will integrate, with respect to energy, the neutrino flux over the 15 year lifetime 
of the IceCube experiment, or T ~ 4.7 x 10^ seconds. At the energies of interest the 
Earth is opaque to neutrinos. Hence, we will only consider neutrinos passing down 
through the Antarctic ice, and we will only accept measurements from this half of the 
available directions, which makes for 27i steradians of solid angle for observation. The 
background rate of non-neutrino events at such high energies is entirely negligible. 
IceCube's effective volume is 1 km^ |104j . which at a density of 900 kg/m^ means the 
number of nucleons available for neutrino interaction targets is ut — 5.4 x 10'^'^. Our 
upper limit of integration is an energy of 10^° GeV, because beyond this the IceCube 
detector will be saturated and unable to resolve all the details of the interaction [99]. 
The total number of black hole signal events over the life of the IceCube experiment 
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Table 5.2: Number of Signal Events. 





A/'bh 


[WB] 


TVbh [AARGHW] 


3 


43 


(19) 


69 (30) 


4 


34 


(15) 


43 (19) 


5 


27 


(12) 


28 (12) 


6 


22 


(9) 


20 (9) 



Table 5.3: 10-Dimensional Planck Mass Sensitivity. 





Mio/TeV [WB] 


Mio/TeV [AARGHW] 


3 


1.5 (1.2) 


1.5 (1.2) 


5 


1.3 (1.1) 


1.3 (1.1) 


7 


1.2 (1.0) 


1.2 (1.0) 


9 


1.1 (1.0) 


1.1 (0.9) 



is 

Mig = 27rnTT j dE, a{E,) ^E,) V,,^ . (5.22) 

In Table 15.21 we calculate the expected number of black hole signals over the 
lifetime of IceCube. With a lower limit of integration of 10^ GeV, we fix Mio = 1 
TeV, but we allow Xmin to vary. We compare the number of events for the WB flux 
to the AARGHW flux. For each flux, we have calculated the number of events using 
both the "new slice," which is given without parentheses; and the "old slice," which 
is given inside parentheses. 

In Table [531 we calculate the maximum 10-dimensional Planck mass for which we 
would expect be able to observe the interaction at the 3(j level. With a lower limit 
of integration of 10*^ GeV, and for differing values of Xmin, we flnd the corresponding 
value of Mio. We do this for both the WB flux and the AARGHW flux. For each 
flux, we have calculated the number of events using both the "new slice," which is 
given without parentheses; and the "old slice," which is given inside parentheses. 

In Fig. 15.71 we plot the TeV-scale discovery reach for both IceCube and the Large 
Hadron Collider (LHC) |120j . assuming a cumulative integrated luminosity of 1 ab~^ 
over the life of the collider. We calculate the maximum value of Mio that could be 
observed at the 5cr level versus Xmin, and we use a lower limit on the energy integral 
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Figure 5.7: IceCube and LHC Discovery Reach. 
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of 10^ GeV. We plot the IceCube discovery reach only in the semi-classical regime of 
a^min > 3; however, the LHC could potentially be focused on superstring resonances 
[1491 1150[ 1151] , and could thus be able to probe the quantum regime [7j. 
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Chapter 6 
Conclusion 



Using the de Sitter Bunch-Davies state for modes of intermediary-g2 and large-g2 is 
valid in the exponentially growing region of the scale factor, but imposing the Bunch- 
Davies state on modes of small-g2 leads to infrared divergences in the dispersion 
spectrum. Maintaining continuity of the scale factor to is necessary to prevent 
ultraviolet divergences of the energy density of particles created during inflation. The 
asymptotically Minkowskian regions of our composite scale factor do not affect the 
near scale-invariance of the intermediate-g2 region of particle production, but it does 
allow for an unambiguous interpretation of the number of particles produced versus 
^2, and it allows for flat-space renormalization. An asymptotically flat scale factor 
segment may be joined continuously to with an exponentially growing segment 
of scale factor, whereas a simple power law such as a{t) oc may not. Both of our 
massive approximations are trustworthy in their respective regimes: little growth of 
the composite scale factor outside of the exponentially growing region for the effective- 
k approach, and with modes not at the interface between the small- and intermediate- 
q2 behavior and not at the interface between the intermediate- and large-q'2 behavior 
for the dominant-term approach. In our model, the average number of particles 
created per mode can be characterized in terms of three parameters: the number of 
e-folds, A^e; the ratio of the mass to the Hubble constant during inflation, ttih', and 



the dimensionless mode number, 52- We find a scale-invariant spectrum when i^infl 
and niH are both constant, provided modes are converted individually into curvature 
perturbations soon after exiting the Hubble radius. The spectral index can be shifted 
towards a blue spectrum if all the curvature perturbations are created around the 
same time or at a time after the end of inflation. The spectral index can be shifted 
towards a red spectrum by taking into consideration a changing value of i/infl or 0. 
We find that an abrupt end to inflation leads to a boosted production of high-energy 
particles and an associated high temperature. If monopoles, or certain other exotic 
particles, were found to be created copiously at low temperatures — at the LHC, for 
instance — it could place rigorous constraints on the characteristics of inflation. 

The predicted energy and angular momentum in the post-Minkowski approxi- 
mation for our binary point mass system with helical symmetry agrees to first post- 
Newtonian order in the case of parametrization-invariant action plus either of the IPN 
correction terms. With g — > 0, we can make a comparison with the Schwarzschild 
solution of General Relativity. Here, both the affine case and the parametrization- 
invariant with a IPN correction term have an Innermost Circular Orbit at about 3M, 
which is outside the event horizon of GR located at 2M. Only the affine case has an 
Innermost Stable Circular Orbit, and it occurs at ~ 7.69 M, which is outside of the 
ISCO predicted by GR located at 6M. These discrepancies may be due to the lin- 
ear order of the post-Minkowski approximation, or they may be due to the radiation 
being pumped into the binary system by the half-advanced plus half-retarded helical 
symmetry. A form of the first law of thermodynamics dE/dv = VLdL/dv is satisfied, 
and this serves as a useful check on the analytical and numerical results. 

With a flux of cosmic neutrinos at the Waxman-Bahcall rate, over its 15 year 
lifespan the IceCube Neutrino Telescope could detect TeV-scale black holes at the 
5(7 level up to a maximum 10-dimensional Planck mass of 1.3 TeV. Our analysis 
shows that PDFs can be approximated well by fits to xf{x) — Ax^, provided the 
range of the parton momentum fraction, a;, for each fit is restricted to a few decades 



of variation on a log^g scale. The fitting of the parameters A and n can best be 
accomphshed by simultaneously varying each, and by sampling a large enough area 
of parameter space to ensure a false minimum deviation is avoided. The integration 
involved in calculating the cross section of the gravitational interaction between a 
parton and a cosmic neutrino is most efficiently carried out with a variable step size 
of integration. Values of the parton momentum fraction closest to the lower limit of 
integration dominate the cross section, so care must be taken to use a small enough 
step size in this range so that these values are not over-weighted in the cross section. 
A convenient way of associating events with a given value of Mbh is to recalculate the 
number of lifetime events for different values of Xmin, and then subtract the difference 
between the events from incremental values of Xmin into bins. 

In the three parts of this dissertation, we have focused on the topics in inflationary 
cosmology and astrophysics described in three papers: O [6l [7] . 
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Appendix A 

Composite History of an Exact 
Reaction-Force Solution 

This Appendix is motivated by and based on the work of [152J. What follows is 
an application of the more general techniques presented in Sec. 13.1.21 for matching 
continuously to second derivative in what could be taken as either a scale factor on 
the one hand or as a particle's velocity on the other. It is hoped that this example 
serves to illustrate some aspects of the self force and radiation reaction mentioned in 
Chapter |H We begin with a charged particle that in its rest frame emits radiation 
when accelerated as given by the Larmor formula (in Gaussian units) of 

P = ~v\ (A.l) 

which leads to, in addition to the external force, a radiation-reaction force of the form 
[1531 [IMl [155] 

• 2 ■■ 

^applied = mv- -— (A. 2) 



165 

as perceived by the particle in its momentarily-comoving rest frame. In this example, 
we will consider only rectilinear motion, so we rewrite this as 



-P'applied = mV - TTITV, 



(A.3) 



where 



2 e2 



3 mc^ 



We thus define the reaction force, or self force, as 



(A.4) 



^seif = -rriTV. 



(A.5) 



For constant acceleration, we have 



Vcit) 
Vcit) 
Vc{t) 
Fc seu(t) 
Fc applied (^) 
Pc radiated (^) 



ac{t - to) 

0, 
0, 

mac, 
mra^ . 



(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.IO) 
(A.ll) 



We then introduce a two similar velocity histories given by a hyperbolic tangent in 
analog with Sec. 13.1.21 



V^{t) 



voi + Ai tanh ■ 



Vf{t) = t>o/ + A/tanh 



t-tf 



(A.12) 
(A.13) 



where A, is twice the difference between early- and late-time velocities for the first 
velocity history, ti is the time at which Vi{t) = Voi, and Si is a throttling parameter 
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that decreases the change in velocity with respect to time as it increases in magnitude; 
and where A/ is twice the difference between early- and late-time velocities for the 
final velocity history, t/ is the time at which Vf{t) = vqj, and s/ is a throttling 
parameter that decreases the change in velocity with respect to time as it increases 
in magnitude. We will take both Si and s/ to be > 0. Then we have 



and 



Vi{t) 
Vi{t) 
Vi(t) 
Fi self(^) 
Fi applied (^) 



i radiated 



Voi + Ajtanh , 

Si 

^fl-tanh^^ 

Si V s,: 

Ai f t-t 
2—^ tanh 



^)[( 



0.2 \ 



Ai / , t-ti 
— zmr — ^ I tanh 

Si \ Si 

A. / , .2^- U 
m — I 1 — tanh 

Si \ Si 



tanh" 



Si 



tanh' 



2mr^(tanh 



t-ti 



t - u 



it) = 



mr 



2 t ti 



— ( 1 - tanh' 

Si \ Si 



tanh 
V Si 

2 



(A. 14) 
(A.15) 
(A.16) 
(A.17) 



(A.18) 
(A.19) 



7 self 



/ applied 



it) 



it) 



f radiated 



it) 



t-tf 

^fl-tanh^^ 

Af I ,t — tf\ ( , ^2t~tf 



Vf{t) — I'oz + Aj^tanh 

Vfit) 
Vfit) 



Af / 
2^ tanh ■ 



tanh 



^ Af ( , t-tf 
—2mT — - tanh 



)[( 



tanh^ 



t-tt 



t - t 



m— ( 1 - tanh^ 

Af / 
—2mT — - tanh 



t-tf 



tanh 



2t-tf 



mr 



^fl-tanh^^-^^ 



- 1 



(A.20) 
(A.21) 
(A.22) 
(A.23) 



(A.24) 
(A.25) 



At times tj and tf, respectively, we have 



and 







(A.26) 


Vi{U) 


A, 

Si 


(A.27) 


Vi{ti) 


= 0, 


(A.28) 


Pi self(^j) 


= 0, 


(A.29) 


Pi applied (^j) 


A, 

= m — , 

Si 


(A.30) 


Pi radiated (^) 


/A \ ^ 


(A.31) 






(A.32) 


Vfitf) 




(A.33) 




= 0, 


(A.34) 


Pf self (tf) 


= 0, 


(A.35) 


Pf applied (^/) 


= m — 


(A.36) 


Pf radiated (^) 


= rriT — - 
\Sf J 


(A.37) 



We then specify a composite velocity history by matching the velocity histories of 
Vi{t) to Vc{t) to Vf{t). We can maintain joining conditions — meaning the velocity, 
acceleration, and radiation-reaction force are all kept continuous — by joining the 
initial segment to the start of a region of constant acceleration at t = tj, and by 
joining the final segment to the end of a region of constant acceleration aX t = tj. 

See Fig. lA.ll where we plot a dimensionless example of a composite velocity where 
Aj = Af = Si = Sf = Qc = I- In this example we take tt = 0, tj = 10, Voi = 1, and 
Vof = 11. 

Maintaining continuity of the velocity history up to its second derivative imposes. 
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in addition to the two conditions of matching times, the following boundary conditions 

A, 



Si 



voi = lim [v(t)] + Ai, 

t— >— oo 



vof = lim [v{t)]-Af. 

t — >- + CxD 



(A.38) 

(A.39) 
(A.40) 



We find that = ti — VQi/a^ and the duration of constant acceleration is = tf — ti. 
See Fig. IA.2t where we plot a dimensionless example of the applied force necessary 




■0.2 - 







5 
t 



10 



Figure A. 2: External Force versus Time. 



15 



to maintain the motion of the particle shown in Fig. lA.ll for two different dimensionless 
values of r. In the small-r limit, we get a Newtonian 2nd Law of F = ma. In the 
large-r limit, we note some peculiarities of the self-force. To initiate the acceleration, 
a force must initially be applied opposite to the direction of motion — this is to be 
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compared with the pre-acceleration found for radiation-reaction forces that ehminates 
runaway-acceleration solutions. To end the period of constant acceleration, the force 
must be increased in the direction of motion. As will be shown below, this additional 
work is needed to compensate for the energy dissipated by the radiation emitted. 
The total change in kinetic energy of the particle is 



AKE = <( ( lim [v{t)] ) - ( lim [v{t)] 



-m 



{{voi + acta + Aff - {voi - Aif] 



m Qc Voi {ta + Si + Sf) + -mac^ {ta'^ + 2taS f + s/ - Si^) . ( A.41) 



The total power radiated is 



radiated total 



-OO / \^ ti 

2 Ai^ 2 2 Ai^ 
-mr h mrCc ta + -ttit 

O Si O Si 



mrac' [ -Si + ta + 



Pf dt 



(A.42) 



The total work done on the particle by the external force is 



total 



(/: 



applied 



Vi{t) dt 



c applied 



Vc{t) dt 



+ 00 



+ I ^ Ff e^pplied Vf{t) dt 



— ma 



I 2 1 

Voi{Si - t) + SjOc <i -r - -s, 



m ac Voi ta + ta^ 



+ mar 



VQi{Sf + r)+ac[^s/ + Sfta +'^Sfr + taT 



= macVQi {si + ta + Sf) + -mac^ {ta^ + 2taSf + s/^ - Si^) 



^mrac^ \ -Si + ta + g^/ 



(A.43) 
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We find that 



W^total - Eradiated - E = 0, 



(A.44) 



and thus energy is conserved at early and late times. See Fig. IA.31 for the case of 
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Figure A. 3: Energy versus Time. 



energy conservation between early and late times. The velocity history is given in 
Fig. lA.ll and we choose r = 1. The energy deficit that develops is primarily due 
to the energy dissipated through the emitted radiation during the phase of constant 
acceleration. This negative energy must be balanced by an additional amount of work 
applied to the particle to end the acceleration. If additional energy is not provided to 
the system, Wiseman has proven that the kinetic energy of the particle decreases to 
compensate |152j . In the limit of Sj — and Sj — 0, we see that the work associated 
with overcoming the reaction force at the initial and final joining points is —macVoi t 
and macVofT, respectively. Because in this velocity history v = if t ^ ti and 
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t ^ tf, and because H^eif = / Fseiiify'^if) dt, we see that in the instantaneous hmit, 
-^seif(^) = mOcT [5{t — tf) — S{t — ti)], where S{t) is the Dirac delta-function. 
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